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In  this  dissertation,  peak-to-average  power  ratio  (PAPR)  suppression  schemes  for 
digital  modulations  are  investigated  and  developed.  The  PAPR  of  a  signal  is  the  ratio  of 
the  peak  signal  power  over  the  average  signal  power  and  a  high  PAPR  signal 
substantially  degrades  the  efficiency  of  the  linear  power  amplifier. 

In  portable  wireless  communications,  bandwidth  efficiency  becomes  more 
important  than  power  efficiency,  thus;  a  bandwidth  efficient  linear  amplifier  and  7i/4 
QPSK  have  become  popular.  Several  PAPR  suppression  schemes  were  suggested  to 
increase  the  power  efficiency  in  t:/4  QPSK  but  they  degraded  the  bandwidth  efficiency. 
Therefore,  an  adaptive  peak  suppression  (APS)  algorithm  is  introduced  to  suppress  the 
PAPR  without  any  bandwidth  efficiency  degradation.  It  adjusts  the  amplitudes  of 


symbols  before  pulse  shaping.  Then,  the  PAPR  of  n/4  QPSK  with  root  raised  cosine 
(RRC)  (roll-off  0.2)  pulse  shape  can  be  suppressed  from  4.8dB  to  1.9dB. 

Orthogonal  frequency  division  multiplexing  (OFDM)  is  an  effective  and  popular 
modulation  for  high  data  rate  transmission  but  causes  very  high  PAPR.  Most  of  the 
previous  PAPR  suppression  schemes  have  weak  performance  or  expand  the  bandwidth. 
Therefore,  a  limiter  plus  APS  scheme,  which  is  the  combination  of  an  amplitude  limiter 
and  the  APS  algorithm,  is  developed.  The  scheme  suppresses  the  PAPR  of  128-QPSK- 
OFDM  with  RRC  (roll-off  0.35)  pulse  shapes  from  12.7dB  to  5.2dB  without  any 
bandwidth  efficiency  degradation  when  IdB  average  power  loss  is  allowed  in  an  AWGN 
channel. 

Synchronous  muti-carrier  code-division  multiple-access  (MC-CDMA)  with  Golay 
complementary  spreading  sequence  has  similar  PAPR  properties  as  OFDM.  Thus,  the 
limiter  plus  APS  scheme  is  also  applied  to  synchronous  MC-CDMA  and  suppresses  the 
PAPR  of  128-MC-CDMA  with  RRC  (roll-off  0.35)  pulses  from  12.1dB  to  3.9dB  when 
IdB  average  power  loss  is  allowed.  The  PAPR  of  asynchronous  MC-CDMA  is 
determined  by  its  spreading  sequence.  Zadoff-Chu  sequences  are  currently  known  as  the 
best  spreading  sequence  for  low  PAPR.  However,  a  new  proposed  PAPR  suppressed 
Zadoff-Chu  (PSZ)  sequence  causes  2-3dB  lower  PAPR  than  Zadoff-Chu  sequences 
without  any  cost. 

The  performances  of  the  PAPR  suppression  schemes  of  this  dissertation  are  very 
competitive  when  compared  to  previous  schemes;  moreover,  these  schemes  do  not 
sacrifice  any  bandwidth  efficiency. 


CHAPTER  1 
INTRODUCTION 

In  this  chapter,  the  concept  of  bandwidth  and  power  efficiency  of  a  modulation 
format  will  be  presented.  The  effect  of  Peak-to-Average  Power  Ratio  (PAPR)  on  the 
power  efficiency  in  a  wireless  communication  system  will  be  described.  This  chapter  will 
also  summarize  the  previous  work  on  finding  a  bandwidth  and  power  efficient 
modulation  and  introduce  three  popular  modulation  schemes  that  require  PAPR 
suppression.  An  overview  of  the  dissertation  will  conclude  this  introductory  chapter. 

1 . 1  Bandwidth  and  Power  Efficient  Modulation 
Throughout  the  history  of  wire  or  wireless  communications,  there  has  been  a  lot 
of  search  for  better  modulation  schemes.  Such  an  improvement  means  better  bandwidth 
efficiency,  better  power  efficiency,  more  robustness  to  noise  and  fading,  improved  cost 
efficiency,  easier  of  implementation,  and  so  on. 

Bandwidth  efficiency  is  one  of  the  most  important  factors  among  them.  The 
bandwidth  is  the  difference  between  the  highest  and  lowest  frequencies  of  a  data 
transmission  channel  and  is  proportional  to  the  capacity  of  the  transmission  channel. 
Therefore,  if  a  larger  bandwidth  is  provided,  a  larger  amount  of  data  can  be  transmitted  or 
received  in  a  fixed  time  interval.  As  an  information  society  develops,  more  and  more  data 
need  to  be  transmitted  to  provide  high  fidelity  and  larger  user  capacity.  However,  there  is 
a  finite  bandwidth  in  any  wire  or  wireless  channel.  Therefore,  the  bandwidth  availability 
becomes  severely  limited  and  a  more  bandwidth  efficient  modulation  technique  becomes 
necessary.  The  bandwidth  efficiency  is  defined  as  the  amount  of  transmitting  data  for  a 
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given  time  and  for  a  given  bandwidth,  so  the  general  unit  of  the  bandwidth  efficiency  is 
bits/sec/Hz. 

In  a  portable  wireless  communication  system  such  as  a  cellular  phone  or  a 
portable  radio,  small  size  and  long  operation  time  are  desirable.  The  size  of  a  portable 
unit  is  determined  by  the  battery  size  because  the  battery  is  often  bigger  than  the  rest  of 
the  unit.  Therefore,  the  battery  size  is  limited  in  order  to  make  the  unit  portable.  To 
achieve  long  operation  time  for  a  limited  battery  size,  a  power  efficient  system  is 
required.  The  power  efficient  system  can  be  provided  by  using  a  power  efficient 
modulation  technique. 

There  are  two  definitions  of  power  efficiency  in  digital  modulations.  One  is  how 
much  average  power  is  required  to  support  a  given  quality  of  service,  which  is  measured 
in  terms  of  bit  error  rate  (BER).  The  other  is  called  a  temporal  efficiency,  which  is 
determined  by  the  width  of  the  amplitude  variations  of  a  transmitted  signal,  and  can  be 
measured  by  peak-to-average  power  ratio  (PAPR).  This  PAPR  determines  the  efficiency 
of  a  power  amplifier,  which  consumes  the  major  power  in  a  portable  communication 
system.  Two  modulation  formats  may  have  different  temporal  efficiencies  (PAPR)  while 
having  the  same  average  power  efficiency  (BER).  The  power  efficiency  of  a  system  is 
defined  as  the  ratio  of  transmitted  signal  power  to  the  power  consumed  in  the  system. 
This  power  efficiency  of  a  system  is  affected  by  both  the  average  power  efficiency  of  the 
modulation  format  (BER)  and  the  efficiency  of  the  power  amplifier  that  is  determined  by 
the  temporal  efficiency  (PAPR). 

Generally,  bandwidth  efficient  modulation  formats  are  not  power  efficient  and 
power  efficient  modulation  formats  are  not  bandwidth  efficient.  For  example,  M-ary 
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Figure  1.1:  Average  power  efficiency  vs.  spectral  efficiency  of  each  modulation 


phase-shift  keying  (MPSK)  or  M-ary  quadrature  amplitude  modulation  (MQAM)  is  more 
bandwidth  efficient  and  less  power  efficient  as  M  (>2)  increases,  and  M-ary  frequency- 
Shift  Keying  (MFSK)  is  more  power  efficient  and  less  bandwidth  efficient  as  M  (>2) 
increases,  as  can  be  seen  in  Figure  1.1.  The  spectral  efficiency  is  another  name  for 
bandwidth  efficiency.  The  average  power  efficiency  can  be  calculated  by  using  the 

relative  power  efficiency  formula  [rj^  -  \  0\og^Q(d^^^J  /4£'^)]  where  dmin  is  the  minimum 

distance  between  points  in  the  signal  constellation  and  Eb  is  energy  per  bit.  The  relative 
power  efficiency  is  comparing  the  required  signal  power  to  achieve  a  fixed  BER  in  an 
additive  white  Gaussian  noise  (AWGN)  channel.  The  spectral  efficiency  is  measured  in 
terms  of  null-to-null  bandwidth  in  Figure  1.1  and  Figure  1.3.  In  those  figures,  the 
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bandwidth  of  MQAM  or  MPSK  is  measured  when  the  Nyquist  root  raised  cosine  (RRC) 
pulse  with  roll-off  factor  a=.35  is  adopted.  Generally,  quadrature  phase-shift  keying 
(QPSK)  and  16QAM  represent  a  good  trade-off  among  bandwidth  efficiency  and  average 
power  efficiency. 

1 .2  The  Relation  between  Power  Amplifiers  and  Modulation 
In  portable  wireless  communications,  long  operation  time  and  small  unit  size  are 
required  and  these  can  be  achieved  by  a  power  efficient  system.  The  power  efficiency  of 
a  system  is  decided  by  the  efficiency  of  the  power  amplifier  since  the  power  amplifier  is 
the  major  source  of  the  power  consumption. 

There  are  two  kinds  of  power  amplifiers  in  wireless  communication  systems.  One 
is  a  linear  amplifier  and  the  other  is  a  nonlinear  amplifier.  A  linear  amplifier,  such  as  a 
class  A,  a  class  AB  or  a  class  B  amplifier,  is  less  power  efficient  than  a  nonlinear  one. 
However,  a  linear  amplifier  is  necessary  for  any  modulation  format  with  an  envelope 
which  is  not  constant.  Otherwise  the  signal  will  be  distorted,  destroying  its  bandwidth 
efficiency.  The  efficiency  of  a  linear  amplifier  is  highly  dependent  on  the  PAPR  of  a 
signal,  as  can  be  seen  in  Figure  1.2  [1].  Therefore,  a  low  PAPR  signal  is  required  for 
improving  the  power  efficiency  of  a  linear  power  amplifier.  On  the  other  hand,  a 
nonlinear  amplifier,  such  as  a  class  C  amplifier,  is  much  more  power  efficient  than  a 
linear  one.  However,  the  nonlinear  amplifier  distorts  the  non-constant  envelope  signal  of 
a  linear  modulation.  The  distorted  signal  causes  sidelobe  regeneration  in  its  spectrum  and 
impairs  the  bandwidth  efficiency  that  is  the  major  merit  of  a  linear  modulation. 
Therefore,  a  nonlinear  amplifier  is  usually  used  for  a  constant  envelope  modulation 
which  has  no  amplitude  fluctuation,  while  a  linear  amplifier  is  usually  used  for  a  linear 
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Figure  1 .2:  The  theoretical  efficiency  of  Hnear  power  amplifier  [1] 


modulation  which  has  some  amplitude  fluctuations. 

As  mentioned  above,  there  are  also  two  kinds  of  modulation  formats,  a  linear 
modulation  and  a  constant  envelope  modulation.  Linear  modulation  schemes,  such  as 
phase-shift  keying  (PSK)  or  quadrature  amplitude  modulation  (QAM),  are  generally 
spectrally  efficient;  that  is,  they  have  high  bits/sec/Hz,  which  can  be  seen  in  Figure  1.3. 
However,  the  signals  of  these  formats  typically  have  high  PAPR,  which  reduces  the 
efficiency  of  a  linear  power  amplifier.  PSK  is  a  constant  envelope  modulation  when  a 
rectangular  pulse  is  used.  However,  the  signal  of  PSK  has  an  amplitude  fluctuation  when 
a  non-constant  overlapped  pulse  shape  such  as  a  time  domain  raised  cosine  (TRC)  pulse 
or  a  Nyquist  RRC  pulse  is  used.  On  the  other  hand,  a  constant  envelope  modulation,  such 
as  frequency-shift  keying  (FSK)  or  minimum-shift  keying  (MSK),  is  usually  spectrally 
less  efficient  than  linear  modulation.  However,  the  signal  of  the  constant  envelope 
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Figure  1 .3 :  PAPR  vs.  spectral  efficiency  of  each  modulation 


Modulation  format  has  no  PAPR  (OdB  PAPR),  i.e.,  it  is  more  power  efficient  than  that  of 
linear  modulation.  2FSK,  4FSK  and  8FSK  in  Figure  1.3  are  the  examples  of  constant 
envelope  modulation  formats.  A  good  compromise  for  both  temporal  efficiency  (PAPR) 
and  spectral  efficiency  is  n/4  shifted  QPSK  {n/4  QPSK)  or  offset  QPSK  (OQPSK). 

1 .3  Previous  Work  on  Bandwidth  and  Power  Efficient  Modulation 
Research  on  bandwidth  and  power  efficient  modulation  has  been  performed  since 
the  inception  of  digital  modulations.  In  the  early  days,  power  efficiency  was  a  bigger 
issue  than  bandwidth  efficiency.  So  constant  envelope  modulations  and  nonlinear 
amplifiers  were  usually  used.  In  papers  written  in  the  1970"s  such  as  Foschini  et  al.  [2]  or 
Simon  and  Smith  [3],  peak  power  efficient  linear  modulation  schemes,  which  have  low 
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PAPR,  were  studied.  However,  bandwidth  was  not  considered  and  the  resuhs  were 
presented  assuming  square  pulse  shapes. 

In  papers  written  in  the  early  to  mid  1980's,  a  nonlinear  amplifier  was  a  default 
power  amplifier  not  only  for  a  constant  envelope  modulation  but  also  for  a  linear 
modulation  because  of  high  power  efficiency.  For  example,  in  the  paper  written  by 
Morals  and  Feher  [4],  a  nonlinear  amplitude  limiting  amplifier,  such  as  a  traveling  wave 
tube  (TWT)  or  an  injection  locked  amplifier  (ILA),  was  used  not  only  for  constant 
envelope  MSK  but  also  for  non-constant  envelope  QPSK  and  OQPSK.  A  QPSK  or 
OQPSK  signal  with  a  rectangular  pulse  was  pre-filtered  by  a  Chevyshev  low  pass  filter 
(LPF)  and  amplified  by  a  hard  limiting  amplifier.  Then,  the  signal  envelope  becomes 
non-constant  after  pre-filtering,  but  becomes  almost  constant  again  after  the  hard  limiting 
amplifier.  So  the  sidelobe  is  reduced  after  pre-filtering  but  regenerated  after  a  hard 
limiting  amplifier.  However,  the  sidelobe  of  the  filtered  and  limited  signal  is  lower  than 
that  of  the  previous  rectangular  pulsed  signal,  so  there  is  some  bandwidth  benefit  from 
the  filtering  and  limiting.  According  to  the  paper  of  Morals  and  Feher  [4],  OQPSK  (linear 
modulation)  is  more  bandwidth  efficient  than  MSK  (constant  envelope  modulation)  in  a 
nonlinear  channel  environment  in  spite  of  the  sidelobe  regeneration  of  OQPSK.  The 
spectrum  of  QPSK  signals  and  OQPSK  signals  becomes  different  in  that  paper  by  the 
effect  of  a  nonlinear  amplifier. 

In  the  paper  written  by  Chang  and  Austin  [5],  a  Nyquist  raised  cosine  (RC)  pulse 
filter  was  adopted  as  a  pre-filter  for  better  bandwidth  efficiency  and  a  TWT  was  used  as  a 
nonlinear  amplifier.  They  showed  that  RC  pre-filtered  QPSK  is  more  power  efficient  than 
RC  pre-filtered  OQPSK.  In  another  paper  by  Chang  and  Ausfin  [6],  a  TRC  pulse  filter 
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was  used  as  a  pre-filter  and  a  TWT  was  also  used  as  a  nonlinear  amplifier.  They 
designated  the  TRC  pulse  shaped  (pre-filtered)  QPSK  as  quadrature  overlapped  raised 
cosine  (QORC)  and  the  TRC  pulse  shaped  staggered  QPSK  (SQPSK)  as  staggered 
QORC  (SQORC).  SQPSK  is  the  other  name  for  OQPSK.  They  showed  the  QORC  is 
more  power  efficient  and  SQORC  is  more  bandwidth  efficient  than  others  (QPSK, 
OQPSK,  MSK)  in  a  nonlinear  environment.  In  a  similar  paper  by  Le-Ngoc,  Feher  and 
Van  [7],  they  called  QORC  as  inter-symbol  interference  and  jitter-free  QPSK  (IJF- 
QPSK)  and  SQORC  as  IJF-OQPSK.  Seo  and  Feher  introduced  superposed  QAM 
(SQAM)  in  their  paper  [8].  SQAM  is  an  OQPSK  with  a  modified  TRC  pulse,  which  is 
more  bandwidth  efficient  than  a  TRC  pulse.  The  bandwidth  efficiency  of  each  pulse 
shape  can  be  compared  in  Figure  1 .4.  According  to  the  paper,  SQAM  is  more  spectrally 
efficient  than  OQPSK,  QBL  [9]  and  MSK  and  it  is  more  power  efficient  than  MSK, 
SQORC  (IJF-OQPSK)  and  TFM  [10]  in  a  nonlinear  hard  limiting  channel.  A  SQAM 
signal  is  very  bandwidth  efficient  in  a  nonlinear  hard  limiting  channel  because  of  very 
low  sidelobes  as  in  Figure  1 .4.  However,  the  mainlobe  of  a  modified  TRC  pulse  is  almost 
twice  the  size  as  that  of  a  RRC  pulse.  So  it  is  not  so  bandwidth  efficient  if  a  linear 
channel  is  provided,  i.e.,  when  a  linear  amplifier  is  adopted. 

In  papers  written  in  the  late  1980's,  researchers  began  to  consider  linear 
amplifiers  for  better  bandwidth  efficiency.  Akaiwa  and  Nagata  used  k/4  QPSK  with  a 
class  AB  linear  amplifier  in  their  paper  [11].  QPSK  was  excluded  in  their  research 
because  the  dynamic  range  of  the  signal  envelope  is  infinite  and  its  PAPR  is  high.  The 
dynamic  range  is  defined  as  the  maximum  signal  power  over  the  minimum  signal  power 
measured  in  decibels  (dB).  The  infinite  dynamic  range  and  the  high  PAPR  degrades  the 
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PSD  of  RRC(.20)  length  12Ts  PSD  of  RRC(.35)  length  8Ts 


Figure  1.4:  Power  spectral  density  (PSD)  of  each  pulse  shape 
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power  efficiency  of  a  system.  Generally,  the  irreducible  BER  is  higher  when  coherent 
detection  is  used  in  a  fast  fading  channel  because  of  random  frequency  modulation  (FM) 
noise.  In  that  case,  differential  detection  is  usually  preferred;  however,  the  differential 
detection  of  OQPSK  is  inferior  to  that  of  QPSK.  Therefore,  OQPSK  was  also  excluded  in 
their  research  because  it  is  not  robust  against  a  fast  fading  channel  even  though  it  has  a 
lower  PAPR  than  QPSK.  Thus,  n/4  QPSK  was  selected  because  it  has  lower  PAPR  than 
QPSK  and  can  be  demodulated  efficiently  by  a  differential  detector  in  a  fast  fading 
channel.  Liu  and  Feher  [12]  also  chose  n/4  QPSK  for  the  same  reasons  for  a  cellular 
digital  mobile  system  in  their  paper.  They  used  a  differential  detector  and  a  linear  power 
amplifier  in  a  co-channel  interference  (CCI)  and  AWGN  environment.  Yamao  et  al.  [13] 
compared  the  performance  of  QPSK  and  n/4  QPSK  for  a  newly  developed  high 
efficiency  linear  power  amplifier  for  digital  mobile  radio.  They  showed  that  n/4  QPSK 
was  slightly  better  than  QPSK  in  a  linear  channel. 

Recently,  Miller  and  O'Dea  [14]  introduced  an  adaptive  peak  suppression  (APS) 
algorithm,  which  suppresses  the  PAPR  of  PSK.  There  were  some  trials  conducted  in  an 
attempt  to  reduce  the  PAPR  by  adopting  an  LPF  when  there  is  a  peak  in  the  signal. 
However,  the  filtering  method  causes  bandwidth  efficiency  loss  because  it  regenerates  the 
sidelobe.  The  adaptive  peak  suppression  (APS)  algorithm  suppresses  the  peak  signal  by 
multiplying  an  amplitude  adjustment  factor  to  each  data  symbol  before  pulse  shaping.  So 
the  pulse  shape  is  not  changed  and  there  is  no  sidelobe  regeneration.  Therefore,  the  APS 
algorithm  increases  the  temporal  efficiency  of  PSK  without  any  bandwidth  efficiency 
loss.  Thus,  the  PAPR  suppressed  PSK  by  the  APS  algorithm  becomes  a  more  power 
efficient  modulation  than  the  original  PSK. 
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1 .4  PAPR  Suppression  in  Popular  Modulation  Schemes 
The  power  efficiency  of  a  modulation  can  be  improved  by  PAPR  suppression 
because  the  PAPR  suppression  improves  the  efficiency  of  a  linear  power  amplifier. 
General  PAPR  suppression  methods  lose  bandwidth  efficiency,  while  they  gain  power 
efficiency.  In  this  dissertation,  bandwidth  efficient  PAPR  suppression  schemes,  which  do 
not  lose  any  spectral  efficiency,  are  developed  and  applied  in  the  following  popular 
modulation  schemes. 

As  mentioned  in  the  last  section,  in  the  late  1980's,  researchers  began  to  adopt  a 
linear  amplifier  for  better  bandwidth  efficiency  because  the  bandwidth  efficiency  became 
more  important  as  the  amount  of  transmitting  data  increased,  n/4  QPSK  was  one  of  the 
best  modulations  for  the  linear  amplifier  in  a  portable  wireless  communication  system. 
The  reason  is  that  n/4  QPSK  is  a  bandwidth  efficient  linear  modulation  and  it  has  good 
average  power  efficiency  and  good  temporal  efficiency.  The  PAPR  of  n/4  QPSK  is  about 
3-5dB  when  a  Nyquist  RRC  (roll-off  a=0.2-0.35)  pulse  is  adopted.  If  the  PAPR  is 
suppressed  by  a  PAPR  suppression  scheme,  the  efficiency  of  the  power  amplifier  will  be 
increased  according  to  Figure  1 .2. 

Orthogonal  frequency  division  multiplexing  (OFDM)  has  been  found  to  be  an 
effective  data  transmission  scheme  in  a  multipath  fading  channel.  It  is  currently  one  of 
the  most  popular  modulation  schemes  for  the  transmission  of  high  data  rates  such  as 
multimedia  communications.  One  of  the  serious  problems  of  OFDM  is  very  high  PAPR, 
which  substantially  degrades  the  efficiency  of  a  power  amplifier.  The  PAPR  of  OFDM  is 
more  than  lOdB  in  most  of  cases.  Therefore,  PAPR  suppression  in  OFDM  is 
indispensable  for  the  power  efficiency  of  an  OFDM  system. 


12 

Multi-carrier  code  division  multiple  access  (MC-CDMA)  is  a  combination  of 
OFDM  and  direct  sequence  code  division  multiple  access  (DS-CDMA).  It  has  the 
advantages  of  both  schemes,  i.e.,  high  data  rate  and  high  multiple  access  capacity.  Thus, 
it  becomes  a  popular  multiple  access  technique  for  muhimedia  services  in  portable  radio 
communications.  Unfortunately,  the  PAPR  of  MC-CDMA  is  as  high  as  that  of  OFDM 
because  the  basis  of  MC-CDMA  is  the  same  with  OFDM.  Therefore,  PAPR  suppression 
in  MC-CDMA  is  also  required. 

These  three  modulation  schemes  are  currently  very  popular,  but  have  high  PAPR 
problems.  New  PAPR  suppression  schemes  will  be  developed  to  suppress  the  PAPR  of 
these  modulations  to  improve  the  power  efficiency  of  their  systems.  The  PAPR 
suppressed  modulations  will  be  new  bandwidth  and  power  efficient  modulations  that 
have  better  power  efficiency  than  the  original  modulations. 

1 .5  The  Outline  of  the  Dissertation 
This  dissertation  is  organized  as  follows.  Chapter  2  describes  the  APS  algorithm 

and  its  application  to  n/4  QPSK.  Section  2.1  introduces  n/4  QPSK  and  its  PAPR 
problem.  The  PAPR  problem  is  investigated  when  a  Nyquist  RRC  pulse  (roll-off  =  0.2, 
0.35)  is  adopted.  Section  2.2  introduces  the  APS  algorithm  which  suppresses  the  PAPR 
without  losing  any  bandwidth  efficiency.  However,  there  are  some  practical  problems  to 
implement  the  algorithm  for  n/4  QPSK.  Secfion  2.3  shows  the  previous  work  of  former 
researchers  on  the  practical  problems.  They  are  the  look-up  table,  artificial  neural 
networks  (ANN)  and  hybrid  methods.  The  consecutive  research  for  the  APS  algorithm 
and  the  look-up  table  is  conducted  and  shown  in  Section  2.4  and  Section  2.5.  In  Section 
2.4.1,  block  processing,  which  is  a  problem  of  the  APS  algorithm,  is  solved  by  using  a 
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fixed  long-term  normalization  factor.  The  major  complexity  reduction  of  APS  algorithm 
is  done  by  choosing  proper  parameters  as  in  Section  2.4.2.  Choosing  appropriate  number 
of  samples  also  reduces  the  complexity  as  in  Section  2.4.3.  In  Section  2.4.4,  the 
computational  complexity  is  investigated  in  detail  for  each  step  of  the  APS  algorithm. 
The  performance  of  the  algorithm  is  shown  in  Section  2.4.5  and  supplementary  schemes 
for  complexity  reduction  are  shown  in  Section  2.4.6.  The  major  problem  of  the  look-up 
table  is  that  it  requires  large  memory  for  n/4  QPSK.  However,  the  look-up  table  size  can 
be  reduced  by  factor  of  1/16  or  1/32  in  Section  2.5.1.  Section  2.5.2  explains  the  practical 
look-up  table  generation  and  Section  2.5.3  closes  the  chapter  comparing  the 
performances  of  the  new  look-up  table,  the  complexity  reduced  APS  algorithm  and  the 
new  hybrid  method. 

Chapter  3  investigates  PAPR  suppression  schemes  in  OFDM.  Section  3.1 
introduces  a  popular  modulation,  OFDM,  which  transmits  data  through  many  sub- 
channels. It  is  a  very  effective  modulation  format  but  has  a  high  PAPR  problem.  Section 
3.2  investigates  the  previous  schemes  on  OFDM  PAPR  suppression.  Most  of  the  schemes 
cause  bandwidth  efficiency  degradation  that  is  reluctant  to  be  sacrificed.  Thus,  an 
amplitude  limiter  is  suggested  as  a  bandwidth  efficient  PAPR  suppression  scheme  in 
Section  3.3  and  the  performance  is  shown  in  Section  3.4.  The  suppressed  PAPR  in  terms 
of  discrete  PAPR  is  increased  3-4dB  in  terms  of  continuous  PAPR  after  a  pulse  shaping 
filter  such  as  a  RRC.  Therefore,  the  APS  algorithm  is  applied  after  the  amplitude  limiter 
to  suppress  PAPR  caused  by  pulse  shaping  in  Section  3.5.  Then,  the  performance  of  this 
limiter  plus  APS  scheme  is  compared  with  the  previous  schemes  in  Section  3.6. 
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Chapter  4  describes  the  PAPR  suppression  schemes  on  synchronous  MC-CDMA. 
Section  4.1  introduces  MC-CDMA  and  its  PAPR  property.  There  are  synchronous  MC- 
CDMA  (down-Hnk)  and  asynchronous  MC-CDMA  (up-Hnk)  and  the  synchronous  MC- 
CDMA  signal  is  a  similar  PAPR  problem  to  an  OFDM  signal.  Section  4.2  investigates 
the  previous  schemes  to  reduce  the  PAPR  of  synchronous  MC-CDMA.  Among  them, 
MC-SS  looks  like  a  pretty  good  scheme  but  it  doubles  the  bandwidth.  The  characteristic 
of  Wash  Hadamard  (WH)  and  Golay  complementary  (GC)  spreading  sequences  for 
synchronous  MC-CDMA  is  studied  in  Section  4.3.  The  PAPR  property  of  a  synchronous 
MC-CDMA  signal  with  GC  sequences  is  similar  to  that  of  an  OFDM  signal.  Therefore, 
an  OFDM  PAPR  suppression  scheme,  limiter  plus  APS,  is  applied  to  the  synchronous 
MC-CDMA  system  with  GC  sequences  in  Section  4.4.  The  performance  of  the  scheme  is 
also  shown  in  that  section.  The  limiter  plus  APS  shows  effective  PAPR  suppression  and 
no  bandwidth  expansion  as  in  OFDM. 

Chapter  5  shows  the  PAPR  suppression  schemes  on  asynchronous  MC-CDMA. 
The  PAPR  property  of  asynchronous  MC-CDMA  is  described  in  Sections.  1.  The  PAPR 
depends  on  the  spreading  sequence.  So  the  previous  works  on  PAPR  suppression  try  to 
find  a  good  spreading  sequence.  The  Zadoff-Chu  sequence  was  the  best  one  for  low 
PAPR  and  good  BER.  Section  5.3  introduces  the  spreading  sequence  and  describes  its 
property.  Its  PAPR  is  OdB  in  terms  of  discrete  PAPR  but  6.4dB  in  terms  of  continuous 
PAPR.  So  the  APS  algorithm  is  applied  on  the  Zadoff-Chu  sequence  and  a  new  PAPR 
suppressed  Zadoff-Chu  (PSZ)  sequence  is  developed  in  Section  5.4.  The  PSZ  generates  a 
2-3dB  lower  PAPR  signal  and  has  good  BER  performance. 
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Chapter  6  contains  the  conclusion  of  this  dissertation  in  Section  6.1  and  the  future 
work  in  Section  6.2.  The  amplitude  limiter  and  the  APS  algorithm  can  be  applied  to  most 
of  the  signals  that  include  pulse  shaping  processes  in  digital  communication.  They  are 
effective  and  bandwidth  efficient  PAPR  suppression  schemes. 


CHAPTER  2 
PAPR  SUPPRESSION  IN  n/4  QPSK 

The  n/4  QPSK  modulation  is  a  popular  bandwidth  and  power  efficient  modulation 

format  for  portable  wireless  communications.  Furthermore,  the  power  efficiency  of  the 

71/4  QPSK  can  be  improved  by  PAPR  suppression.  For  bandwidth  efficient  PAPR 

suppression,  an  APS  algorithm  is  introduced.  The  major  problem  of  the  APS  algorithm  is 

high  computational  complexity.  The  problem  will  be  solved  by  an  iteration  reduction 

method  and  a  look-up  table  method.  The  complexity  of  the  look-up  table  will  be  also 

reduced  by  a  look-up  table  size  reduction  method.  For  maximal  PAPR  suppression,  the 

hybrid  scheme  that  is  the  combination  of  the  APS  algorithm  and  the  look-up  table  will  be 

proposed.  It  is  the  goal  of  this  chapter  to  suppress  the  PAPR  of  n/4  QPSK  for  power 

efficiency  improvement  without  any  bandwidth  efficiency  degradation. 

2.1  Introduction  to  n/4  QPSK  Modulation  Format 
The  n/4  QPSK  modulation  format  was  first  introduced  by  Barker  [15]  in  1962  as 
a  compromise  between  QPSK  and  OQPSK.  OQPSK  is  a  modified  QPSK  that  has  Ts/2 
time  offset  between  in-phase  channel  (I-channel)  symbols  and  quadrature-phase  channel 
(Q-channel)  symbols,  where  Ts  is  a  symbol  interval.  The  I-channel  is  the  real  part  of 
complex  envelope  representation  and  the  Q-channel  is  the  complex  part  of  it.  If  OQPSK 
is  used  with  a  half  sinusoid  pulse,  then  it  is  identical  to  MSK  that  is  a  constant  envelope 
modulation.  Let  a  data  symbol  be  ak+jbk.  where  k  is  a  symbol  number.  As  can  be  seen  in 
Figure  2.1,  the  data  symbol  {ak+jbk}  is  {±1,  ±j}  or  {e-^'"\  e-j^"'^}  for  QPSK.  The  signal 
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Figure  2.2:  Signal  Trajectory 
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Figure  2.3:  The  PAPR  of  7t/4,  offset  and  normal  QPSK  with  RRC 


constellation  of  OQPSK  looks  the  same  as  that  of  QPSK,  but  the  symbol  lasts  only  the 
half  of  a  symbol  interval  (Ts/2)  because  ak  changes  at  time  t=  0,  Ts,  2Ts,  ...  and  bk 
changes  at  t=  TJ2,  3TJ2,  5TJ2,  ....  For  n/4  QPSK,  the  data  symbol  {ak+jbk}  is  {±1,  ±j} 
for  odd  k  and  {e*^"''*,  e*^''"''*}  for  even  k.  So  the  maximum  phase  change  of  n/4  QPSK  is 
135°,  while  that  of  QPSK  is  180°.  The  maximum  phase  change  of  OQPSK  is  only  90°, 
but  it  changes  two  times  during  a  symbol  interval.  Therefore,  there  is  no  near  zero 
amplitude  signal  for  k/4  QPSK  and  OQPSK  as  seen  in  Figure  2.2.  The  dynamic  range  of 
OQPSK  is  the  smallest  and  that  of  n/4  QPSK  is  also  small.  However,  the  dynamic  range 
of  QPSK  is  infinite  because  there  are  zero  amplitude  signals. 
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According  to  the  paper  of  Ariyavisitakul  and  Liu  [16],  a  linear  amplifier,  such  as 
class  A  or  class  AB  or  class  B,  does  not  have  an  exact  linear  characteristic.  It  can  not 
amplify  a  high  peak  signal  or  a  low  amplitude  signal  without  any  distortion  because  of  its 
nonlinearity.  Therefore,  a  signal  with  small  dynamic  range,  which  means  low  PAPR  and 
large  minimum  signal  amplitude,  is  desirable  even  for  a  linear  amplifier.  As  can  be  found 
in  Figure  2.2  and  Figure  2.3,  OQPSK  and  n/4  QPSK  have  no  near  zero  amplitude  signals 
and  low  PAPR.  So  OQPSK  and  n/4  QPSK  are  better  suited  than  QPSK  for  a  linear 
amplifier.  The  power  spectral  density  (PSD)  and  the  BER  of  7t/4  QPSK  and  OQPSK  are 
the  same  with  those  of  QPSK  when  they  adopt  the  same  pulse  shape  in  a  linear  channel. 
Therefore,  they  have  all  the  same  bandwidth  and  average  power  efficiency.  The 
performance  of  n/4  QPSK  is  better  than  OQPSK  in  a  fast  fading  channel  [1 1]  [12]  [13], 
as  reviewed  in  Section  1.3. 

Therefore,  we  conclude  that  n/4  QPSK  is  the  best  compromise  of  power  and 
bandwidth  efficiency  in  a  general  portable  wireless  system.  So  the  APS  algorithm  will  be 
applied  to  n/4  QPSK  modulation  format  to  suppress  the  PAPR.  Motorola,  our  sponsor, 
wanted  to  apply  the  peak  suppression  scheme  to  a  n/4  QPSK  signal  with  RRC  (roll-off 
factor,  a=  0.35,  0.20)  pulse.  The  window  length  of  this  RRC  pulse  is  8Ts  for  a=  0.35  and 
12Ts  for  a=  0.20  to  satisfy  the  -50dB  adjacent  channel  interference.  As  can  be  found 
from  the  PSD  of  RRC  in  Figure  1 .4.  the  -50dB  interference  can  be  measured  at  the  center 
of  an  adjacent  channel,  where  a  channel  bandwidth  is  (l+a)/Ts.  If  a  bandwidth  efficient 
window  such  as  a  Blackman  window  is  applied  instead  of  the  rectangular  window,  the 
adjacent  channel  interference  can  go  down  to  below  -60dB. 
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2.2  Introduction  to  APS  Algorithm 

As  mentioned  in  Section  1 .2,  PSK  has  considerable  PAPR  when  it  is  used  with 
overlapped  RRC  pulses.  For  example,  as  can  be  seen  in  Figure  2.4,  the  peak  signal  is 
constituted  by  two  adjacent  overlapped  pulses  of  identical  symbols  in  binary  PSK 
(BPSK)  modulation.  A  signal  processing  such  as  a  Hamming  window  or  a  butterworth 
filter  can  be  applied  to  suppress  the  peak  signal.  Then,  power  efficiency  and  bandwidth 
efficiency  will  be  degraded  because  of  the  signal  distortion  and  sidelobe  regeneration.  To 
reduce  the  amount  of  these  efficiency  degradations,  the  signal  processing  can  be  applied 
where  there  is  a  peak  signal  and  it  leaves  the  rest  of  the  signal  alone.  However,  there  is 
still  considerable  sidelobe  regeneration,  which  degrades  bandwidth  efficiency.  To 
suppress  the  PAPR  of  the  signal  without  loosing  any  spectral  efficiency,  the  APS 
algorithm  was  introduced  by  Miller  and  O'Dea  [1]  [14]. 

The  bandwidth  of  the  signal  is  only  dependent  on  the  pulse  shape  in  PSK  as  long 
as  the  data  symbols  are  uncorrelated  and  symbol  interval  is  not  changed.  So,  if  the  pulse 
shape  is  not  changed,  the  signal  has  the  same  PSD  and  there  is  no  bandwidth  efficiency 
degradation.  The  APS  algorithm  changes  the  amplitude  of  each  symbol  that  causes  the 
peak  signal.  The  PAPR  of  the  signal  is  suppressed  by  multiplying  an  amplitude 
adjustment  factor  to  the  symbol.  This  scheme  does  not  change  any  pulse  shape  but  only 
changes  the  amplitude  of  pulses.  For  example,  as  can  be  seen  in  Figure  2.5,  the  symbols 
[+1,  +1]  are  multiplied  by  the  amplitude  adjustment  factors  [.8,  .8];  then  the  symbols 
become  [+.8,  +.8]  and  the  peak  signal  is  suppressed.  If  the  peak  is  higher,  smaller 
amplitude  adjustment  factors  will  be  multiplied  to  the  symbols  that  cause  the  higher  peak. 
On  the  other  hand,  to  the  symbol  that  does  not  cause  a  peak  signal,  an  amplitude 
adjustment  factor  1  will  be  multiplied,  so  it  does  not  change  the  amplitude  of  the  symbol. 
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Figure  2.4.  Peak  signal 


Figure  2.5.  Suppressed  peak  signal 


There  is  no  spectral  efficiency  loss  in  this  PAPR  suppression  scheme;  however,  there  is 
slight  loss  in  average  power  efficiency.  The  loss  is  about  0.1 4dB  for  n/4  QPSK  in  a 
Rayleigh  fading  channel  [1].  The  APS  algorithm  is  developed  to  find  the  amplitude 
adjustment  factors  to  suppress  the  PAPR  of  the  signal. 

A  summary  of  the  APS  algorithm  is  as  follows.  Generally,  a  PSK  signal  is 
represented  by  the  following  equation. 

S(t)=  X   dkp(t-kTs)  (2.1) 

k=-<» 

The  PAPR  suppressed  PSK  signal  by  the  APS  algorithm  can  be  expressed  as 

00 

S'(t)=  X  akdkp(t-kTs) 


(2.2) 


where  dk  is  a  data  symbol,  p(t)  is  a  Nyquist  RRC  pulse  shape  and  Ts  is  a  symbol  interval. 
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Figure  2.6:  Input  and  output  of  the  APS  algorithm 


The  input  of  the  APS  algorithm  is  a  data  symbol  dk  and  the  output  is  an  adjusted  data 
symbol  akdk,  as  seen  in  Figure  2.6. 

The  APS  algorithm  finds  the  amplitude  adjustment  factor  ak  for  each  data  symbol 
dk  through  following  steps. 
Step  1 :  Initially  set  ak  =1  for  all  k. 

Step  2:  Construct  S'(t)  from  the  data  dk  according  to  (2.2). 
Step  3:  Compute  the  peak  of  S'(t)  in  each  symbol  interval. 

Calculate  the  maximum  peak  value  (Sk)  and  its  location  (5k). 

Sk  =  max(  I  S'(t)  I  )  (k-l)Ts  <  t  <  kTs  (2.3) 

tk :  The  time  where  a  peak  is  located  in  (k-l)Ts  <  t  <  kTs 
Step  4:  Adjust  ak  according  to 

ak  ^  ak[  1 -pSkf(Sk)-P(l -6k.i)f(sk,i)]  (2.5) 
f(sk)  =  (Sk-Y)U(sk-y)  (2.6) 
where  U(.)  is  a  unit  step  function. 
P:  A  parameter  that  controls  the  suppression  rate.  It  determines  how  much  the 
algorithm  suppresses  a  peak  for  a  given  peak  value  and  its  location. 
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y:  A  parameter  that  determines  a  suppression  reference  level.  The  level 
determines  if  a  peak  is  high  enough  to  be  suppressed  or  not. 
Step  5:  Normalize  the  ak  to  make  E[ak  ]  =  1 . 

ak    <-  Sija  where  a  =  a^-  (2.7) 

Step  6:  Repeat  step  2  through  step  5  until  we  get  desirable  ak. 
Step  7:  Calculate  the  new  adjusted  data  akdk. 

In  step  3,  the  peak  amplitude  and  its  position  is  detected  in  each  symbol  interval. 
The  amplitude  adjustment  factor  ak  is  determined  by  these  two  values  Sk  and  5k.  In  step  4, 
the  factor  ak  is  smaller  as  the  peak  is  higher.  The  peak  is  usually  in  between  the  two 
symbols.  To  suppress  the  peak,  both  symbols  should  be  multiplied  by  the  proper  factors. 
To  suppress  the  peak  efficiently,  the  smaller  factor  should  be  multiplied  to  the  symbol 
that  is  closer  to  the  peak.  It  means  the  constant  ak  becomes  smaller  as  the  peak  is  close  to 
the  symbol  in  step  4.  Then,  the  total  power  will  be  decreased  through  this  PAPR 
suppression  because  every  value  of  ak  is  less  than  equal  to  one.  So  every  factor  ak  needs 
to  be  divided  by  the  root  mean-squared  (RMS)  value  of  the  ak  to  normalize  the  total 
signal  power. 

Through  this  APS  algorithm,  the  PAPR  of  the  signal  can  be  reduced  by  up  to  3dB. 
For  example,  the  PAPR  of  a  BPSK  signal  is  reduced  from  5.6dB  to  2.5dB  and  that  of  a 
7r/4  QPSK  signal  is  reduced  from  4.8dB  to  1 .9dB  for  RRC  (a=0.20)  pulse.  For  RRC 
(a=0.35)  pulse,  the  PAPR  of  a  BPSK  signal  is  reduced  from  3.9dB  to  2.2dB  and  that  of 
7t/4  QPSK  signal  is  reduced  from  3.2dB  to  1.7dB.  The  average  power  efficiency 
degradation  at  a  receiver  is  less  than  0.2dB  in  a  Rayleigh  fading  channel.  That  is  much 
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less  than  the  power  gain  at  a  transmitter.  The  3dB  gain  in  PAPR  corresponds  to  around  1- 
2dB  gain  in  a  linear  power  amplifier. 

However,  there  are  two  main  practical  problems  to  implement  the  APS  algorithm. 
Generally,  it  takes  15-30  iterations  to  get  an  appropriate  amplitude  adjustment  factor  ak. 
This  huge  computational  complexity  is  the  first  major  problem  of  the  algorithm.  A 
portable  communication  system  usually  has  a  digital  signal  processor  (DSP)  chip  that 
may  implement  the  algorithm  and  other  signal  processes.  However,  the  DSP  chip 
probably  can  not  afford  this  huge  computational  complexity  because  of  the  capacity  limit 
of  the  DSP  chip.  Another  problem  is  the  block  processing  in  step  5.  If  the  RMS  value  of 
ak  is  calculated  from  500  symbols,  then  the  time  delay  of  500Ts  will  be  caused.  This 
delay  may  make  the  real  time  processing  difficult. 

2.3  Previous  Work  on  Practical  Problems 
Fedon-Rovira  [17]  and  Miller  suggested  a  look-up  table  method  and  an  artificial 
neural  networks  (ANN)  method  to  avoid  the  above  two  problems  of  the  APS  algorithm. 
A  PSK  signal  is  constituted  by  the  sum  of  overlapped  RCC  pulses,  so  a  peak  signal  is 
constructed  by  a  certain  combination  of  adjacent  symbols.  If  a  truncated  RRC  pulse 
(a=0.35)  of  length  8Ts  is  used,  then  9  adjacent  symbols  (N=9)  need  to  be  considered  to 
find  the  peak  value  of  the  signal.  It  means  the  9  adjacent  symbols  give  the  major 
contribution  to  calculate  the  amplitude  adjustment  factor  ak,  which  will  be  multiplied  to 
the  center  symbol  (dk)  of  the  9  adjacent  symbols.  The  symbols  outside  of  the  9  symbols 
give  only  negligible  contribution  to  the  peak  value  and  the  ak.  For  a  RRC  pulse  (a=  0.20) 
of  length  12Ts,  13  adjacent  symbols  need  to  be  considered.  The  look-up  table  method 
uses  a  pre-calculated  value  of  ak  for  each  given  adjacent  symbol  stream.  The  peak  value 
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Figure  2.7:  Peak  suppression  by  using  look-up  table  or  ANN  (N=5) 


and  the  location  can  be  calculated  in  the  APS  algorithm  when  the  adjacent  symbols  are 
given.  Then,  an  appropriate  ak  can  be  pre-calculated  for  every  possible  adjacent  symbol 
stream  through  the  algorithm  with  sufficiently  large  number  of  iterations.  The 
computational  complexity  or  the  long  time  delay  is  not  a  problem  because  the  calculation 
is  not  performed  in  a  real  time  in  a  portable  communication  system.  Every  possible  ak  is 
ready  in  the  look-up  table  of  a  system  and  can  be  picked  up  when  a  corresponding 
symbol  stream  is  detected  as  in  Figure  2.7.  Then,  the  selected  ak  is  multiplied  to  the 
middle  symbol  of  the  symbol  stream. 

However,  this  look-up  table  method  also  has  a  practical  problem.  The  N=9  look- 
up table  means  that  9  symbols  are  considered  to  calculate  the  ak  of  the  look-up  table. 
Then,  there  are  2^  possible  streams  of  9  symbols  for  BPSK  and  2*4^  possible  streams  of  9 
symbols  for  n/4  QPSK.  The  look-up  table  size  is  reduced  to  1/2  for  BPSK  and  1/8  for  n/4 
QPSK  because  the  rotation  of  the  symbol  phase  does  not  change  the  peak  value.  So  2 
symbol  streams  of  BPSK  and  8  symbol  streams  of  n/4  QPSK  can  share  the  same  ak  value. 
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Figure  2.8:  The  structure  of  an  artificial  neural  networks  (ANN) 


Then,  the  look-up  table  size  is  2^"'  for  BPSK  and  4^"'  for  nlA  QPSK.  The  look-up  table 
size  of  2*  that  requires  256  bytes  memory  is  acceptable,  but  the  look-up  table  size  of  4* 
that  requires  64,000  bytes  memory  is  probably  too  big  for  a  portable  communication 
system. 

The  ANN  method  is  similar  to  the  look-up  table  method.  In  this  case,  a  nonlinear 
filter  that  is  constructed  by  a  feed-forward  ANN  is  used  instead  of  a  look-up  table  in 
Figure  2.7.  The  ANN  nonlinear  filter  is  pre-trained  by  using  the  best  look-up  table.  The 
filter  needs  only  a  small  memory  to  store  the  values  of  weights  (Wj/''*)  and  biases  (by)  in 
the  ANN,  as  can  be  seen  in  Figure  2.8.  However,  it  needs  a  considerable  computational 
complexity  for  a  large  size  ANN  filter.  Thus,  the  ANN  method  has  a  similar  problem  to 
the  look-up  table.  For  BPSK,  the  ANN  requires  an  acceptable  computational  complexity. 
However,  for  7i/4  QPSK,  the  ANN  requires  very  large  computation  because  such  a  big 
ANN  is  necessary  to  suppress  the  peak  sufficiently  in  Tr/4  QPSK. 
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Table  2.1 :  Suppressed  PAPR  by  hybrid  methods 


XT           1_  V 

Number  oi 

1  ^  C   1    A  XTXT 

IN— 4  1  auie 

IN — U  1  aUlC 

Iterations 

(71  bytes) 

(64  bytes) 

(1024  bytes) 

0 

2.40  dB 

2.46  dB 

2.11  dB 

1 

2.19  dB 

2.19  dB 

1.97  dB 

2 

2.06  dB 

2.01  dB 

1.89  dB 

3 

1.97  dB 

1.91  dB 

1.83  dB 

Therefore,  both  a  look-up  table  method  and  an  ANN  method  work  pretty  well  for  BPSK, 
but  they  require  large  memory  or  high  computational  complexity  for  n/4  QPSK  that  is  of 
greatest  interest.  Thus,  they  suggested  a  hybrid  method  that  is  the  combination  of  the 
look-up  table  and  the  APS  algorithm  or  the  combination  of  the  ANN  and  the  APS 
algorithm.  In  this  method  they  used  a  small  size  look-up  table,  whose  window  length  N  = 
4  or  6  symbols,  or  a  small  size  ANN  (12-5-1)  plus  the  algorithm  with  1-3  iterations.  The 
results  of  peak  suppression  performance  are  shown  in  Table  2.1  for  7i/4  QPSK  with  RRC 
(roll-off  =  0.35).  The  original  PAPR  of  the  n/4  QPSK  signal  is  3.18dB  and  the  suppressed 
PAPR  by  the  algorithm  with  30  iterations  is  1 .71  dB. 

It  was  a  pretty  good  result;  however,  it  was  difficult  to  say  which  one  is  the  best 
choice  because  we  didn't  know  how  much  calculation  was  needed  for  each  iteration  of 
the  algorithm.  Therefore,  it  is  necessary  to  calculate  how  much  computational  complexity 
is  required  for  each  algorithm  iteration. 

2.4  Research  on  APS  Algorithm 

2.4.1  Removing  Block  Processing 

As  mentioned  in  Section  2.2,  it  is  desirable  to  remove  the  block  processing  in  the 
APS  algorithm  because  it  causes  a  time  delay  which  makes  real  time  processing  difficult. 
In  step  5  of  the  APS  algorithm,  the  RMS  value  of  the  amplitude  adjustment  factor  is 
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required  to  be  calculated  to  normalize  the  total  signal  power.  Thus,  the  past  and  future 
values  of  ak  are  needed  to  calculate  each  current  value  of  ak.  This  problem  is  addressed  by 
using  block  processing  whereby  the  amplitude  adjustment  factors  for  all  symbols  in  a 
block  of  length  N  are  calculated  simultaneously.  The  block  processing  causes  a  time 
delay  of  N  (100-10,000)  symbol  intervals  (2.7),  which  may  be  a  big  problem  in  a  real 
time  communication. 

This  problem  can  be  solved  by  using  a  fixed  a,  that  is  a  very  long-term  RMS 

value  of  ak.  This  a,  is  pre-calculated  by  using  a  long  random  training  sequence.  We  can 

find  out  the  fixed  a,  by  averaging  a  large  number  of  a  which  is  calculated  by  large 
random  blocks  of  data.  For  example, 

—       -       1  - 

a,  =E[a]=  —  ^  [-9484  .9628  .9742  .9821  ],  (2.8) 

20  „=i 

where   a=    —V  a.^  .  Four  iterations  of  the  algorithm,  6=1.8,  7=1.22,  N=5,000  and 

RRC  (a=.35,  9  samples/symbol)  are  used  for  this  calculation.  The  long-term  RMS  value 

a,  has  a  standard  deviation  of  only  about  .002,  which  means  that  a,  is  a  stable  value. 
The  power  difference  of  Pj,  which  the  total  signal  power  after  the  i"^  algorithm  iterafions 
by  using  a, ,  is  negligible.  That  means  that  a,  provides  good  power  normalization.  For 
example, 

I  10  log  (Po  /  P4  )  I  <  .001  dB,        where  Pi  =  —  f  [S,{t)fdt .  (2.9) 

AT  * 

Therefore,  the  total  power  can  be  normalized  by  using  this  fixed  a  without  block 
processing.  The  average  power  loss  caused  by  this  method  is  about  0.03dB  with  0.02dB 
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Standard  deviation.  Now,  the  delay  of  this  algorithm  is  only  half  of  the  length  of  RRC 
pulse  (4Ts  or  6Ts)  for  each  iteration.  So  the  delay  is  only  16Ts  or  24Ts  for  the  algorithm 
with  4  iterations  instead  of  the  long  block  length  (100-10,OOOTs).  Then,  there  should  be 
no  problem  for  a  real  time  communication. 

2.4.2  Choosing  Parameters  P  and  y 

Through  the  research,  choosing  appropriate  values  of  the  parameters  y  and  P  turns 
out  to  be  the  best  way  to  reduce  the  high  computational  complexity  of  the  APS  algorithm. 
The  APS  algorithm  suppresses  the  peak  signal  when  the  normalized  value  of  the  peak  is 
larger  than  the  parameter  y  as  in  step  4  of  the  algorithm.  The  parameter  P  controls  the 
suppression  rate,  that  is,  how  much  the  algorithm  suppresses  the  peak  for  the  same  peak 
value  and  position  given.  Therefore,  choosing  the  proper  P  and  y  is  the  major  key  for 
efficient  PAPR  reduction.  By  efficient  PAPR  suppression,  the  number  of  iterations  can  be 
reduced  and  hence  the  computational  complexity  is  reduced.  The  process  to  find  proper 
values  of  P  and  y  is  as  follows. 

Step  1 :  Generate  the  maximum  peaks.  When  the  n/4  QPSK  symbol  is 


I  =  [1  0  0  1],  Q  =  [1  0  1  0]  generates  a  maximum  peak  of  2.75dB  PAR,  and  1  = 
[1  1  1  0  0  1  1  1],  Q  =  [0  0  1  0  1  0  1  1]  generates  a  maximum  peak  of  3.15dB  PAR. 
Step  2:  While  changing  the  p  and  y  values,  see  how  much  the  above  maximum  peaks  can 
be  suppressed.  It  is  important  to  vary  the  P  and  y  values  up  to  enough  range.  Then, 
find  some  candidates  of  the  P  and  y  values  which  suppress  both  peaks  very  well. 


(2.10) 
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Step  3:  Test  the  candidates  of  the  p  and  y  values  for  the  peak  signal  that  is  constituted  by 
1,000-10,000  random  data  and  choose  the  best  one. 

In  the  original  paper  (1994)  of  the  APS  algorithm  that  is  a  draft  of  [1],  the 
parameter  P  was  restricted  to  be  in  the  range  [.5,  1]  for  stable  and  adaptive  peak 
suppression  and  typically,  p  =  .6  and  y  ^  1.2  were  used  for  RRC  (a=.35)  without 
considering  the  number  of  iterations  of  the  algorithm.  However,  a  larger  value  of  P  than 
the  restriction  is  required  to  suppress  the  peak  effectively  for  the  algorithm  with  a  small 
number  of  iterations.  For  example,  a  good  choice  of  P  and  y  is  1.8  and  1 .22  respectively 
for  the  algorithm  with  4  iterations,  while  a  good  choice  of  p  and  y  is  0.9  and  1.19 
respectively  for  the  algorithm  with  1 5  iterations.  However,  the  parameter  P  should  not  be 
too  large  because  a  large  value  of  p  causes  BER  degradation,  which  means  the  system 
loses  power  efficiency.  The  PAPR  can  be  reduced  from  3.12  to  1.95  dB  by  the  algorithm 
with  only  4  iterations  when  p  and  y  is  1.8  and  1.22  respectively.  However,  the  PAPR  can 
only  be  reduced  from  3.12  to  2.32  dB  by  the  same  algorithm  with  the  same  number  of 
iterations  when  P  and  y  is  0.9  and  1.19  respectively.  It  also  takes  more  than  10  iterations 
to  reduce  the  PAPR  from  3.12  to  1.95dB  when  P  and  y  is  0.9  and  1.19  respectively. 
Therefore,  by  using  a  proper  selection  of  the  parameters,  the  PAPR  can  be  suppressed 
efficiently,  so  the  number  of  iterations  and  the  computational  complexity  can  be  reduced. 

Another  important  factor  to  be  considered  to  find  proper  values  of  p  and  y  is  the 
power  efficiency  loss  r|.  The  power  efficiency  loss  is  measured  by  the  difference  between 
the  value  of  Eb/No  of  the  suppressed  signal,  ((Eb/No)sp),  and  the  value  of  Eb/No  of  the 
original  signal,  (Eb/No),  to  achieve  the  same  BER.  This  is  calculated  easily  in  a  Rayleigh 
fading  channel  model  because  the  required  Eb/No  is  approximately  a  linear  function  of 
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BER,  Pe.  Then,  the  power  loss  r]  is  the  same  at  all  BER's  in  a  Rayleigh  fading  channel 
[14]  and  can  be  written  by 

JI^^(X^P^^^_l^  (2.11) 

If  the  PAPR  of  a  signal  is  reduced,  then  the  power  efficiency  gain  will  be  achieved  from 
Figure  1.2.  However,  the  noise  sensitivity  power  loss  r\  degrades  the  power  gain  that  is 
achieved  by  PAPR  suppression.  Therefore,  parameters  P  and  y  that  make  the  suppressed 
PAPR  (SP)  small  do  not  need  to  be  found,  but  the  parameters  that  make  the  suppressed 
PAPR  plus  the  power  loss  r\  times  a  constant  c  (SP+r|c)  small  need  to  be  found.  The 
constant  c  is  decided  by  the  efficiency  characteristic  of  the  power  amplifier  in  Figure  1 .2. 
Let  c  equal  to  1  for  simplicity.  For  example,  for  the  RRC  (a=.20)  and  an  algorithm  with  4 
iterations,  (P,  y)  =  (1.9,  1.22)  gives  SP  =  2.222dB,  r\  =  .114dB,  then  SP+r]  =  2.336dB, 
while  (p,  y)  =  (1.7,  1.18)  gives  SP  =  2.255dB,  ri  =  .046dB,  then  SP+ti  =  2.301dB.  The  SP 
of  the  first  choice  of  parameters  is  smaller  than  that  of  the  second  one,  but  the  second  one 
should  be  selected  because  the  SP+r|  of  the  second  one  is  smaller  than  that  of  the  first 
one.  Good  choices  of  P  and  y  are  in  Table  2.2.  In  Figure  2.9,  the  new  suppression 
performance  by  these  good  choices  of  parameters  is  compared  with  the  previous 
suppression  performance  by  the  fixed  parameters  p  and  y  which  are  0.6  and  1.2 
respectively.  The  solid  line  is  for  chosen  parameters  and  the  dotted  line  is  for  fixed 
parameters.  The  notation  'x'  is  for  the  RRC  (a  =  .20)  pulse  and  'o'  is  for  the  RRC  (a  = 
.35).  It  shows  that  the  performance  of  the  4  iterations  with  good  choice  of  parameters  is 
better  than  that  of  10  iterations  with  fixed  parameters. 
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Table  2.2:  Good  choices  of  P  and  y  for  various  iteration  numbers 


1  1 

2 

3 

4 

1 5  iterations 

a  =  .35 

P 

1.7 

1.8 

2.0 

1.8 

0.9 

Y 
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1.23 

1.23 

1.22 

1.19 

a  =  .20 

P 
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2.1 

2.2 

1.9 

1.2 

Y 

1.18 

1.25 

1.26 

1.235 

1.225 

5 


1.5'  '  '  ' 

0  5  10  15 

Number  of  Iterations 


Figure  2.9:  The  PAPR  suppression  performance  of  the  APS  algorithm 


2.4.3  Choosing  Number  of  Samples 

Another  way  to  reduce  the  computational  complexity  is  by  using  a  small  number 
of  samples  per  symbol  since  a  large  number  of  samples  causes  high  computational 
complexity  in  the  APS  algorithm.  In  the  original  paper  of  the  APS  algorithm,  they  used 
25  samples  per  symbol  interval  when  the  signal  S'(t)  of  (2.2)  is  constructed  in  step  2  of 
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the  algorithm.  Then  the  discrete  signal  has  an  accurate  representation  of  almost  every 
peak  value  of  the  analog  signal.  However,  using  that  many  samples  causes  a  lot  of 
computational  complexity.  If  5  samples  are  used  instead,  then  the  number  of 
computations  can  be  substantially  reduced,  but  some  peak  values  may  be  missed  when 
the  peak  values  falls  between  sampling  points.  This  may  cause  improper  reduction  of  the 
peak  signal  and  the  reduction  capability  of  the  algorithm  will  be  degraded.  Therefore, 
there  is  a  trade  off  between  computational  complexity  and  PAPR  suppression 
performance  in  choosing  a  proper  sampling  rate. 

There  is  another  factor  involved  in  choosing  a  proper  sampling  rate.  The 
maximum  peak  often  falls  at  the  middle  of  the  symbol  interval  Ts.  So  there  should  be  a 
sampling  point  at  the  middle  of  symbol  interval  to  sample  the  maximum  peak  value 
accurately.  If  a  sampling  point  is  located  at  the  middle  with  an  even  number  of  samples 
per  symbol,  then  the  sampling  is  not  symmetric.  This  may  be  a  problem  when  a  small 
number  of  samples  is  used  because  the  peak  reduction  rate  of  adjacent  left  symbol  will  be 
different  from  the  peak  reduction  rate  of  right  symbol.  This  also  may  degrade  the 
performance  of  the  algorithm.  Therefore,  it  is  better  to  use  an  odd  number  of  samples  to 
reduce  the  peak  effectively. 

According  to  some  performance  tests,  7  or  9  samples  per  symbol  interval  are  the 
best  trade  off  between  computational  complexity  and  PAPR  suppression  performance,  as 
can  be  seen  in  Figure  2.10.  In  this  simulation,  the  parameter  SP+r]  is  used  instead  of  SP 
for  the  same  reason  as  in  Section  2.4.2.  The  algorithm  with  15  iterations  and  parameters 
(P=.6,  y=1.2)  is  used  for  this  simulation. 
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Figure  2.10:  The  effect  of  sampling  number  on  peak  suppression  performance 


It  should  be  noted  that  a  large  enough  number  of  samples  should  be  used  to 
construct  the  peak  suppressed  S'(t)  after  the  step  7  to  measure  the  SP  accurately  in  a 
computer  program  such  as  Matlab.  Otherwise,  a  smaller  SP  value  than  a  real  SP  value 
will  be  detected.  However,  the  smaller  SP  is  a  fake  since  the  real  peak  value  may  be 
missed  because  of  the  small  sampling  rate. 

2.4.4  Computational  Complexity  of  the  Algorithm 

Our  sponsor,  Motorola,  wanted  to  apply  the  APS  algorithm  to  a  portable 
communication  system.  However,  the  system  could  only  afford  about  2  or  3  million 
instructions  per  second  (Mips).  So  we  needed  to  know  how  much  complexity  is  required 
for  the  algorithm  with  N  iterations  and  P  samples/symbol.  To  measure  the  complexity  for 
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Table  2.3:  Computational  complexity  in  terms  of  Flops  in  Matlab 


Step 

Flops/bit  (for  1000  random  data) 

A 

365  (22.3%) 

B 

72  (4.4%) 

C 

28(1.7%) 

D 

2(0.12%) 

E 

1167  (71.2%) 

F 

4  (0.24%) 

Total 

1638  (100%) 

a  given  value  of  N  and  P,  how  many  computations  are  required  for  each  step  of  the 
algorithm  was  investigated.  The  steps  are  as  follows. 

Step  A:  Generating  a  discrete  signal  of  P  samples/symbol  from  the  data  (dk).  [=  Step  2] 
Step  B:  Finding  maximum  values  and  their  locations  in  a  symbol  interval  (Ts).  [=  Step  3] 
Step  C:  Calculating  a^.  [=  Step  4] 
Step  D:  Normalizing  a^.  [=  Step  5] 

Step  E:  Generating  a  suppressed  signal  of  P  samples/symbol. 

[=2nd  or  later  iteration  of  Step  2] 
Step  F:  Multiplying  a^  and  dk.  [=  Step  7] 

First,  the  complexity  was  calculated  in  floating-point  operations  (Flops)  in  the 
Matlab  program.  However,  the  counting  of  Flops  is  a  little  bit  different  from  the  counting 
of  Mips  in  a  DSP  chip.  The  results  for  the  4-iteration  algorithm  with  9  samples/symbol 
are  shown  in  Table  2.3.  It  is  noted  that  most  of  the  complexity  comes  from  steps  A  and  E, 
which  involve  generating  a  discrete  signal  by  pulse  shaping.  Since  there  are  redundant 
calculations  in  steps  A  and  E  in  the  Matlab  program,  we  tried  to  count  the  necessary 
number  of  calculations  theoretically  as  follows. 
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The  truncated  RRC  (a=.35)  pulse  of  length  8Ts  was  adopted  when  a  signal  is 
generated  by  pulse  shaping  for  peak  suppression.  Therefore,  the  8  individual  time  shifted 
pulses,  p(t-kTs)  for  k  =  -3,-2,  ...  ,3,4,  are  added  during  a  symbol  interval,  Tj.  Then,  the 
equation  (2.2)  will  become  the  following  equation. 

4 

S(t)=  ^  akdkp(t-kTs),  where  0<t<Ts.  (2.12) 

If  the  signal  is  sampled  when  the  sampling  rate  is  P  samples  per  symbol  interval,  then  the 
discrete  signal  equation  will  become  the  following  equation. 

S(mTs/P)=  X  akdkp(mTs/P-kTs),         m=l,2,...,P.  (2.13) 

S(mTs/P)  can  be  calculated  by  7  additions  and  16  multiplications  in  (2.13).  However,  the 
1 6  multiplications  are  actually  24  calculations  because  2  multiplications  (akxdkxp(t))  of  1 
complex  number  (dk)  with  2  real  numbers  (ak,  p(t))  become  3  calculations.  Then,  7+24  = 
31  calculations  are  needed  for  each  sample.  It  means  3 IP  calculations  are  required  per 
every  symbol  interval  Ts  in  step  E.  In  step  A,  only  (7+16)P  =  23P  calculations  are 
required  because  every  ak  is  set  to  one.  In  step  B,  we  have  to  calculate  the  absolute  values 
of  the  discrete  signal  to  find  a  maximum  value  and  its  location  in  every  symbol  interval. 
It  requires  4  calculations  for  every  sample  because  each  signal  component  is  complex. 
Then,  4P  calculations  are  needed  per  every  symbol  interval  in  step  B.  It  should  be  noted 
that  there  must  be  some  additional  complexity  that  is  used  to  find  the  biggest  peak  value 
and  its  location.  But  the  additional  complexity  is  excluded  here.  In  step  C,  there  are  2 
decisions,  3  addifions  and  7  multiplications  per  symbol,  where  the  I  decision  is  2 
calculations.  So  a  total  of  14  calculations  are  needed  per  every  symbol  interval  in  step  C. 
There  are  1  mukiplication  that  is  1  calculation  in  step  D  and  1  complex  multiplication 
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Table  2.4:  Minimum  number  of  calculations  for  each  step 


Step 

A 

B 

C 

D 

E 

F 

Total 

Calculations/  bit 

11. 5P 

2P 

7 

0.5 

15. 5P 

1 

29P+8.5 

that  requires  2  calculations  in  step  F.  Hence,  the  minimum  necessary  calculations  in  each 
step  are  as  in  Table  2.4.  It  should  be  noted  that  1  symbol  is  2  bits  in  the  n/4  QPSK 
modulation  format. 

For  example,  the  APS  algorithm  with  9  samples  per  symbol  and  4  iterations  has  1 
step  A  and  F,  3  steps  E  and  4  steps  B,  C  and  D.  So  the  algorithm  requires  (1 1.5*9+1)  + 
(15.5*9)*3  +  (2*9+7+0.5)*4  =  625  calculations  per  bit  to  implement.  If  the  transmission 
speed  is  9,600  bit/sec,  it  requires  at  least  6.000  Mips  to  implement  the  algorithm  in  a  DSP 
chip.  If  the  algorithm  is  used  with  7  samples  per  symbol  and  1  iteration,  then  (1 1.5*7)  + 
(1)  +(2*7+7+0.5)  =  103  calculations  are  required.  Thus,  the  algorithm  requires  0.989 
Mips  to  be  implemented  in  a  DSP  chip.  Practically,  the  algorithm  must  require  more  than 
this  number  of  calculations  when  it  is  implemented  in  a  real  DSP  chip. 

2.4.5  The  Performance  of  the  APS  algorithm 

Table  2.5  shows  simulation  results  for  random  data  of  5000  bits.  The  algorithm 

uses  9  samples  per  symbol  and  proper  parameters  p  and  y.  These  parameters  are  decided 
experimentally  for  each  number  of  iterations  respectively.  So  the  parameters  are  different 
for  the  APS  algorithm  with  different  number  of  iterations.  Therefore,  it  should  be  noted 
that  these  results  are  not  the  results  of  the  APS  algorithm  with  the  same  parameters  while 
the  number  of  iterations  changes.  The  original  peak  is  3.09  dB  for  RRC  (a  =  .35)  and 
4.75dB  for  RRC  (a  =  .20)  and  the  suppressed  peak  (SP)  is  as  shown  in  Table  2.5. 
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Table  2.5:  The  example  of  suppression  results  by  using  proper  (3,  y 


1 
1 

L 

J 

A 
H- 

1  S  it*^t*!3tir\Tic 
IJ  llCldLlUllo 

a  =  .20 

Loss  r| 

u.zv 

U.Z7 

KJ.JZ. 

U..J  / 

SP 

2.47 

2.31 

2.23 

2.16 

2.03  dB 

SP  +  Tl 

2.76 

2.60 

2.55 

2.54 

2.47  dB 

a  =  .35 

Loss  r| 

0.03 

0.07 

0.10 

0.11 

0.15  dB 

SP 

2.38 

2.12 

2.03 

1.95 

1.79  dB 

SP  +  Tl 

2.41 

2.19 

2.13 

2.06 

1.94  dB 

Complexity  (a  =  .35) 

1.25 

2.83 

4.41 

6.00 

23.4  Mips 

4.5 


^  A 

w 
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Figure  2.1 1 :  The  Maximum,  average  and  minimum  suppressed  PAPR 


It  should  be  noted  that  SP+r|  is  a  more  important  parameter  than  SP  to  measure 
the  performance  of  the  algorithm.  This  result  may  vary  slightly  according  to  the  random 
data  since  the  PAPR  suppression  performance  also  depends  on  the  random  data. 
However,  it  is  evident  that  the  algorithm  with  4  iterations  gives  pretty  competitive 
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reduction  performance  with  much  smaller  complexity  than  that  of  the  algorithm  with  15 
iterations.  The  algorithm  with  2  or  3  iterations  also  gives  reasonable  reduction 
performance  with  smaller  complexity.  For  a  RRC  of  roll-off  factor  a  =  .20,  the  reduction 
performance  of  1  iteration  is  about  85%  of  that  of  15  iterations,  while  the  calculation 
amount  of  1  iteration  needs  only  5.3%  of  that  of  15  iterations.  Furthermore,  the 
performance  of  the  algorithm  with  a  small  number  of  iteration  is  better  than  it  looks 
because  power  loss  r\  is  decreasing  while  reduction  performance  is  decreasing.  Figure 
2.1 1  shows  the  suppressed  PAPR  of  ti/4  QPSK  from  50  realizations  of  the  APS  algorithm 
for  a  random  data  block  of  5000  bits.  The  solid  line  is  the  average  of  the  peaks,  the  dash- 
dot  line  is  the  maximum  of  the  peaks,  and  the  dotted  line  is  the  minimum  of  the  peaks. 

2.4.6.  Other  schemes  for  Computational  Complexity  Reduction 

As  seen  in  Table  2.4,  step  A  and  step  E  have  over  80%  of  the  total  algorithm 
computation.  So  the  steps  A  and  E  are  the  places  where  the  computation  reduction  should 
be  conducted  to  reduce  the  total  algorithm  computation  effectively.  The  computations  of 
steps  A  and  E  are  multiplying  ak  to  RRC  pulses  and  adding  the  overlapped  pulses.  The 
results  of  left  half  side  multiplications  and  those  of  right  half  side  multiplications  are  the 
same  because  the  truncated  RRC  pulse  p(mTs/P)  in  (2.13)  is  symmetric.  Therefore,  the 
results  of  the  multiplications  for  the  left  side  pulse  can  be  used  again  for  the  right  side 
pulse.  Then,  50%  of  the  multiplication  computation  and  37%  of  the  total  algorithm 
computation  can  be  reduced  for  the  APS  algorithm  with  9  samples  and  4  iterations.  This 
scheme  does  not  cause  any  performance  degradation. 

Another  possible  scheme  to  reduce  the  algorithm  computation  is  a  varying 
window  method.  The  original  window  length  for  truncated  RRC  is  determined  to  provide 
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-50dB  adjacent  channel  interference.  So  the  window  lengths  of  12Ts  and  8Ts  are  required 
for  RRC  pulses  (a=.20  and  a=.35)  respectively.  However,  these  big  window  lengths  are 
not  necessary  during  the  initial  algorithm  stages  because  the  -50dB  adjacent  channel 
interference  is  determined  at  the  stage  of  converting  digital  signal  to  analog  signal  after 
the  algorithm.  Using  a  small  window  reduces  the  algorithm  computation,  but  it  generates 
less  accurate  ak.  However,  very  accurate  ak  are  not  necessary  in  the  first  few  iterations  of 
the  algorithm,  while  the  last  few  iterations  need  very  accurate  ak  for  best  peak 
suppression.  For  example,  the  algorithm  with  3  iterations  and  RRC  (a=.20),  whose 
window  lengths  are  8Ts,  lOTs  and  12Ts  instead  of  12Ts,  12Ts  and  12Ts,  has  a  37%  of 
calculation  reduction  with  only  O.OSdB  less  peak  suppression.  The  algorithm  with  4 
iterations  and  RRC  (a=.35),  whose  window  lengths  are  4Ts,  4Ts,  6Ts  and  8Ts  instead  of 
8Ts,  8Ts,  8Ts  and  8Ts,  has  a  31%  calculation  reduction  with  only  O.OSdB  less  peak 
suppression. 

The  combination  of  both  calculation  reduction  schemes  makes  the  total 
calculation  about  half.  For  example,  the  minimum  calculations  for  the  APS  algorithm 
with  3  iterations,  9  samples  per  symbol  and  RRC  (a=.20)  is  reduced  from  8.47  Mips  to 
3.03  Mips  with  O.OSdB  degradation  in  peak  suppression.  The  minimum  calculations  for 
the  APS  algorithm  with  4  iterations,  9  samples  and  RRC  pulse  (a=.3S)  is  reduced  from 
6.00  Mips  to  3.11  Mips  with  O.OSdB  degradation  in  peak  suppression.  Therefore,  the 
computational  complexity  is  reduced  from  23  Mips  to  3  Mips  for  RRC  (a=.3S)  at  a  cost 
of  less  than  0.2dB  in  peak  suppression  loss  by  using  proper  parameters,  the  symmetric 
property  and  the  varying  window  method. 
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2.5  Research  on  Look-up  Table  Method 
2.5.1  Look-up  Table  Size  Reduction 

A  look-up  table  can  be  a  substitute  for  the  APS  algorithm  in  a  portable  wireless 
communication  system.  There  is  relatively  very  small  computational  complexity  in  the 
look-up  table  method,  but  a  lot  of  memory  is  required  to  suppress  the  peak  signal 
maximally.  As  mentioned  in  Section  2.3,  an  N=9  look-up  table  is  necessary  to  suppresses 
the  peak  signal  maximally  for  the  truncated  RRC  pulse  of  length  8Ts.  The  N=9  look-up 
table  has  2*4^  entries  for  n/4  QPSK.  The  table  size  can  be  reduced  to  4*  (=2*4^/8) 
because  the  streams  of  9  symbols  whose  phases  are  rotated  by  a  multiple  of  the  7i/4  have 
the  same  peak  value.  Let's  call  this  the  first  reason  of  table  size  reduction.  For  a  RRC 
pulse  (a=.20)  with  length  12Ts,  the  N=13  look-up  table  has  2*4*^  table  entries  for 
maximum  peak  suppression;  however,  the  table  size  can  be  reduced  to  4'^  (=2*4'V8) 
because  of  the  first  reason.  Due  to  the  memory  limitation  in  a  DSP  chip,  the  N=4-7  look- 
up table  is  mainly  considered  even  though  these  can  not  suppress  the  signal  maximally. 
Let's  explain  1/8  size  reducfion  through  an  example.  Let's  denote  a  n/4  QPSK  symbol 
exp(j7rw/4)  as  w  for  simplicity  where  we 0,1, 2,  ...  ,7  as  can  be  seen  in  Figure  2.14.  Then, 
the  symbol  stream  w  =  [0121]  indicates  exp(j7:w/4)  =  [1  exp(j7r/4)  exp(j7i/2)  exp(j7r/4)]  in 
the  N=4  look-up  table.  The  8  symbol  streams  rotated  by  multiples  n/4  are  [0121]  [1232] 
[2303]  ...  [7010].  These  symbol  streams  generate  the  same  peak  value,  so  they  need  only 
one  ak  value  in  the  look-up  table.  Then,  the  look-up  table  size  can  be  reduced  by  factor  of 
1/8  because  of  the  first  reason  of  the  table  size  reduction. 
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Figure  2.12.  The  signal  constellation  of  a  n/4  QPSK  signal 
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Figure  2.13:  The  symbol  streams  that 
generate  the  same  peak  because  of 
the  second  reason 
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Figure  2.14:  The  symbol  streams  that 
generate  the  same  peak  because  of 
the  third  reason 
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Table  2.6:  The  required  memory  and  the  suppression  performance  of  look-up  tables 


Table  size 

Original 
Memory  (bytes) 

Reduced 
Memory  (bytes) 

Original 
Peak  (dB) 

Suppressed 
Peak  +  loss  (dB) 

N=4 

2*4- 
=  512 

2*4" /16 
=  32 

4.75  (a=  0.20) 
3.09  (a=  0.35) 

3.49  (a=0.20) 
2.45  (a=  0.35) 

N=5 

2*4^ 
=  2048 

2*4'  /32 
=  64 

4.75 
3.09 

3.24 
2.24 

N=6 

2*4'^ 
=  8192 

2*4V16 
=  512 

4.75 
3.09 

3.02 
2.10 

N=7 

2*4' 
=  32768 

2*4' /32 
=  1024 

4.75 
3.09 

2.90 
2.07 

Furthermore,  through  new  table  size  reduction  schemes,  the  table  size  of  an  N=4, 
6  (even)  look-up  table  can  be  reduced  by  the  factor  of  1/16  and  that  of  an  N=5,  7  (odd) 
look-up  table  can  be  reduced  by  the  factor  of  1/32.  This  possibility  can  be  confirmed  by 
detecting  2  or  4  common  ak  values  in  the  1/8  size  reduced  look-up  table.  The  1/16 
reduction  is  possible  because  a  clockwise  symbol  stream  in  the  signal  constellation  and 
the  corresponding  counter  clockwise  symbol  stream  have  the  same  peak.  Let's  call  this 
the  second  reason  of  table  size  reduction.  For  example,  the  symbol  streams  [5616]  and 
[5414]  have  the  same  peak,  as  can  be  seen  in  Figure  2.13.  Then,  the  16  symbol  streams 
[5616],  [6727],  ...  ,[4505]  and  [5414],  [6525],  ...  ,[4303]  have  all  the  same  peak  value 
because  of  the  first  and  the  second  reason.  The  1/32  reduction  is  possible  because  the 
flipped  symbol  stream  has  the  same  peak  as  the  original  one.  Let's  call  this  the  third 
reason  of  table  size  reduction.  It  should  be  noted  that  1/32  size  reduction  is  only  applied 
to  an  N=5,  7  (odd)  look-up  table  because  the  flipped  symbol  stream  of  an  even  number 
does  not  have  the  same  peak  value  in  n/4  QPSK.  For  example,  the  data  streams  [56167] 
and  [76165]  have  the  same  peak,  as  can  be  seen  in  Figure  2.14.  So  the  32  symbol  streams 
[56167],  [67270],  ...  and  [54143],  [65254],  ...  and  [76165],  [07276],  ...  and  [70501], 
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[01612],  ...  all  have  the  same  peak  value  because  of  the  first,  the  second  and  the  third 
reasons.  Table  2.6  shows  the  required  memory  and  suppression  capability  for  n/4  QPSK 
signal  with  a  RRC  pulse  (a=  0.20,  0.35). 

2.5.2  Look-up  Table  Generation 

There  is  a  technical  tip  for  implementing  a  look-up  table.  The  amplitude 
adjustment  factor  ak  depends  on  both  the  value  of  the  left  side  peak  and  the  right  side 
peak  when  N  is  even  number.  But,  the  previous  look-up  table  of  Section  2.3  was 
implemented  by  calculating  ak  only  from  the  right  side  peak  value.  So  the  APS 
suppression  result  of  the  look-up  table  was  not  stable  because  the  left  side  peak  was  not 
considered.  Therefore,  two  ak  which  is  left  side  ak  (kk)  and  right  side  ak  (rak)  are  needed 
to  be  considered  and  the  average  of  them  should  be  used  for  a  new  ak.  The  new  ak  is 
defined  as  (lak.i+rak)/2  where  N  is  an  even  number.  Then,  the  stable  PAPR  suppression 
can  be  achieved.  This  problem  does  not  happen  for  the  look-up  table  with  N  of  odd 
number  because  there  is  only  one  ak  at  the  middle  of  the  symbol  stream  instead  of  two  ak. 

The  look-up  table  is  made  from  the  pre-calculation  of  the  APS  algorithm.  If 
proper  parameters  P  and  y  are  used  to  implement  a  look-up  table,  then  more  PAPR 
suppression  can  be  achieved  as  was  done  in  the  algorithm  previously.  The  computational 
complexity  of  the  algorithm  is  not  a  problem  any  more  because  only  one  time 
computation  is  required  for  the  look-up  table  method,  while  the  algorithm  method 
requires  the  computation  at  every  time  of  signal  transmission.  Therefore,  25  samples  and 
50  iterations  instead  of  9  samples  and  4  iterations  can  be  considered  for  the  APS 
algorithm  for  a  look-up  table.  However,  the  algorithm  with  15  samples  and  30  iterations 
will  be  used  to  save  the  simulation  time.  The  numbers  are  good  enough  since  there  is  no 
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Table  2.7:  Maximally  suppressed  SP+Ti(dB)  for  each  number  of  iterations  and  samples 


1 5  iterations 

30  iterations 

50  iterations 

RRC  (a=.20) 
Original  peak 
4.75  dB 

9  samples 

2.43 

2.43 

2.44 

15  samples 

2.38 

2.37 

2.37 

25  samples 

2.38 

2.35 

2.34 

RRC  (a=.35) 
Original  peak 
3.09  dB 

9  samples 

1.97 

1.96 

1.96 

1 5  samples 

1.90 

1.89 

1.88 

25  samples 

1.90 

1.87 

1.86 

big  performance  difference  as  in  Table  2.7.  A  look-up  table  is  generated  through 
following  procedures.  First,  find  proper  values  of  P  and  y  for  the  algorithm  with  15 
samples  and  30  iterations  to  suppress  the  peak  maximally.  Then  find  the  RMS  values  of 
ak,  that  is  a  ,  for  each  iteration  to  normalize  the  value  of  ak.  Finally,  make  look-up  tables 
by  getting  ak  values  through  the  algorithm  for  each  data  entry  of  symbol  streams. 

2.5.3  The  Performance  of  the  Look-up  Table  and  Hybrid  Methods 

To  suppress  the  peak  maximally,  the  combination  method  of  the  look-up  table 
plus  the  APS  algorithm  can  be  considered.  The  combination  method  suppresses  the  peak 
maximally  by  the  algorithm  method  after  the  look-up  table  method.  The  PAPR 
suppression  performances  of  the  look-up  table  method  and  the  combination  method  with 
1  or  2  iterations  are  compared  in  terms  of  PS+t)  as  shown  in  Table  2.8. 

Let's  compare  the  performance  of  the  look-up  table  method  and  the  combination 
method  with  the  previous  APS  algorithm  method.  The  algorithm  method  is  much  better 
than  the  look-up  table  method  for  RRC  (a=.20).  The  algorithm  method  of  one  iteration 
can  suppress  the  PAPR  from  4.75dB  to  2.76dB,  while  the  look-up  table  method  can 
suppress  it  to  2.90-3.49dB,  as  shown  in  Table  2.8  and  Table  2.9.  The  combination 
method  for  RRC  (a=.20)  is  just  a  little  bit  better  than  the  APS  algorithm  method.  These 
are  because  the  peak  signal  of  RRC  (a=.20)  needs  an  N=13  look-up  table,  so  the  N==7 
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Table  2.8:  The  PAPR  suppression  performance  of  look-up  table  or  combination  method 


Pulse  shape 

Original  Peak 

1  aOIC  bIZC 

jr^T\  lor  loOK-up  laDie  dnu  u, i,z  ueraiioii 

RRC 
(a=.20) 

4.75  dB 

M— /I 
IN— 4 

1  AO    0  71  9 
j.Hy  ,  Z.  /  1  ,  Z.J  J 

IN— J 

J.Z'I  ,  Z.  /  /  ,  Z.JO 

IN — O 

3  09   9  69   9  48 

N=7 

2.90  ,  2.58  ,  2.45 

RRC 
(a=.35) 

3.09  dB 

N=4 

2.45  ,2.15  ,  1.98 

N=5 

2.24,2.06,  1.97 

N=6 

2.10, 1.98,  1.95 

N=7 

2.07, 1.94,  1.94 

Table  2.9:  The  performance  of  APS  algorithm 


Pulse  shape 

Original  peak 

SP  +  Tj  for  the  algorithm  with  1,  2,  3  ,4  ,15  iterations 

RRC  (a=.20) 

4.75  dB 

2.76  ,  2.60  ,  2.55  ,  2.54  ,  2.47 

RRC  (a=.35) 

3.09  dB 

2.41  ,2.19,2.13  ,2.06,  1.94 

look-up  table  is  not  big  enough.  However,  the  peak  signal  of  RRC  (a=.35)  needs  an  N=9 
look-up  table  so  the  N=7  look-up  table  is  close  enough.  Therefore,  the  PAPR  suppression 
performance  of  the  N=7  look-up  table  for  RRC  (a=.35)  is  as  good  as  that  of  the  APS 
algorithm  with  4  iterations.  The  N=7  look-up  table  requires  only  1  Kbytes  memory,  so 
the  look-up  table  method  is  very  practical  for  RRC  (a=.35). 

Therefore,  the  look-up  table  method  will  be  adopted  for  the  PAPR  suppression  of 
a  relatively  small  windowed  (<10Ts)  RRC  pulse  signal  and  the  algorithm  method  will  be 
adopted  for  a  relatively  large  windowed  (>10Ts)  RRC  pulse  signal.  In  other  words,  the 
look-up  table  method  will  be  adopted  for  the  RRC  pulse  signal  of  relatively  high  roll-off 
factor  (a>.25)  and  the  algorithm  method  will  be  adopted  for  the  RRC  pulse  signal  of 
relatively  low  roll-off  factor  (a<.25).  The  combination  method  will  be  adopted  when 
maximal  PAPR  suppression  is  necessary.  However,  these  are  subject  to  change  according 
to  practical  factors  such  as  complexity  and  implementation  cost. 
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2.6  Summary 

It  is  known  that  7i/4  QPSK  is  one  of  the  most  effective  modulations  in  terms  of 
power  efficiency  and  bandwidth  efficiency.  The  power  efficiency  of  n/4  QPSK  system 
can  be  increased  by  suppressing  its  PAPR  because  the  efficiency  of  a  linear  power 
amplifier  increases  when  the  PAPR  decreases.  The  APS  algorithm  is  introduced  to 
suppress  PAPR  without  loosing  any  bandwidth  efficiency.  However,  the  APS  algorithm 
has  high  computational  complexity,  which  makes  it  difficult  to  implement.  The  problem 
is  solved  by  the  complexity  reduction  method  in  the  algorithm  and  the  look-up  table 
method.  The  hybrid  method  of  the  look-up  table  and  the  APS  algorithm  is  also  adopted  to 
suppress  the  PAPR  more  effectively.  Hence,  the  PAPR  suppressed  n/4  QPSK,  which  is  a 
more  effective  modulation  than  the  original  n/4  QPSK,  is  introduced. 


CHAPTER  3 
PAPR  SUPPRESSION  IN  OFDM 

In  this  chapter,  orthogonal  frequency  division  multiplexing  (OFDM)  will  be 
introduced.  OFDM  is  a  very  popular  modulation  for  high  data  rate  transmission,  but  has 
serious  PAPR.  So  the  PAPR  problem  of  OFDM  will  be  discussed  and  the  previous  PAPR 
suppression  schemes  will  be  investigated.  The  major  problem  of  the  previous  schemes  is 
bandwidth  efficiency  degradation.  Thus,  a  new  bandwidth  efficient  PAPR  suppression 
schemes  will  be  proposed  and  the  performance  of  the  new  schemes  will  be  compared 
with  that  of  previous  schemes. 

3.1.  Introduction  to  OFDM 
OFDM  has  been  found  to  be  an  effective  data  transmission  scheme  in  fading  and 
multipath  channels  because  it  reduces  inter-symbol  interference  (ISI)  and  makes 
equalization  simple.  The  concept  of  using  parallel  data  transmission  and  frequency 
division  multiplexing  (FDM)  was  published  in  the  mid  60' s  [18].  These  schemes  are 
adopted  to  avoid  the  use  of  high  speed  equalization  and  to  combat  impulsive  noise  and 
multipath  distortion  as  well  as  to  use  the  available  bandwidth  efficiently.  In  the  early 
70' s,  the  application  of  discrete  Fourier  transform  (DFT)  on  FDM  was  introduced  to 
eliminate  arrays  of  sinusoidal  generators  and  coherent  demodulation  making  the 
implementation  cost-effective  [19].  Recently,  International  Mobile  Telecommunications- 
2000  (IMT-2000)  chose  MC-CDMA,  which  is  an  OFDM  application,  as  a  future  CDMA 
standard,  so  that  the  popularity  of  OFDM  is  growing  rapidly.  In  spite  of  many 
advantages,  a  major  drawback  of  OFDM  is  that  it  has  a  high  PAPR.  The  high  peak  power 
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signal  requires  a  linear  amplifier  with  a  large  dynamic  range  to  maintain  spectral 
efficiency.  So,  the  large  peaks  in  the  signal  degrade  the  efficiency  of  a  linear  power 
amplifier  seriously  [1].  Therefore,  a  good  PAPR  suppression  scheme  is  required. 

The  idea  of  OFDM  is  to  use  parallel  data  and  frequency  division  multiplexing 
with  overlapping  sub-channels.  Parallel  data  transmission  combats  impulsive  noise  by 
spreading  the  effects  over  many  symbols.  Then,  there  is  only  slight  damage  on  many 
symbols  instead  of  serious  damage  on  a  few  adjacent  symbols.  In  Figure  3.1,  the  PSK  or 
QAM  symbol  interval  is  Ts  in  the  usual  serial  data  transmission.  However,  in  parallel  data 
transmission,  the  PSK  or  QAM  symbol  interval  is  T  (=DTs)  where  D  is  the  number  of 
subcarriers.  These  symbols  have  rectangular  pulses  with  length  T  or  Ts.  Let  the  PSK  or 
QAM  symbols  of  symbol  interval  T  be  A(n),  which  is  called  the  frequency  domain 
symbols  in  OFDM,  and  let  the  subcarriers  be  exp(j2Trnt/T)  where  n=l,2,  ...  ,D.  Then, 
there  are  D  subsignals  which  are  ai(t)=A(l)exp(j27rt/T),  a2(t)=A(2)exp(i47it/T),  ...  , 
aD(t)=A(D)exp(j2D7rt/T)  where  t  is  continuous  time  (0<t<T).  We  refer  to  Ts  as  a  data 
symbol  interval  and  T  as  an  OFDM  frame  interval.  One  OFDM  frame  contains  D  PSK  or 
QAM  symbols.  The  D  subsignals  have  equally  spaced  carrier  frequencies,  i.e.,  fi=l/T, 
f2=2/T,  . . .  ,  fD=D/T.  These  carriers  are  orthogonal  to  one  another  over  the  duration  of  an 
OFDM  frame,  since 


where  a  and  b  are  any  different  integers  between  1  and  D.  Thus,  this  frequency  division 
multiplexing  method  is  called  orthogonal  FDM  (OFDM).  The  OFDM  signal  is 
constituted  by  the  time  domain  summation  of  these  D  subsignals  for  0<t<T. 


(3.1) 
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Time  Domain  Frequency  Domain 


Figure  3.1 :  Time  and  frequency  representation  of  OFDM  symbols 
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a(0  =  -^f^a„(t)  -  -^f^A{n)QxpU2mit/T)  (3.2) 

A  discrete  time  signal  a(m)  can  be  sampled  from  the  continuous  signal  a(t)  at  t=mTs 
where  m=  1,2,  ...  ,D.  Then, 

1 

aim)  =  a(0|,^„„.,  =  ^Y,A{n)QxpU2mm/  D) .  (3.3) 

The  a(m)  is  called  the  time  domain  symbols  in  OFDM  and  can  be  easily  implemented  by 
using  the  inverse  discrete  Fourier  transform  (IDFT)  because 

1 

a{m)  =  IDFT[A{n)]  =  -j=Y^Ain)exp{j2mim/  D) .  (3.4) 

Then,  the  time  domain  OFDM  symbols  A(n)  and  the  frequency  domain  OFDM  symbols 
a(m)  are  DFT  pairs,  so 

1 

A{n)  =  DFT[a{m)]  = -j=Y^a{m)Q\p{-j27mm  I D) .  (3.5) 

Hence,  the  OFDM  signal  a(t)  can  be  easily  implemented  by  IDFT  and  a  LPF  without 
arrays  of  sinusoidal  generators. 

a{t)  =  LPF[a{m)]  =  LPF{IDFT[A{n)]) .  (3.6) 
According  to  (3.3),  a  discrete  signal  a(m)  is  the  sampled  signal  of  an  OFDM  signal  a(t)  at 
time  interval  Ts.  Therefore,  if  the  LPF  is  ideal,  then  the  OFDM  signal  a(t)  can  be 
successfully  generated  in  (3.6)  according  to  Nyquist  sampling  theorem. 

The  transmitting  process  of  this  DFT-based  OFDM  is  as  follow  [20].  First,  h- 
length  binary  data  streams  are  modulated  onto  the  complex  symbols  A(n)  n=l,2,  ...  ,D 
through  a  PSK  or  a  QAM  modulation  scheme,  which  has  h  bits  per  symbol.  The  D 
frequency  domain  symbols  constitute  a  symbol  frame.  Second,  the  IDFT  of  A(n)  is  taken 
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Figure  3.2:  FFT  based  OFDM  system  [20] 


and  a  new  symbol  frame  which  contains  the  time  domain  symbols  a(m)  m=l,2,  ...  ,D  are 
sent.  Even  though  several  a(m)  are  completely  lost,  the  A(n)  can  be  recovered  from  the 
rest  of  a(m)  with  some  BER  degradation  since  each  a(m)  contains  the  information  of  all 
A(n)  as  (3.4).  Therefore,  OFDM  transmits  hxD  bits  simultaneously  during  one  OFDM 
frame  via  D  subchannels.  For  instance,  if  32-QPSK-OFDM  (D=32,  h=2)  is  used,  the  32 
point  IDFT  of  32  QPSK  symbols,  A(n),  is  taken  and  the  32  IDFT  symbols,  a(m),  are  sent. 
Third,  parallel  data  is  converted  to  serial  data  as  in  Figure  3.2,  and  then  guard  intervals 
are  made  by  placing  a  time  interval,  Tg,  between  OFDM  symbol  frames  as  in  Figure  3.3. 
This  prevents  ISI  caused  by  multipath  distortion  in  a  frequency  selective  fading  channel. 
The  guard  interval  makes  equalization  simple  because  of  the  low  ISI.  The  guard  interval, 
Tg,  is  usually  less  than  T/4  and  larger  than  the  channel  impulse  response  or  the  multipath 
delay.  Fourth,  let  a(k,m)  k=  1,2,3,  ...  be  time  domain  symbols  in  the  k"'  OFDM  symbol 
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Figure  3.3:  Time  and  frequency  representation  of  OFDM  with  guard  interval  [20] 

frame.  Then,  the  discrete  signal  a{k,m),  with  Tg  is  converted  to  the  analog  signal,  a{t), 
through  a  LPF  as  follow. 

a{t)  =  LPF[a{m)]  =  LPF[Y^a{k,m-  k{T^  +  T^))] .  (3.7) 

k 

Last,  the  analog  signal  is  converted  to  a  bandpass  signal  through  an  up-converter  and 
transmitted  as  in  Figure  3.2.  In  Figure  3.2,  the  DFT  is  replaced  by  the  fast  Fourier 
transform  (FFT),  which  is  a  mathematically  equivalent  method  but  is  much  faster  to 
calculate. 

The  ideal  spectrum  of  an  OFDM  sub-channel  is  a  sine  function  because  the  time 
domain  pulse  shape  of  A(n)  is  rectangular  with  length  T  as  found  in  Figure  3.1.  Each  sub- 
channel is  spaced  orthogonally  in  frequency  such  that  the  composite  spectrum,  that  is  the 
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sum  of  individual  spectrums,  is  flat,  as  shown  in  Figure  3.1.  So  the  ideal  null-to-null 
bandwidth  (BWnuii)  is  [D+l]/T,  which  is  close  to  1/Ts  when  D  is  big  enough.  However, 
this  efficient  bandwidth  is  only  valid  when  the  LPF  is  ideal.  In  practice,  if  a  Nyquist  RRC 
pulse  shaping  LPF  that  spans  many  symbols  is  adopted,  then  the  spectrum  of  OFDM  will 
be  exactly  same  as  the  spectrum  of  the  general  PSK  or  QAM.  Both  spectrums  are  the 
same  because  the  pulse  shape  is  the  same,  the  symbol  interval  is  the  same  and  the 
symbols  are  independent  to  one  another.  Thus,  the  BWnuii  of  them  is  (l+a)D/T  or 
(l+a)/Ts  when  the  a  is  a  roll-off  factor  of  RRC  pulse. 

The  PAPR  of  OFDM  can  be  measured  in  two  ways.  The  one  is  general 
continuous  PAPR  which  is  the  PAPR  of  a  continuous  OFDM  signal  a(t)  and  the  other  is 
discrete  PAPR  which  is  the  PAPR  of  a  discrete  OFDM  signal  a(m)  before  pulse  shaping 
LPF.  The  PAPR  of  a  signal  is  measured  in  terms  of  continuous  PAPR  in  the  most  of 
papers,  but,  in  the  papers  of  selective  mapping  scheme  [24]-[26],  the  PAPR  of  a  signal  is 
measured  in  terms  of  discrete  PAPR.  The  continuous  PAPR  of  OFDM,  that  is  general 
PAPR,  is  defined  by 

p^pj^  =    Peak_Power    ^  Max{a^{t)\ 
A  verage  _  Power    Mean[a  ^  (/)] 

where  t  is  continuous  time.  The  discrete  PAPR  is  based  on  the  time  domain  symbols 
a(m),  which  is  the  sampled  value  of  a(t)  and  defined  by 

^.  Discrete    Peak    Power  Max\a^(m)] 

Discrete  _  PAPR  =  =  =  =  ^  ,  (3.9) 

Discrete  _  Average  _  Power     Mean[a'  (m)] 

where  m  is  discrete  time.  The  difference  between  the  continuous  PAPR  and  the  discrete 
PAPR  depends  on  the  digital  to  analog  (D/A)  LPF.  If  the  D/A  LPF  is  just  a  smoothing 
filter  such  as  a  butterworth  filter,  then  the  difference  will  be  minimal.  However,  the 
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OFDM  signal  will  not  be  bandwidth  efficient  because  the  smoothing  filter  is  not  designed 
for  bandwidth  efficiency.  If  an  ideal  LPF  or  a  RRC  LPF  is  adopted  to  provide  good 
bandwidth  efficiency,  then  there  will  be  a  considerable  difference  between  the  continuous 
PAPR  and  discrete  PAPR.  It  is  not  practical  to  use  an  ideal  LPF  as  a  pulse  shaping  filter 
even  though  they  assume  such  a  filter  in  many  papers.  Therefore,  in  this  chapter,  a  RRC 
LPF  (a=0.35)  is  adopted  as  a  default  pulse  shaping  D/A  LPF  for  OFDM.  Then,  the 
continuous  PAPR  of  a  peak  suppressed  OFDM  signal  is  about  3dB  higher  than  the 
discrete  PAPR  of  the  signal  and  the  continuous  PAPR  of  an  original  OFDM  signal  is 
about  IdB  higher  than  the  discrete  PAPR  of  the  signal.  Generally,  The  PAPR  of  OFDM 
signal  a(t)  is  around  10-13dB,  while  the  PAPR  of  7r/4  QPSK  is  around  3-4dB.  This  high 
PAPR  of  OFDM  decreases  the  efficiency  of  a  power  amplifier  much  more  than  the  PAPR 
of  7i/4  QPSK  does,  as  can  be  found  in  Figure  1.2.  Several  PAPR  suppression  schemes 
have  been  introduced  as  can  be  seen  in  the  next  section  and  a  new  scheme  will  be 
introduced  in  this  chapter  to  suppress  the  PAPR  more  effectively. 

3.2  Previous  Work  on  OFDM  PAPR  Suppression 
There  have  been  three  kinds  of  peak  suppression  methods  presented  in  previous 
work.  One  involves  coding  [21]-[22]  another  involves  phase  rotation  [23]-[27]  and  the 
other  involves  clipping  [28]-[32]. 

In  a  paper  written  by  Wilkinson  and  Jones  [21],  a  block  coding  method  is 
presented  which  employs  only  the  frames  that  have  the  lowest  PAPR.  For  instance,  there 
are  4  symbols  in  an  OFDM  frame  in  4-BPSK-OFDM  (D=4.  h=l).  As  a  result,  there  are 
16  (=2'')  possible  symbol  frames.  Four  of  the  frames  have  a  PAPR  of  6.0dB,  4  others 
have  a  PAPR  of  3.7dB  and  the  remaining  8  have  a  PAPR  of  2.5dB.  Their  scheme  only 
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employs  the  8  symbol  frames  that  have  a  PAPR  of  2.5dB  to  avoid  high  PAPR.  Then  the 
number  of  meaningful  data  symbols  in  an  OFDM  frame  is  only  3  (rate  3/4  code).  So  it 
reduces  the  PAPR  from  6dB  to  2.5dB  at  a  cost  of  33.3%  in  overhead.  This  scheme  also 
reduces  the  PAPR  of  8-BPSK-OFDM  (D=8,  h=l)  from  9dB  to  4.4dB  at  a  cost  of  14% 
overhead  (rate  7/8  code).  As  a  side  benefit,  this  scheme  may  provide  some  error 
correction  capability.  In  a  paper  written  by  Van  Nee  [22],  an  advanced  version  of  the 
block  coding  technique  is  presented.  This  approach  uses  complementary  codes  and 
reduces  the  PAPR  of  8-8PSK-OFDM  (D=8,  h=3)  to  3dB  at  a  cost  of  100%  overhead  (rate 
1/2  code).  It  should  be  noted  that  this  scheme  has  a  coding  gain  of  3dB  for  this  case.  If  a 
conventional  error  correction  coding  gain  is  4.77dB  (rate  1/2  code),  they  sacrifice  1.77dB 
of  coding  gain  to  allow  for  PAPR  suppression.  However,  the  achievable  code  rate 
decreases  for  an  increasing  number  of  carriers.  For  example,  the  code  rate  is  only  3/16 
(433%  overhead)  to  suppress  the  PAPR  of  32-8PSK-OFDM  to  3dB.  So  they  propose  to 
use  four  of  8-8PSK  codes  (rate  1/2  code)  in  that  case.  But  the  results  of  using  these 
several  short  codes  are  not  mentioned  clearly.  Therefore,  a  good  coding  solution  for  a 
large  number  (D>32)  of  subcarriers  has  not  yet  been  proposed. 

Another  approach  is  the  phase  rotation  method.  In  a  paper  written  by  Friese  [23], 
an  iterative  optimization  algorithm  of  generating  a  low  crest- factor  (=  low  PAPR)  in  an 
OFDM  signal  is  adopted.  An  OFDM  frame  A(n)  n=l,2,  ...  ,D  is  divided  into  several 
blocks  and  the  phase  of  the  symbols  in  each  block  are  rotated  before  the  IDFT  to  have  the 
lowest  PAPR  after  IDFT.  The  first  symbol  in  each  block  is  not  a  data  symbol  but 
specifies  the  rotated  phases  in  each  block  for  demodulation.  Then  the  overhead  is  1 00/(B- 
1)%  where  B  is  the  number  of  symbol  in  each  block.  So  the  overhead  of  this  coding 
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scheme  increases  as  the  number  of  divided  blocks  increases.  For  example,  the  32-OFDM 
divided  into  8  blocks  (B=4)  has  33%  overhead  and  the  32-OFDM  divided  into  16  blocks 
(B=2)  has  100%  overhead.  An  optimization  algorithm  is  exploited  to  find  the  best  values 
of  the  rotated  phase  to  reduce  the  PAPR.  This  scheme  can  be  used  for  any  number  of 
subcarriers  and  reduces  the  PAPR  of  32-BPSK-OFDM  to  3.3dB  at  a  cost  of  100% 
overhead.  In  a  paper  written  by  Bauml  et  al.  [24],  the  selected  mapping  (SLM)  method  is 
presented.  The  phase  rotated  A'(n)  =  A(n)expGwk),  whose  IDFT  a'(m)  has  the  lowest 
PAPR,  is  selected  and  the  a'  (m)  is  transmitted.  There  is  no  block  division  in  this  method 
so  the  phase  rotation  values  exp(jwk)  in  a  frame  are  all  the  same.  The  integer  k  is  also 
sent  as  additional  side  information  to  reconstruct  A(n)  from  A'(n)  in  the  receiver.  This 
side  information  adds  only  0.4%  overhead  and  the  scheme  reduces  the  discrete  PAPR 
from  11.4dB  to  8.6dB  in  128-16QAM-OFDM  when  a  4"'  order  SLM  is  employed.  In  a 
paper  written  by  Muller  and  Huber  [25],  the  partial  transmit  sequence  (PTS)  technique, 
which  is  an  advanced  version  of  the  SLM  method,  is  presented.  However,  the  basis  of 
these  two  methods  is  the  same.  PTS  has  larger  computational  complexity  and  higher 
overhead  than  SLM,  but  suppresses  the  PAPR  more.  It  suppresses  the  discrete  PAPR  of 
128-16QAM-OFDM  from  11.4dB  to  6.3dB  at  a  cost  of  1.6%  overhead  by  using  a  S"' 
order  PTS.  The  cost  increases  when  the  number  of  subcarriers  decreases.  According  to 
Muller  and  Huber  [26],  the  discrete  PAPR  suppression  performance  of  QPSK-OFDM  is 
almost  the  same  as  16QAM-0FDM  when  SLM  or  PTS  is  used.  If  continuous  PAPR  that 
is  the  PAPR  of  a(t)  is  considered  in  these  schemes,  then  much  less  PAPR  reduction 
performance  is  expected.  So  a  new  phase  optimization  criteria  based  on  continuous 
PAPR  is  suggested  in  a  paper  written  by  Tellambura  [27].  However,  the  PAPR 
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suppression  performance  of  PTS  in  continuous  PAPR  (11.8-^9.6dB  for  128-QPSK- 
OFDM)  is  not  attractive  even  with  the  new  optimization  criteria.  The  discrete  OFDM 
signal  is  oversampled  by  a  factor  of  8  to  measure  continuous  PAPR  in  the  paper  [27].  The 
above  PAPR  of  SLM  or  PTS  are  found  from  those  papers  [24]-[27]  where  the  probability 
that  the  real  PAPR  is  higher  than  the  measured  values  is  about  10"''. 

The  other  approach  is  clipping  or  limiting  an  OFDM  signal  to  reduce  the  PAPR 
[28]-[32].  In  the  paper  written  by  Li  and  Cimini  [28],  an  analog  clipping  and  filtering 
scheme  is  used  for  PAPR  suppression.  A  discrete  OFDM  signal  is  oversampled  by  a 
factor  of  8  to  approximate  the  analog  OFDM  signal  a(t)  and  the  amplitude  of  samples  are 
clipped  at  a  certain  level  to  remove  the  peak  signal  portion  of  OFDM  signal.  The  PAPR 
is  reduced  a  lot  but  it  causes  sidelobe  regeneration  losing  spectral  efficiency.  So  a  filter  is 
employed  to  remove  the  regenerated  sidelobes  caused  by  the  clipper,  but  the  filter  causes 
additional  PAPR  again.  Consequently,  the  PAPR  of  128-QPSK-OFDM  is  reduced  from 
13dB  to  9dB  at  a  cost  of  IdB  power  loss  at  10"^  BER  and  minimal  spectral  loss.  If  more 
spectral  loss  is  allowed,  then  more  PAPR  reducfion  must  be  achieved.  Van  Nee  and  Wild 
introduces  a  peak  windowing  instead  of  clipping  in  their  paper  [29].  A  Cosine,  Kaiser  or 
Hamming  window  is  exploited  only  where  there  is  a  peak  signal  to  remove  the  peaky  part 
of  the  analog  OFDM  signal.  The  PAPR  of  48-16QAM-OFDM  is  reduced  to  4dB  with 
0.5dB  power  loss  at  10"^  BER.  But,  almost  double  bandwidth  is  required  to  satisfy  -60dB 
adjacent  channel  interference.  However,  they  insist  that  this  spectral  loss  is  minimal  if  the 
spectral  impairment  caused  by  the  non-linearity  of  linear  amplifier  is  considered.  In  the 
papers  written  by  O'Neill  and  Lopes  [30],  Dinur  and  Wulich  [31],  and  Wulich  and 
Goldfeld  [32],  the  amplitudes  of  discrete  symbols  a(m)  are  limited  to  suppress  the  peak 
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Table  3.1 :  The  performance  comparison  of  previous  PAPR  suppression  techniques 


Suppression 
Method 

Modulation 
Used 

Original 
PAPR 
(dB) 

Reduced 
PAPR 
(dB) 

Spectral  and 
Power  Cost 

Block 
Coding  [21] 

8-BPSK-OFDM 

9.0dB 
9.0dB 

4.45dB 
3.0dB 

14%  overhead 
33%  overhead 

Complementary 
coding  [22] 

8-8PSK-OFDM 
32-8PSK-OFDM 

Any 

3.0 
3.0 

100% 
433% 

Low 
Crest-Factor 
[23] 

32-BPSK-OFDM 

Any 

6.5 
5.1 
3.3 

33% 
50% 
100% 

SLM  [24]  [26] 
("4"'  order) 

128-16QAM-OFDM 
128-QPSK-OFDM 
512-QPSK-OFDM 

11.4* 
11.4* 
11.8* 

8.6* 
8.6* 
9.3* 

0.4% 
0.8% 
0.2% 

PTS  [251  [261 
(4"'  order) 

128-16QAM-OFDM 
128-QPSK-OFDM 
512-QPSK-OFDM 

11.4* 
11.4* 
11.8* 

6.8* 
6.8* 
8.1* 

1.2% 
2.4% 
0.6% 

Clipping  and 
Filtering  [28] 

128-QPSK-OFDM 

13.0 

9.1 

2.0  dB  at  10  SER 
Small  spectral  loss 

Peak 
Windowing  [29] 

48-16QAM-OFDM 

Any 

4.0 

0.6  dBat  10"  SER 
Big  spectral  loss 

Amplitude 
Limiter 

128-QPSK-OFDM 

11.3* 
(«12.7) 

4.0* 
(«7.3) 

0.74dBat  10"'  SER 
No  Spectral  loss 

(*  Discrete  PAPR) 


signal.  But  the  performance  of  amplitude  limiter  is  not  showed  clearly  in  those  papers  to 
compare  with  other  schemes.  However,  the  simulation  results  of  this  paper  shows  that  an 
amplitude  limiter  suppresses  the  discrete  PARP  of  128-QPSK-OFDM  from  11.3dB  to 
4.0dB  at  the  cost  of  0.74dB  power  loss  at  10"^  SER.  In  summary,  the  coding  schemes 
[21]-[22]  substantially  suppress  the  PAPR  but  increase  the  spectral  cost  and  they  can't  be 
used  for  high  number  of  subcarrier.  The  phase  rotation  schemes  [24]-[27]  cost  much  less 
but  cannot  suppress  the  PAPR  as  much.  The  analog  clipping  [28]-[29]  can  suppress  the 
PAPR  a  lot  but  the  spectral  loss  is  substantial.  An  amplitude  limiter  [30]-[32]  can 
suppress  the  discrete  PAPR  effectively  but  there  is  an  additional  PAPR  increase  (about  2- 
4dB)  after  pulse  shape  filtering.  The  comparison  in  PAPR  reduction  performance  of 
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Normalized  Amplitude  of  a(m) 
Figure  3.4:  An  amplitude  distribution  of  a(m) 

these  previous  methods  is  summarized  in  Table  3.1.  The  probability  that  the  PAPR  is 
higher  than  the  measured  values  is  about  lO'"*. 

3.3.  Introduction  to  Amplitude  Limiter 
Generally,  the  discrete  OFDM  signal  a(m)  has  high  peaks  but  the  high  peaks 
happen  rarely.  So  only  8.0%  of  a(m)  has  the  discrete  PAPR  that  is  greater  than  4dB  and 
only  0.66%  of  a(m)  has  the  discrete  PAPR  that  is  greater  than  7dB  in  128-QPSK-OFDM. 
The  4dB  discrete  PAPR  corresponds  to  normalized  amplitude  of  1.58  and  the  7dB 
discrete  PAPR  corresponds  to  normalized  amplitude  of  2.24  in  Figure  3.4.  Intuitively,  we 
can  expect  that  even  if  the  signal  a(m)  is  passed  through  an  amplitude  limiter  that 
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suppresses  the  peak  values,  there  will  not  be  a  serious  distortion  because  the  most  part  of 
of  the  signal  a(m)  is  not  changed. 

It  is  true  that  clipping  an  analog  signal  causes  a  significant  out-of-band  signal 
regrowth  that  causes  bandwidth  efficiency  degradation.  However,  limiting  the  amplitudes 
of  a  discrete  sequence  does  not  change  any  of  its  spectral  characteristics.  That  is  because 
the  PSD  is  only  dependent  on  the  pulse  shape  of  a  signal,  the  symbol  interval  and  the 
correlation  between  the  time  domain  symbols  a(m).  The  pulse  shape  of  a  signal  is  the 
impulse  response  of  the  LPF  for  D/A  conversion.  It  can  be  written  in  the  following 
equation  [33]. 

PXf)  =  ^-^^^  f^Rik)cxpU27rkfTJ ,  (3.10) 

«  / 
where  s{t)  =  ^  aj{t  -  nT^ ) ,  R{k)  =  E[a„a„^,  ]  =  ^  {a„a„^^  . 

where  Ps(f)  is  a  PSD  of  a  signal  s(t),  F(f)  is  the  Fourier  transform  of  the  pulse  shape,  Ts  is 
a  symbol  interval,  R(k)  is  an  autocorrelation  of  the  data,  an  is  the  n'^  symbol,  f(t)  is  the 
pulse  shape  of  a  signal  s(t),  and  Pi  is  the  probability  of  having  the  i"^  apan+k  product.  It  is 
clear  that  the  symbol  interval  and  the  pulse  shaping  LPF  are  the  same  for  the  both 
amplitude  limited  signal  and  the  original  signal.  The  question  is  whether  the  correlation 
properties  of  both  signals  are  the  same. 

From  equation  (3.4),  it  is  known  that  a(m)  =  IDFT  [A(n)].  Let's  express  this  IDFT 

process  in  matrix  notation.  Then,  a=  ^ A,  where  a  =  [a(l),  a(2),  ...  ,  a(D)]\  A  = 
[A(l),  A(2),  ...  ,A(D)]  and  W  is  a  DxD  complex  exponential  matrix.  The  matrix  W 
turns  out  a  unitary  matrix  so  that  WW*^  =  Id,  where  W"  is  the  Hermitian  matrix  of  W  and 
Id  is  DxD  identity  matrix.  Generally,  It  can  be  assume  that  the  PSK  or  QAM  symbols 


62 

A(n)  are  uncorrelated  to  one  another  because  they  are  random  data  symbols,  random  data 
is  independent  and  the  independent  data  are  uncorrelated.  Then  the  expected  value  of 

A  a"  becomes  an  identity  matrix  {E[A  A^]  =  g'Io).  Also  the  expected  value  of  a  a"  is 

derived  to  an  identity  matrix  (E[a  5"  ]  =  E[W^  ^'^  W"]  =  WE[^  A"]W^  =  a^WW"  = 
a^Io).  Therefore,  we  can  conclude  that  the  time  domain  symbols  a(m)  are  uncorrelated  to 
one  another. 

There  may  be  a  possibility  that  the  limiting  process  may  create  a  correlation  that 
was  not  there  before.  However,  since  the  data  symbols  are  independent,  the  suppression 
of  each  symbol  is  also  independent,  and  hence  we  would  expect  the  suppressed  symbols 
to  be  uncorrelated.  Furthermore,  we  have  observed  that  simulation  results  confirm  this. 
So  the  possibility  of  creating  new  correlation  is  negligible  and  we  can  conclude  that  the 
amplitude  limited  time  domain  symbols  are  also  uncorrelated  to  one  another.  Therefore, 
the  PSD  of  an  original  signal  and  the  PSD  of  a  suppressed  signal  are  virtually  identical 
and  there  is  no  spectral  regeneration  caused  by  the  PAPR  suppression.  Therefore,  it  is 
safe  to  say  the  limiter  does  not  cause  any  bandwidth  efficiency  degradation. 

Various  amplitude  limiting  functions  that  suppress  the  high  PAPR  time  domain 
symbols  a(m)  are  investigated  for  efficient  peak  suppression.  The  PAPR  suppressed 
symbols  by  an  amplitude  limiting  function  are  denoted  as  a'(m).  Let's  find  the  best 
limiting  function  that  gives  the  least  distortion  on  A'(n)  that  are  recovered  PSK  or  QAM 
symbols  form  a'{m)  through  DFT  process  (A'(n)  =  DFT  [a'(m)]).  Intuitively,  it  is 
reasonable  that  less  amplitude  change  in  a'(m)  causes  less  distortion  in  A'(n).  So  the 
amplitudes  of  a(m)  are  suppressed  by  the  limiting  functions  without  changing  any  phase 
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Soft  Limiter  Fuction 


3  Point  Limiter  Fuction 


Input 


Input 


2  Point  Limiter  Fuction 


Simple  Limiter  Fuction 


Input 


Input 


Figure  3.5:  Various  nonlinear  amplitude  limiting  functions 


of  a(m).  Ultimately,  a  function  that  suppresses  the  PAPR  with  less  power  efficiency 
degradation  in  the  BER  or  symbol  error  rate  (SER)  plot  is  wanted.  The  various  PAPR 
suppression  functions  in  Figure  3.5  were  tested  and  a  simple  limiter  function  turned  out 
to  suppress  PAPR  more  than  the  other  functions  when  the  same  power  efficiency  loss  in 
SER  plot  is  allowed.  Therefore,  the  simple  limiter  function  will  be  adopted  for  the 
limiting  function  of  the  amplitude  limiter  for  OFDM  PAPR  suppression. 
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Another  issue  involves  whether  or  not  the  simple  limiter  is  the  best  nonlinear 
function  to  suppress  the  discrete  PAPR  and  if  not,  how  much  different  it  is.  Therefore,  the 
peak  suppression  problem  is  formulated  as  a  mathematical  optimization  problem.  The 
power  loss  is  mostly  dependent  on  the  variance  of  the  symbol  errors  which  is  E[d(nr] 
=E[|A(k,n)-A'(k,n)|  It  is  assumed  that  fixing  the  variance  of  the  symbol  errors  that  is 
mean-squared  error  (MSE)  is  nearly  equivalent  to  fixing  the  power  loss.  The  first 
constraint  of  this  problem  is  that  the  MSE  is  fixed.  The  MSE  of  10'^  corresponds  to  about 
0.7dB  power  loss  at  1 0"^  SER.  The  other  constraint  is  that  the  power  of  the  signal  is  fixed 
and  normalized  to  1 .  Then  the  following  equation  is  minimized  for  an  arbitrarily  given 
time  domain  symbols  a(m). 

Minimize  J  =Max  I  a'{rn)  I'  =  Max  \  —p=rV        -  i/(«))exp(7 — nm) 

"'     4D7t^  D 

where  I' =0.01,  fl'(w) |' =  1 .  (3.19) 

The  above  optimization  problem  was  solved  using  the  'minimax'  function  in  the 
Matlab  Toolbox  for  128-QPSK-OFDM.  For  arbitrarily  given  a(m)  m=l,2,  ...  ,128,  the 
errors  d(n)  n=l,2,  ...  ,128  that  create  the  lowest  discrete  PAPR  are  attained  and  the 
reduced  discrete  PAPR  is  measured  in  that  case.  Then  a  limiter,  which  creates  the  same 
MSE,  is  found  and  the  reduced  discrete  PAPR  by  the  limiter  is  measured  and  compared. 
Both  cases  have  the  same  power  loss  as  expected.  This  indicates  that  our  initial 
assumption  regarding  the  equivalence  between  MSE  and  power  loss  is  good.  But  this 
does  not  mean  that  the  assumption  is  always  valid.  If  the  two  reduced  PAPR  by  different 
suppression  functions  are  different  enough,  then  the  power  losses  are  not  the  same  even 
though  the  MSE  are  the  same.  By  the  way,  the  reduced  discrete  PAPR  of  the  simple 
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Figure  3.6:  The  result  of  amplitude  limiter  for  a(m),  m=l,2,  ...  ,128. 


limiter  is  10"^-10'''dB  better  than  that  of  the  computer-based  method.  This  tells  us  that  it  is 
not  likely  to  find  another  nonlinear  function  that  suppresses  the  PAPR  substantially  more 
than  the  simple  limiter. 

As  can  be  expected  from  the  simple  limiter  function  of  Figure  3.5,  the  amplitude 
limiter  just  reduces  the  amplitudes  of  a(m)  m=l,2,3,  ...  D  to  a  certain  level  if  the 
amplitudes  are  above  the  level  5  (normalized  amplitude),  otherwise  it  leaves  the  symbol 
alone.  So  the  amplitude  of  the  complex  signal  a(m)  is  limited  as  seen  in  Figure  3.6,  where 
a(m)  are  time  domain  symbols  in  an  arbitrary  OFDM  frame  of  128-QPSK-OFDM.  The 
power  of  a  amplitude  limited  symbols  ai(m)  should  be  normalized  so  that  the  average 
power  of  the  amplitude  limited  symbols  is  identical  to  that  of  the  original  symbols  a(m). 
So  the  ai(m)  is  divided  by  the  RMS  of  the  ai(m)  in  an  OFDM  frame  for  normalization.  If 
this  normalization  is  performed  frame  by  frame,  then  the  amplitude  limiter  method 
requires  some  computational  complexity  at  transmitting.  However,  if  the  RMS  of  ai(m)  is 
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pre-calculated  and  used  as  a  fixed  constant  for  all  OFDM  frames,  then  it  does  not  require 
any  computational  complexity  at  transmitting.  Therefore,  as  performed  for  n/4  QPSK  in 
Section  2.4.1,  the  long-term  RMS  value  of  ai(m)  is  pre-calculated  and  taken  for  the 
normalization  constant.  The  constant  is  a  fixed  value  for  a  certain  limiter  and  a  certain 
OFDM  modulation.  The  long-term  RMS  normalization  performs  much  better  than  the 
fame  by  frame  normalization  in  OFDM,  while  the  reverse  is  true  in  n/4  QPSK.  That  is 
because  multiplying  different  normalization  constants  to  each  OFDM  frame  causes 
additional  PAPR.  Hence,  a'(m)  is  the  amplitude  limited  and  normalized  time  domain 
symbols  of  OFDM  by  the  amplitude  limiter  method.  The  process  of  the  amplitude  limiter 
can  be  written  simply  as  follows. 

ai(m)   =a(m),  where     |a(m)|  <  5,  (3.11) 

=  6  Za(m),  where  |a(m)|  >  5. 
a'(m)  =  a,(m)  /  RMS[a,(m)].  (3.12) 
If  a  noise  free  channel  is  assumed,  the  frequency  domain  symbols  A'(n)  are 
recovered  from  the  PAPR  suppressed  time  domain  symbols  a'(m)  by  the  DFT  at  the 
receiver.  The  superimposed  signal  constellations  of  A(n)  and  A'(n)  for  128-QPSK- 
OFDM  are  shown  in  Figure  3.7.  The  original  frequency  domain  symbols  A(n)  are  general 
QPSK  symbols.  There  are  some  errors  on  the  signal  constellation  of  the  PAPR 
suppressed  frequency  domain  symbols  A'(n).  However,  the  symbols  A'(n)  does  not 
include  any  noise  term,  rather  the  noise-like  behavior  is  due  to  the  reduction  of  the  signal 
amplitude  in  the  time  domain. 
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Figure  3.7:  The  result  of  amplitude  limiter  for  A(n)  n=l,2,  ...  ,128 


3.4  Performance  of  Amplitude  Limiter 
As  mentioned  in  Section  3.3,  we  can  expect  there  will  not  be  a  serious  distortion 
by  the  amplitude  limiter  because  the  high  peaks  of  OFDM  happen  rarely  so  that  most  of 
the  signal  a(m)  is  not  changed.  As  expected,  there  is  small  SER  degradation  in  Fig.  3.8 
when  an  amplitude  limiter  suppresses  the  discrete  PAPR  of  128-QPSK-OFDM  from 
11.2dB  to  4.0dB.  The  power  efficiency  degradation  is  0.75dB  at  10"^  SER  and  0.46dB  at 
a  10'^  SER  in  this  case.  The  amplitude  limiter  that  is  used  in  this  simulation  is  called  1.5 
limiter  because  the  amplitude  limiting  level  5  is  1 .5  in  terms  of  normalized  amplitude. 

The  1.5  limiter  forces  the  normalized  amplitude  of  a(m)  to  be  1.5  if  the  amplitude 
is  higher  than  1 .5.  These  SERs  are  calculated  by  a  Monte  Carlo  simulation  and  by  a  semi- 

3  2 

analytic  simulation.  However,  both  simulation  results  are  consistent  around  10'  -10' 
SER.  The  performance  of  this  limiter  is  excellent  for  BPSK-OFDM  or  QPSK-OFDM. 
However,  it  is  not  so  excellent  for  16QAM-0FDM.  Semi-analytic  simulation  results  for 


68 


10 


10 


liJ  10 

CO 


10 


10 


10 


II 1 1 1  1  1 

II  1  M  1  1 

nil  1  1  1 
II 1 1 1  1  1 
II 1 1 1  1  1 
II 1 1 1  1  1 

Mil  1  1  1 

II 1 1 1  1 
II 1 1 1  1 
II 1 1 1  1 
II 1 1 1  1 
II 1 1 1  1 
II 1 II  1 

tH  H  1- 
II  1  1  1  1 
IIIII  1 
Mill  1 
IIIII  1 
IIIII  1 
Mil  1  1 

mil  1  1 
II 1 1 1  1  1 

Mill  ~ 
IIIII 
Mill  1 
IIIII 
IIIII 

IIIII 
IIIII  1 
IIIII  1 
IIIII 
nil  1 
II 1 1 1 

U-l  u 

II 1 1 1 

Mill 
Mill  1 

II 1 1 1  1  1 
nil  1  1  1 
MUM  1 

MM  1  M 

1                              1  _ 

1                              1  1 

*  *  *  Monte  Carlo 
=  =  =  -=  000  Semi-analytic    =  = 

Mil 
MM 

MM 

1  1 

tilE E  EEeEEEEEiiEEEEEEE 

1  1 

X^SufDpressed  

1  1  1  1  III        1  1 
1  1  IIIII        1  1 

1  1  IIIII        1  1 
1  1  1  1  III        1  1 

+  I4mi  — 1  -t 

1  1  1  mil  1 
1  1  1  mil  1 
1  1  1  mil  1 

1  1  1 1 1  III  1 
1  1  1  mil  1 

 II  1 

1  1  1 1 1  III  1 
1  1  1  mil  1 

OriginaiX  \i 

I  1  1  1  IN         1  I 

1  1  mil     1  1 
1  1  mil     1  1 

1  1  1 1  Ml        1  1 
1  1  1  1  Ml        1  1 

-1  H  1411-  -1  -1 

1  1  mil     1  1 

II  

1  1 1 1  III  1 
1  1  mil  1 
1  1  mil  1 
1  1 1 1  III  1 
J  umi  1  - 

1  1  mil  1 
1  1  mil  1 

1  1 1 1  III  1 

1  1 
1  i 

1                  1                 1  - 

4  6 

Eb/No  (dB) 


8 


10 


12 


Figure  3.8:  SER  degradation  by  the  amplitude  limiter  for  128-QPSK-OFDM  in  an 

AWGN  channel 


various  OFDM  modulations  are  shown  in  Table  3.2.  The  SER  is  calculated  for  DxlOOO 
random  data  bits  and  the  power  efficiency  losses  are  measured  by  eye  balling  in  SER 
plots.  These  results  turn  out  to  be  very  good  in  terms  of  discrete  PAPR  suppression. 

In  the  Monte  Carlo  simulation,  SER  is  a  number  of  received  error  symbols 
divided  by  a  number  of  transmitted  symbols.  The  received  symbols  are  transmitted 
symbols  plus  complex  random  AWGN.  The  Monte  Carlo  simulation  is  simple  but  takes 
long  simulation  time  because  it  requires  very  large  data  to  get  an  accurate  result.  It 
requires  at  least  1 0^  data  symbols  to  get  an  accurate  1 0"''  SER.  Therefore,  a  semi-analytic 
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Table  3.2:  The  simulation  results  of  the  amplitude  limiter  in  discrete  PAPR 


Limiter  level 
Maximum  PAPR 
Average  PAPR 
Power  loss  (at  1 0'^) 

BPSK-OFDM 

QPSK-OFDM 

16QAM-0FDM 

D  =  8 

1 .5  limiter 
9.0  ^  4.0dB 
3.9  ^  3.6dB 
0.97dB  loss 

1 .5  limiter 
8.3  ^  3.9dB 
4.3  ->  3.7dB 
0.74dB  loss 

1 .9  limiter 
8.1  ^  5.6dB 
4.3  4.2dB 
0.57dB  loss 

D  =  32 

1 .5  limiter 
10.5  ->4.0dB 
5.5  ^  4.0dB 
0.70dB  loss 

1.5  limiter 
10.1  ^4.0dB 
6.1  ->4.0dB 
0.74dB  loss 

1 .9  limiter 
9.7  ^  5.7dB 
6.0  ^  5.5dB 
0.80dB  loss 

D=128 

1.5  limiter 
10.5  ^4.0dB 
6.7  ->  4.0dB 
0.64dB  loss 

I .  5  limiter 

II.  2  ^4.0 
7.3  ^4.0 

0.75dB  loss 

1 .9  limiter 
10.4^  5.7dB 
7.3  ^  5.7dB 
0.75dB  loss 

D  =  512 

I .  5  limiter 

II.  0^4.0 
7.8-^4.0 

0.66dB  loss 

I .  5  limiter 

II.  2^4.0 
8.3  ^4.0 

0.73dB  loss 

1 .9  limiter 
11.1  ^5.7dB 
8.3  ^  5.7dB 
0.73dB  loss 

simulation  is  adopted  to  save  the  simulation  time  and  to  confirm  the  results  of  the  Monte 
Carlo  simulation.  The  semi-analytic  simulation  uses  SER  expression  of  the  adopted 
modulation,  which  is  a  function  of  Eb/No  where  Eb  is  energy  per  bit  and  No  is  noise  PSD. 
It  does  not  need  random  noise  term  because  the  SER  expression  already  includes  the 
noise  effect.  For  example,  SER  of  BPSK-OFDM  is 


SER  =  Q 


2£. 


(3.13) 


Then,  the  SER  for  a  PAPR  suppressed  frequency  domain  symbols  A'(k,n),  which  does 
not  include  noise  term,  is 


SER{k,n)  =  Q 


\A{k,n)real{A\k,n)) 


IE. 


N 


(3.14) 
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where  A(k,n)  is  n""  frequency  domain  symbol  in  OFDM  frame,  SER(k,n)  is  the  SER 
for  a  PAPR  suppressed  frequency  domain  symbol  A'(k,n),  k=l,2,  ...  and  n=l,2,  ...  ,D. 
Therefore,  the  SER  of  PAPR  suppressed  BPSK-OFDM  is  the  average  of  the  individual 
SER  of  each  PAPR  suppressed  frequency  domain  symbol,  which  is 


f 

mean 

mean 

Q 

(J 

k 

II 

V 

2E 

\A{k,n)real{A\Kn))-^ 


<•  / 


(3.15) 


On  the  other  hand,  SER  of  QPSK-OFDM  is 


SER  =  2Q 


IE. 


0' 


2E. 


(3.16) 


Then,  the  SER  of  PAPR  suppressed  QPSK-OFDM  is 


SER  ~  mean 

k 


mean 


Q\  +  Q2 


(3.17) 


where  Q\  =  Q 


{real[^I2A'(k,n)]y- 


2  2E, 


N.. 


Q2  =  Q 


iimag[^^2A'{k,n)]) 


2  2E, 


N 


The  (Tea\[yf2  A'ik,n)]f  is  a  normalized  distance  to  the  one  decision  boundary  (y  axis) 

and  (imag[  ^/2  A'(k,n)])^  is  a  normalized  distance  to  the  other  decision  boundary  (x  axis) 
where  real[]  measures  real  amplitude  of  a  symbol  and  imag[]  measures  imaginary 
amplitude  of  a  symbol.  For  16QAM-0FDM,  SER  expression  is 


SER  =  3Q 


4£, 


4 


4£. 


(3.18) 


Then,  the  SER  of  PAPR  suppressed  16QAM-0FDM  is 


SER  «  mean 

k 


mean 


Q\  +  Q2  +  Q3  +  Q4 


Q\  +  Q2  +  Q3  +  Q4 


,  (3.19) 
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where  Q\  =  Q 


\-^\Oreal[A'ik,n)-A{k,n)] 


2F 
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The  difference  between  Monte  Carlo  simulation  result  and  semi-analytic  simulation  result 
is  less  than  a  few  hundredth  of  dB  in  term  of  power  loss  at  1 0'^  SER.  Therefore,  we  can 
say  both  simulation  results  are  consistent  as  can  be  found  in  the  Figure  3.8. 

Figure  3.8  is  SER  in  an  AWGN  channel.  The  power  efficiency  degradation 
caused  by  the  amplitude  limiter  increases  while  the  SER  decreases  in  the  figure. 
However,  the  power  efficiency  degradation  is  nearly  constant  at  every  SER  in  a  fading 
channel.  A  1.4  limiter  suppresses  the  discrete  PAPR  of  128-QPSK-OFDM  from  11.2dB 
to  3.6dB  in  a  Rayleigh  fading  channel  where  fdTs=10"''.  In  the  simulation  result,  there  is 
only  about  0.3dB  power  efficiency  degradation  at  every  SER  as  in  Figure  3.9.  It  is 
interesting  that  an  error  floor  appears  for  a  relatively  fast  fading  channel  where 
fdTs=2xlO^  which  may  be  a  serious  problem  for  OFDM.  Like  in  AWGN  channel,  a 
semi-analytic  simulation  can  be  used  in  a  fading  channel.  The  SER  of  BPSK-OFDM  is 


(3.20) 


Then,  the  SER  of  PAPR  suppressed  BPSK-OFDM  is 
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Figure  3.9:  SER  degradation  by  the  amplitude  limiter  for  128-QPSK-OFDM  in  a 
Rayleigh  fading  channel  (fdTs=10"'*,  2x10"^) 


SER  =  mean 
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mean  — 
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(3.21) 


where  F  =  A(k,n)real[A'(k,n)]Eb/No.  For  QPSK-OFDM,  the  SER  is 


SER^l- 


EJN 


(3.22) 


Then,  the  SER  of  PAPR  suppressed  QPSK-OFDM  is 


SER  «  mean 
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mean 


F\ 


Vl  +  Fl 
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I  F2 
1  +  F2 


(3.23) 
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Figure  3.10:  Suppressed  discrete  PAPR  vs.  power  efficiency  loss  of  128-QPSK-OFDM 


where  Fl  =  (real[  ^  A'(k,n)])^Eb/No  and  F2  =  (imag[  ^/2  A'(k,n)])^Eb/No. 

The  PAPR  suppression  performance  of  an  amplitude  limiter  depends  on  the 
limiting  level.  When  the  amplitude  limiter  level  6  is  decreased,  more  PAPR  reduction  can 
be  achieved  but  the  power  loss  is  also  more.  So  there  is  a  trade-off  between  suppressed 
discrete  PAPR  and  power  efficiency  loss.  The  original  discrete  PAPR  of  128-QPSK- 
OFDM  is  1 1 .2dB  and  the  suppressed  discrete  PAPR  along  with  the  corresponding  power 
loss  in  an  AWGN  channel  and  a  Rayleigh  fading  channel  are  shown  in  Figure  3.10.  In  an 
AWGN  channel,  the  maximum  discrete  PAPR  is  suppressed  from  1 1.2dB  to  3.2-6.1  dB. 
The  average  discrete  PAPR  is  also  suppressed  from  7.3dB  to  the  same  level.  The 
maximum  suppressed  discrete  PAPR  and  average  suppressed  discrete  PAPR  is  exactly 
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same  for  1.9  limiter  or  the  lower  limiter  in  an  AWGN  channel.  In  a  Rayleigh  fading 
channel  where  fdTs=10"^  the  maximum  discrete  PAPR  is  suppressed  from  1 1.2dB  to  2.0- 
4.9dB  and  the  average  PAPR  is  suppressed  from  7.3dB  to  the  exactly  same  level  (2.0- 
4.9dB).  The  PAPR  suppression  performance  in  a  fading  channel  is  about  1.2dB  better 
than  in  a  Gaussian  noise  channel.  These  are  the  results  of  semi-analytic  simulations  for  a 
set  of  128x1000  random  data  symbols. 

3.5  Limiter  plus  APS  Algorithm  for  OFDM  PAPR  Suppression 
As  mentioned  in  Section  3.1,  the  PAPR  of  OFDM  can  be  measured  in  two  ways. 
The  one  is  the  continuous  PAPR  that  is  the  PAPR  of  a  continuous  OFDM  signal  a(t)  after 
D/A  LPF  and  the  other  is  the  discrete  PAPR  that  is  the  PAPR  of  a  discrete  OFDM  signal 
a(m)  before  D/A  LPF.  The  a(m)  is  called  time  domain  symbols  and  formed  by  IDFT  of 
frequency  domain  symbols  A(n).  The  PAPR  of  an  OFDM  signal,  which  is  suppressed  by 
the  amplitude  limiter,  was  measured  in  terms  of  discrete  PAPR  in  the  last  section.  The 
suppressed  discrete  PAPR  was  pretty  low  for  BPSK-OFDM  and  QPSK-OFDM  as  seen  in 
Table  3.2.  However,  the  continuous  PAPR  of  the  OFDM  signal  is  0-4dB  higher  than  the 
discrete  PAPR.  It  is  the  D/A  LPF  that  determines  this  difference  of  continuous  and 
discrete  PAPR,  which  is  called  an  additional  PAPR.  If  the  LPF  is  a  conventional 
smoothing  filter,  then  the  additional  PAPR  will  be  small,  but  the  bandwidth  of  the  signal 
will  be  wide.  If  a  bandwidth  efficient  LPF  such  as  RRC  pulse  shaping  filter  is  adopted, 
then  the  OFDM  signal  will  be  spectrally  efficient,  but  the  additional  PAPR  will  be  large. 
Generally,  bandwidth  efficiency  is  more  important  than  power  efficiency  and  can  not  be 
sacrificed  in  many  applications  of  wireless  communications.  Thus,  the  RRC  filter  is 
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adopted  as  a  default  pulse  shaping  LPF  in  this  research  and  the  APS  algorithm  is 
exploited  to  suppress  the  additional  PAPR  caused  by  the  pulse  shaping. 

In  Chapter  2,  the  APS  algorithm  was  used  to  suppress  the  peak  signal  of  7c/4 
QPSK,  which  is  constructed  by  pulse  shaping.  It  should  be  noted  that  the  discrete  PAPR 
of  n/4  QPSK  before  the  pulse  shaping  is  OdB  and  7i/4  QPSK  is  a  constant  envelop 
modulation  when  a  rectangular  pulse  shape  is  adopted.  So  the  PAPR  of  7:/4  QPSK  is  high 
in  terms  of  continuous  PAPR,  which  is  the  PAPR  of  a  continuous  signal  that  is 
constructed  by  the  overlapped  pulse  shape  such  as  RRC.  Therefore,  the  APS  algorithm 
was  exploited  to  find  the  amplitude  adjustment  factors  that  minimize  the  high  continuous 
PAPR  of  n/4  QPSK.  The  APS  algorithm  scheme  includes  pre-pulse  shaping  in  its 
algorithm  procedure  but  the  pre-pulse  shaping  is  to  find  the  amplitude  adjustment  factors. 
All  the  algorithm  procedures  and  amplitude  adjustment  are  done  before  the  pulse  shaping 
by  a  D/A  LPF. 

The  basic  idea  of  the  amplitude  limiter  scheme  and  the  APS  algorithm  scheme  is 
the  same,  in  that  they  both  adjust  the  amplitudes  of  time  domain  symbols  before  pulse 
shaping.  The  difference  is  how  to  get  the  amplitude  adjustment  factors.  The  amplitude 
limiter  scheme  gets  them  from  the  simple  limiting  function  and  the  APS  algorithm 
scheme  gets  them  from  the  complex  adaptive  algorithm.  The  major  ftmction  of  the 
amplitude  limiter  scheme  is  to  suppress  the  discrete  PAPR  of  a  signal  and  that  of  the  APS 
algorithm  scheme  is  to  suppress  the  additional  PAPR,  which  is  expected  after  pulse 
shaping.  The  amplitude  limiter  scheme  can  not  suppress  the  PAPR  in  terms  of  continuous 
PAPR  because  the  scheme  does  not  consider  the  pulse  shape  of  a  signal.  The  APS 
algorithm  itself  can  suppress  both  the  discrete  PAPR  and  the  additional  PAPR,  but  it 
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takes  many  iterations  to  get  a  good  result.  The  complexity  of  the  amplitude  limiter  is 
much  smaller  than  the  single  iteration  of  the  APS  algorithm.  On  the  other  hand,  the 
suppressed  symbols  by  the  amplitude  limiter  are  much  closer  to  the  suppressed  symbols 
by  the  APS  algorithm  than  the  original  symbols.  Thus,  a  new  hybrid  scheme,  which  is 
called  limiter  plus  APS,  is  suggested  to  reduce  the  number  of  iterations  in  the  APS 
algorithm.  The  scheme  suppresses  the  amplitude  of  time  domain  symbols  by  the 
amplitude  limiter,  then  apply  the  APS  algorithm  to  the  amplitude  limited  symbols  before 
pulse  shaping  by  a  D/A  LPF.  The  spectrum  of  the  suppressed  signal  by  the  scheme  is  not 
changed  like  the  amplitude  limiter  scheme  because  the  D/A  LPF,  the  symbol  interval  and 
the  correlation  of  the  time  domain  symbols  are  not  changed  in  (3.10). 

A  modification  is  required  to  adopt  the  APS  algorithm  in  OFDM.  That  is  a 
normalization  procedure.  The  power  of  7i/4  QPSK  symbols  are  all  the  same  because  the 
symbols  contains  their  information  in  their  phase.  So  if  the  amplitude  adjustment  factors 
are  normalized,  then  the  total  signal  power  is  normalized  as  seen  in  Step  5  of  the  APS 
algorithm  in  Section  2.2.  Form  (2.7), 


However,  the  powers  of  time  domain  symbols  of  OFDM  deviate  widely  because  of 
IDFT  process.  So  normalizing  amplitude  adjustment  factors  does  not  normalize  the  total 
signal  power.  Therefore,  the  amplitude  adjusted  time  domain  symbols  should  be 
normalized.  Then, 


a'(m)  =  o,.  (m)  I    ,    where      =   —  ^     {mf  ,     (w)  =  a{m)a^ .  (3.25) 


(3.24) 
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The  ak  is  the  ampHtude  adjustment  factors  for  the  time  domain  symbols  a(m)  and  the 
as(m)  is  the  peak  suppressed  symbols  by  an  iteration  of  the  APS  algorithm.  The 
normalized  PAPR  suppressed  symbols  are  calculated  by  multiplying  the  reciprocal  of  the 
RMS  of  the  peak  suppressed  symbols  to  the  peak  suppressed  symbols  as(m).  Then,  the 
power  of  the  original  symbols  a(m)  and  the  power  of  the  normalized  peak  suppressed 
symbols  a'(m)  will  be  the  same.  However,  the  symbol  power  in  each  OFDM  frame  is  not 
exactly  normalized  because  the  RMS  is  the  long-term  RMS,  which  is  pre-calculated  and 
fixed  for  each  case.  The  long-term  RMS  normalize  the  power  of  long-term  time  domain 
symbols,  while  the  RMS  of  the  symbols  in  an  OFDM  frame  normalize  the  power  of  the 
symbols  frame  by  frame.  However,  the  RMS  calculation  in  each  frame  requires 
additional  complexity  and  the  long-term  RMS  increases  the  PAPR  suppression 
performance  as  they  do  for  the  amplitude  limiter  scheme  in  Section  3.3.  Thus,  the 
normalization  factor,  which  is  the  long-term  RMS  of  time  domain  symbols,  is  pre- 
calculated  for  each  iteration  and  used  for  the  normalized  procedure  in  the  limiter  plus 
APS  scheme. 

Table  3.3  shows  the  PAPR  suppression  performance  for  128-QPSK-OFDM  with 
RRC  (roll-off  0.35)  when  IdB  power  efficiency  loss  at  10"^  SER  is  allowed  at  the 
receiver.  This  semi-analytic  simulation  used  (128x1000)  QPSK  data  symbols  and  15 
samples  per  time  domain  symbol  in  an  AWGN  channel  and  a  Rayleigh  fading  channel. 
The  normalization  factor  (RMS  value)  for  the  amplitude  limiter  was  calculated  from 
(128x10000)  random  data  and  the  normalization  factors  for  the  APS  algorithm  is 
calculated  from  (128x1000)  random  data.  The  limiter  plus  APS  scheme  with  4  iterations 
suppresses  the  continuous  PAPR  of  128-QPSK-OFDM  from  12.7dB  to  5.2dB  with  l.OdB 
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Table  3.3:  Suppressed  PAPR  by  limiter,  APS  algorithm  or  limiter  plus  APS 


Suppressed 
PAPR 

Original 

Limiter 
(5=L413) 

APS  alj 

gorithm 

Limiter 

plus  APS 

4 

iterations 

10 

iterations 

5=1.5  & 
1  iteration 

5=1.5  & 
4  Iterations 

Discrete 
PAPR 

11.33dB 

3.63dB 

5.56dB 

5.47dB 

4.1  OdB 

4.1  IdB 

Continuous 
PAPR 

12.69dB 

T.OOdB 

6.34dB 

5.48dB 

5.56dB 

5.17dB 

Power  Loss 
In  Gaussian 

OdB 

l.OOdB 

l.OOdB 

l.OOdB 

1  .OOdB 

l.OOdB 

Power  Loss 
In  Fading 

Ode 

0.32dB 

0.29dB 

0.30dB 

0.32dB 

0.31dB 

power  efficiency  loss  at  10"^  SER  in  an  AWGN  channel.  The  power  efficiency  loss  in  a 
fading  channel  is  about  one  third  of  the  power  efficiency  loss  in  an  AWGN  channel  for 
OFDM  PAPR  suppression 

The  complexity  of  amplitude  limiter  is  much  smaller  than  that  of  one  iteration  of 
the  APS  algorithm.  So  the  complexity  of  limiter  plus  APS  scheme  with  4  iterations  of  the 
APS  algorithm  is  close  to  that  of  the  APS  algorithm  with  4  iterations.  However,  the 
PAPR  suppression  performance  of  the  former  is  much  better  than  that  of  the  later.  Even 
the  limiter  plus  APS  with  one  iteration  shows  similar  performance  to  the  APS  algorithm 
with  1 0  iterations.  However,  the  APS  algorithm  with  very  large  number  of  iterations  will 
suppress  the  continuous  PAPR  as  well  as  the  limiter  plus  APS.  Thus,  the  limiter  plus  APS 
scheme  is  introduced  to  reduce  the  complexity  of  the  APS  algorithm  in  OFDM  PAPR 
suppression. 

3.6  Performance  Comparison  with  Previous  Work 
It  is  not  easy  to  compare  the  performance  of  the  previous  work  with  that  of  the 
amplitude  limiter  or  the  limiter  plus  APS  because  their  primary  loss  is  in  bandwidth 
efficiency  and  our  primary  loss  is  in  power  efficiency.  However,  according  to  Shannon's 
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theorem,  the  bandwidth  efficiency  loss  and  the  power  efficiency  loss  can  be  compared  by 
the  following  equation  in  an  AWGN  channel  [33]. 

C=51og,o(l  +  |),  (3.26) 

where  C  is  the  channel  capacity  (bits/Sec),  B  is  the  channel  bandwidth  (Hz),  S  is  the 
signal  power  and  N  is  the  noise  power.  The  signal  to  noise  ratio  S/N  can  be  substituted 
for  Es/No  here  because  S  =  EJTs  and  N  =  NqB  if  8=1 /Ts  is  assumed.  For  instance,  if 
0.36dB  more  power  is  required  in  terms  of  Eg/No  (from  7.33dB  to  7.69dB)  to  achieve  a 
10"^  SER  in  128-QPSK-OFDM  with  1.7  limiter,  then  it  corresponds  to  a  3.8%  overhead 
in  data  rate.  The  overhead  is  calculated  from  x  =  logio(l+10''^^^),  y  =  logio(l+10''^^^)  and 
x/y  =  1.038.  The  1.7  limiter  suppresses  the  discrete  PAPR  of  128-QPSK-OFDM  from 
1 1.2dB  to  4.9dB  as  in  Table  3.6.  There  are  several  performance  comparisons  for  BPSK- 
OFDM  in  Table  3.5  and  for  QPSK-OFDM  in  Table  3.6.  It  should  be  noted  that  the  PAPR 
is  measured  in  terms  of  a  continuous  PAPR  in  Table  3.5  and  in  terms  of  a  discrete  PAPR 
in  Table  3.6.  The  power  efficiency  loss  is  measured  in  dB  scale  at  10"^  SER  and  the 
equivalent  bandwidth  efficiency  loss  is  calculated  in  %  overhead  from  (3.26). 

The  block  coding  method  looks  better  than  the  limiter  plus  APS  scheme  with  4 
iterations  in  an  AWGN  channel;  however,  this  better  performance  is  valid  only  for  very 
small  number  of  subcarriers  (D<8).  The  performance  of  coding  schemes  for  D>8  is  not  so 
good.  Generally,  it  takes  about  100%  overhead  for  a  block  coding  scheme  or  a 
complementary  coding  scheme  to  suppress  the  continuous  PAPR  to  3dB.  Also,  there  is 
no  reasonable  coding  scheme  for  D>32.  The  low  crest  factor  scheme  is  out-performed  by 
the  limiter  plus  APS  scheme  for  D=32  as  in  Table3.5.  But,  it  is  not  possible  to  suppress 
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Table  3.5:  Performance  comparison  of  PAPR  suppression  techniques  for  BPSK-OFDM 


it  rvf 
ft  01 

carriers 

Suppression  Method 

Suppression 
Performance 

Cost  m  a 
Gaussian  Channel 

Cost  in  a 
Fading  Channel 

D  =  8 

Block  coding 
[21] 

L  J 

9.0  ^  4.5dB 
9  0^3  OdB 

1 4%  overhead 
33%  overhead 

1 4%o  overhead 
33%  overhead 

Limiter  plus  APS 
(4  iterations) 

9.0  ^  4.5dB 
9.0  ^  3.0dB 

1.40dB  <^  16% 
4.86dB  <->57% 

0.5dB 
2.2dB 

D  =  32 

Low  crest  factor 
[23] 

Anv  -4  6  SHB 
Any  5.1dB 
Any  3.3dB 

33%  overhead 
50%)  overhead 
100%  overhead 

33%  overhead 
50%)  overhead 
100%  overhead 

Limiter  plus  APS 
(4  iterations) 

12.6^6.5dB 
12.6->5.1dB 
12.6^3.5dB 

0.34dB  <^4% 
1.22dB  <^  12% 
5.0dB  ^59% 

0.05dB 
0.3dB 
1.7dB 

the  PAPR  to  very  low  level  (3.3dB  in  32-BPSK-OFDM)  by  the  limiter  plus  APS 
algorithm.  Thus,  the  limiter  plus  APS  algorithm  is  found  to  be  very  effective  when  the 
number  of  subcarriers  is  not  very  small  (D>32)  and  the  suppressed  continuous  PAPR  is 
not  very  small  (>5dB).  In  a  Rayleigh  fading  channel,  it  is  difficult  to  compare  the 
performances  numerically.  However,  the  limiter  plus  APS  looks  competitive  even  with 
the  small  number  of  subcarriers  because  the  power  loss  in  a  fading  channel  is  less  than  a 
half  of  the  power  loss  in  an  A  WON  channel. 

For  a  large  number  of  subcarriers  (D>128),  SLM  and  PTS  have  good  discrete 
PAPR  suppression  performance  with  minimal  bandwidth  efficiency  loss.  PTS  has  better 
suppression  capability  than  SLM  for  the  same  condition.  However,  an  amplitude  limiter 
out-performs  PTS  with  comparable  minimal  power  efficiency  loss  in  QPSK-OFDM  as  in 
Table  3.6.  The  performance  of  the  amplitude  limiter  is  not  so  excellent  for  16QAM- 
OFDM  but  it  is  still  competitive.  The  suppressed  PAPR  by  PTS  in  terms  of  continuous 
PAPR  is  about  2.5dB  worse  than  that  by  PTS  in  terms  of  discrete  PAPR  in  128-QPSK- 
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Table  3.6:  Performance  comparison  of  PAPR  suppression  techniques  for  QPSK-OFDM 

and  16QAM-0FDM 


#of 
carriers 

Suppression 
Method 

Suppression  Performance 
m  discrete  PAPR 

Power  loss  at  1 0  SER 
or  %  overhead 

SLM  [26] 

11.4dB->8.2-9.7dB 

<  1 .0%  overhead 

D  =  128 

PTS  [26] 

11.4dB^6.3-8.9dB 

<  3.2%  overhead 

Q 
p 

r 

1 .7  Limiter 

11.2dB^4.9dB 

0.36dBat  10"'<^3.8% 
0.1 2dB  (fadmg  channel) 

s 

K 

SLM  [25] 

ll.SdB     9.2-1 0.3dB 

<  0.2%  overhead 

D  =  512 

PTS  [25] 

11.8dB  8-lOdB 

<  0.6%  overhead 

1.7  Limiter 

11.5dB-^4.9dB 

0.33dB  at  lO""*  o  3.5% 
u.  1  luD  ^laamg  cnanneij 

SLM  [26] 

11.4dB  ^  8.2~9.7dB 

<  0.5%  overhead 

16 

D=  128 

PTS  [26] 

11.4dB  ^6.3~8.9dB 

<  1 .6%  overhead 

Q 

2.0  Limiter 

ll.ldB^6.1dB 

0.45dBat  10"' -0-3.6% 

A 
M 

SLM  [25] 

11.8dB^9.2~10.3dB 

<  0.1%  overhead 

D  =  512 

PTS  [25] 

11.8dB  ->  8-1  OdB 

<  0.3%  overhead 

2.0  Limiter 

12.0dB^6.1dB 

0.47dBat  10"' <^  3.7% 

OFDM  according  to  the  paper  written  by  Tellambura  [27].  On  the  other  hand,  the 
suppressed  PAPR  by  the  limiter  plus  APS  in  terms  of  continuous  PAPR  (5.17dB)  is  about 
1 .4  dB  worse  than  that  by  the  amplitude  limiter  without  the  APS  algorithm  in  terms  of 
discrete  PAPR  (3.77dB)  in  128-QPSK-OFDM  as  seen  in  Table  3.3.  Therefore,  a  limiter 
plus  APS  is  obviously  better  method  than  PTS  in  terms  of  continuous  PAPR,  too.  The 
clipping  and  filtering  scheme  [28]  is  obviously  worse  than  the  limiter  plus  APS  scheme. 
Because  the  former  reduces  the  continuous  PAPR  of  128-QPSK-OFDM  only  to  9.1dB 
with  2. OdB  power  loss  at  10"^  SER,  while  the  later  reduces  the  PAPR  to  5.2dB  with  IdB 
power  loss  as  in  Table  3.1  and  Table  3.3.  The  peak  windowing  scheme  [29]  is  difficult  to 
compare  the  performance  with  given  data  of  the  paper  because  the  spectral  efficiency  loss 
is  difficult  to  measure  accurately.  However,  the  scheme  is  not  so  attractive  since  it  almost 
doubles  the  bandwidth  (100%  overhead)  to  satisfy  -60dB  adjacent  chaimel  interference, 
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while  it  suppresses  the  continuous  PAPR  to  4dB.  Therefore,  the  performance  of  the 
limiter  plus  APS  scheme  is  better  than  any  other  present  PAPR  suppression  schemes  for 
BPSK-OFDM  or  QPSK-OFDM  with  large  number  of  subcarriers  (D>32)  both  in  an 
AWGN  channel  and  in  a  Rayleigh  fading  channel. 

3.7  Summary 

OFDM  is  a  very  popular  and  effective  modulation  but  its  PAPR  is  very  high. 
Therefore,  many  PAPR  suppression  schemes  are  proposed  but  most  of  them  degrade  the 
bandwidth  efficiency.  An  amplitude  limiter  can  suppress  discrete  PAPR  effectively 
without  any  bandwidth  efficiency  loss.  However,  the  suppressed  PAPR  is  increased  3- 
4dB  after  pulse  shaping.  The  additional  PAPR  caused  by  pulse  shaping  can  be  suppressed 
again  by  using  the  APS  algorithm.  Thus,  the  bandwidth  efficient  limiter  plus  APS  scheme 
is  proposed  for  PAPR  suppression  in  OFDM.  The  performance  of  the  scheme  is  excellent 
for  BPSK-OFDM  or  QPSK-OFDM  when  the  number  of  sub-carrier  is  large  (>32)  both  in 
an  AWGN  channel  and  in  a  fading  channel. 


CHAPTER  4 

PAPR  SUPPRESSION  IN  SYNCHRONOUS  MC-CDMA 
In  this  chapter,  multi-carrier  code  division  multiple  access  (MC-CDMA)  is 
introduced.  The  PAPR  problem  of  synchronous  MC-CDMA  is  described  and  the 
previous  schemes  on  the  problem  are  investigated.  The  Golay  complementary  (GC) 
sequence  and  the  limiter  plus  APS  scheme  are  adopted  to  suppress  the  PAPR  of 
synchronous  MC-CDMA  effectively.  The  performance  of  this  scheme  is  found  by  a 
semi-analytic  simulation  both  in  an  AWGN  channel  and  a  fading  channel. 

4.1  PAPR  Problem  in  Synchronous  MC-CDMA 
MC-CDMA  is  a  combination  of  OFDM  and  direct  sequence  code  division 
multiple  access  (DS-CDMA)  [34]-[37].  OFDM  has  been  found  to  be  an  effective  data 
transmission  scheme  in  multipath  fading  channels.  So  it  becomes  a  promising  technique 
for  high  data  rate  transmission  in  mobile  environments.  DS-CDMA  is  one  of  the  most 
popular  multiple  access  techniques  because  of  its  higher  user  capacity  over  conventional 
access  techniques  such  as  time  division  multiple  access  (TDMA)  and  frequency  division 
multiple  access  (FDMA).  MC-CDMA  has  advantages  of  both  techniques,  OFDM  and 
DS-CDMA,  and  is  becoming  a  popular  multiple  access  technique  for  multimedia  service 
in  mobile  radio  communications.  The  spreading  sequence  of  MC-CDMA  is  applied  to  the 
signal  in  the  frequency  domain,  while  the  spreading  sequence  of  DS-CDMA  is  applied  to 
the  signal  in  the  time  domain.  Therefore,  MC-CDMA  is  a  dual  system  of  DS-CDMA 
[34].  Actually,  the  signal  of  a  MC-CDMA  system  with  a  spreading  sequence  and  the 
signal  of  a  DS-CDMA  system  with  the  IDFT  of  the  same  spreading  sequence  are 
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identical.  Therefore,  the  bandwidths  and  the  symbol  intervals  of  both  systems  are  the 
same  for  a  given  data  rate.  However,  both  systems  are  different  because  the  demodulation 
processes  are  different. 

MC-CDMA  has  a  high  PAPR  like  OFDM  because  MC-CDMA  also  has  a  similar 
IDFT  procedure  which  creates  the  high  PAPR.  The  PAPR  of  OFDM  is  not  predictable 
because  the  input  of  each  subcarrier  is  a  different  random  data  symbol.  However,  the 
PAPR  of  MC-CDMA  depends  on  its  spreading  sequence  because  the  input  of  each 
subcarrier  is  a  chip  of  the  spreading  sequence  multiplied  by  one  data  symbol.  So  the 
PAPR  of  MC-CDMA  will  be  fixed  by  the  spreading  sequence.  In  detail,  there  are  two 
different  types  of  MC-CDMA.  One  is  synchronous  MC-CDMA  (down-link)  that 
transmits  a  signal  from  a  base  station  to  subscribers  and  the  other  is  asynchronous  MC- 
CDMA  (up-link)  that  transmits  signals  from  subscribers  to  a  base  station.  An 
asynchronous  MC-CDMA  signal  is  comprised  of  the  signal  of  just  one  user.  So  the  PAPR 
is  fixed  by  the  adopted  spreading  sequence  because  the  spreading  sequence  multiplied  by 
one  data  symbol  does  not  change  the  PAPR.  Therefore,  finding  a  set  of  good  spreading 
sequences  is  the  major  way  to  reduce  the  PAPR  in  asynchronous  MC-CDMA.  On  the 
other  hand,  a  synchronous  MC-CDMA  signal  is  the  sum  of  the  signals  of  all  active  users 
because  the  base  station  has  to  send  the  signals  for  all  active  users  simultaneously.  Then 
the  PAPR  can  be  changed  by  the  number  of  active  users,  the  data  of  the  active  users,  and 
the  spreading  sequences  that  are  assigned  for  the  active  users  because  the  PAPR  of  each 
spreading  sequence  in  a  set  is  not  the  same.  Therefore,  the  PAPR  behavior  of 
synchronous  MC-CDMA  turns  out  to  be  similar  to  that  of  OFDM  in  that  they  have 
unpredictable  high  PAPR  and  the  high  PAPR  is  fundamentally  caused  by  the  IDFT 
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Figure  4. 1 :  The  signal  generator  model  of  synchronous  MC-CDMA 
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Figure  4.2:  The  signal  generator  model  of  IDFT  based  synchronous  MC-CDMA 
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procedure.  Thus,  the  PAPR  suppression  schemes  developed  for  OFDM  can  be  applied  to 
synchronous  MC-CDMA  in  addition  to  finding  a  good  spreading  sequence. 

The  synchronous  MC-CDMA  system  model  (transmitter  part)  that  is  considered 
in  this  chapter  is  shown  in  Figure  4.1.  Instead  of  frequency  domain  data  symbols  A(n)  in 
OFDM,  the  data  dj  of  the  j""  user  is  multiplied  by  the  spreading  sequence  Sj(n)  of  the  j"" 
user,  where  n=l,  2,  ...  ,D  then  they  are  added  for  every  active  user  j.  From  (3.2),  the 
OFDM  signal  is 

a{t)  =  ^f^A{n)expU27rfj),  where  f„=f„.,+l/T.  (4.1) 

Then  the  synchronous  MC-CDMA  signal  is 
1  " 

a{()  =  —j^  ^  ^  d-Sj  (n)  expijlTtfj) ,        where  j  are  the  active  users.  (4.2) 

\D  „=i  j 

If  an  IDFT  is  used  for  OFDM  implementation,  then  the  MC-CDMA  signal  is 

a(0  =  LPFiIDFT[Y,d^S^in)]),  (4.3) 

where  LPF  is  a  pulse  shaping  filter  that  converts  a  discrete  signal  to  a  continuous  signal. 
Figure  4.1  is  the  original  MC-CDMA  signal  generator  including  large  arrays  of 
sinusoidal  generators,  which  is  difficult  to  implement  for  large  number  of  subcarriers. 
Figure  4.2  is  an  equivalent  MC-CDMA  system  based  on  IDFT.  It  is  easy  to  implement 
and  used  widely  but  the  bandwidth  of  the  signal  depends  on  the  D/A  LPF.  The  processing 
gain  N  is  the  same  as  the  number  of  IDFT  points  D  in  this  basic  MC-CDMA  system.  The 
maximum  number  of  users  U  is  the  same  as  N  or  D  for  some  spreading  sequences  such  as 
WH  sequences  but  not  for  others  such  as  Gold  sequences. 
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4.2  Previous  Work  on  Synchronous  MC-CDMA  PAPR  Suppression 
The  PAPR  suppression  schemes  of  OFDM  are  not  applicable  directly  to 
synchronous  MC-CDMA  PAPR  suppression.  The  major  reason  is  that  the  synchronous 
MC-CDMA  signal  is  the  sum  of  the  OFDM  signals  of  many  active  users.  There  are  three 
major  techniques  to  suppress  the  PAPR  in  synchronous  MC-CDMA. 

Aue  [38]  suggested  a  multi-carrier  spread  spectrum  (MC-SS)  modulation  for  a 
synchronous  MC-CDMA  substitute.  There  has  been  a  lot  of  research  on  finding 
sequences  that  generate  the  low  PAPR  in  multi-carrier  systems  [39]-[42].  Schroeder, 
Shapiro-Rudin,  and  Newman  sequences  are  representative.  Their  PAPR  are  around  5dB 
for  real  (±1)  sequence  and  around  3  dB  for  complex  (e'^)  sequence.  But  they  are  just  a 
single  sequence  not  a  set  of  sequences,  hence  they  cannot  be  used  as  a  spreading 
sequence  in  a  conventional  MC-CDMA  system  that  requires  multiple  access.  Thus,  a 
MC-SS  system  was  introduced  by  Aue  [38]  to  use  one  spreading  sequence  for  all  users. 
Each  user  has  a  different  frequency  offset  as  can  be  seen  in  Figure  4.3.  The  spread  signals 
of  each  user  do  not  interfere  with  one  another  because  they  are  orthogonal  to  one  another, 
when  the  spreading  sequence  has  good  auto-correlation  properties.  Aue  found  a  sequence 
that  has  1 .6dB  PAPR  and  6dB  dynamic  range  in  MC-SS  by  a  time-frequency  swapping 
algorithm,  while  the  Newman  or  the  Schroeder  sequence  has  about  2.6dB  PAPR.  In  other 
papers,  the  Newman  sequence  has  4.6dB  PAPR  [40]  and  the  Schroeder  sequence  has 
S.ldB  [41]  in  a  multi  carrier  signal  that  is  the  MC-SS  signal  for  a  single  user.  In  [42],  the 
PAPR  of  Newman  phases  in  the  multi  carrier  signal  is  2.6dB  because  their  definition  of 
PAPR  is  different.  Therefore,  the  PAPR  of  MC-SS  for  all  users  is  lower  than  the  MC-SS 
signal  for  a  single  user  according  to  these  papers.  Popovic  [43]  proposed  a  Huffman 
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Figure  4.3:  Spectrum  allocation  in  MC-SS  (fo=l/T) 
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sequence  that  has  3dB  PAPR  and  0.09dB  dynamic  range  and  a  GC  sequence  [44]  that  has 
6dB  PAPR  in  MC-SS.  However,  the  average  BER  of  asynchronous  MC-SS  with  the 
Huffman  sequence  is  much  worse  than  that  of  MC-SS  with  the  GC  sequence 
(Huffman~2xlO■^  Golay~4xlO"^  at  Eb/No=10dB)  [45].  The  major  drawback  of  MC-SS  is 
bandwidth  expansion.  The  required  bandwidth  of  MC-SS  is  almost  twice  the  bandwidth 
of  MC-CDMA  as  you  can  see  from  Figure  4.3.  It  is  an  open  question  as  to  how  the 
number  of  active  users  has  an  effect  on  the  PAPR  of  MC-SS.  The  PAPR  is  suspected  to 
be  increased  when  the  number  of  active  users  changes. 
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Freiberg  [46]  proposed  a  set  of  quasi-orthogonal  (QO)  spreading  sequences.  The 
set  of  QO  sequences  is  generated  by  sampHng  a  set  of  sinusoid  signals  in  different 
frequency  offset,  which  makes  the  signals  orthogonal  to  one  another,  and  taking  signs 
(±1)  of  the  samples.  Then,  the  set  of  QO  sequences  is  identical  to  the  set  of  Walsh 
Hadamard  (WH)  sequences  for  code  length  N<8.  But,  for  N>16,  the  set  of  QO  sequences 
is  slightly  different  than  the  set  of  WH  sequences.  The  QO  sequences  are  not  orthogonal 
to  one  another  any  longer  for  N>16,  while  the  WH  sequences  are  always  orthogonal  to 
one  another.  This  may  cause  considerable  multiple  access  interference  (MAI)  in  a 
synchronous  MC-CDMA  system.  The  PAPR  of  the  QO  sequences  is  about  5.5dB  while 
that  of  the  WH  sequence  is  about  7dB  when  the  processing  gain  is  32  (N=32),  the 
maximum  number  of  users  is  32  (U=32),  and  the  number  of  active  user  is  32  (K=32). 
Performance  can  be  improved  by  removing  the  limitation  imposed  by  using  only  two 
values  (±1)  in  the  QO  sequence.  The  discrete  PAPR  is  OdB  for  full  active  users 
(K=U=32)  when  the  sampled  values  of  complex  sinusoid  signal  are  used  instead  of  taking 
signs  of  them.  However,  the  resulting  MC-CDMA  signal  of  each  user  is  only  2  delta 
functions  in  time  domain,  then  the  spreading  effect  of  MC-CDMA  may  disappear.  This 
problem  can  be  reduced  by  quantizing  the  sampled  values.  However,  the  major  problem 
of  this  scheme  is  that  a  base  station  is  always  required  to  transmit  the  synchronous  MC- 
CDMA  signal  of  all  users  to  maintain  low  PAPR  regardless  of  the  number  of  active  users. 
This  will  degrade  average  power  efficiency  in  the  MC-CDMA  system. 

Ochiai  [47]  modified  the  selective  mapping  (SLM)  or  partial  transmit  sequence 
(PTS)  scheme  that  was  used  for  OFDM  PAPR  suppression  and  applied  it  to  synchronous 
MC-CDMA.  An  extended  MC-CDMA  (OFDM-CDMA)  that  transmits  several  frequency 
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symbols  in  an  OFDM  symbol  is  proposed  to  adopt  this  PAPR  suppression  scheme.  This 
scheme  can  not  be  applied  to  the  basic  MC-CDMA  system  that  transmits  one  frequency 
symbol  in  an  OFDM  symbol.  The  processing  gain  is  a  factor  of  the  number  of  IDFT 
points  (NxM=D)  in  the  extended  MC-CDMA,  while  the  processing  gain  is  the  same  as 
the  number  of  IDFT  points  (N=D)  in  the  basic  MC-CDMA.  An  orthogonal  set  of  WH 
sequences  and  an  orthogonal  set  of  GC  sequences  are  adopted  in  the  scheme.  The  WH 
sequences  generate  high  PAPR  for  a  small  number  of  active  users  and  moderate  PAPR 
for  a  large  number  of  users,  while  the  GC  sequences  generate  moderate  PAPR  for  a  small 
number  of  active  users  and  higher  PAPR  than  the  WH  sequence  for  a  large  number  of 
users.  This  scheme  uses  a  WH  sequence  or  a  GC  sequence  for  every  frequency  symbol  in 
an  OFDM  symbol  and  measures  the  PAPR  of  the  OFDM  symbol  that  is  composed  by 
every  possible  combination  of  spreading  sequences.  Then,  the  scheme  chooses  the 
combination  that  generates  the  lowest  PAPR  and  transmits  the  low  PAPR  OFDM  symbol 
with  the  side  information,  which  is  the  combination  of  spreading  sequences.  For  example, 
if  there  are  8  frequency  symbols  (M=8)  in  an  OFDM  symbol  (N=16,  D=128),  then  2* 
possible  combinations  of  WH  and  GC  sequence  can  be  assigned  to  each  frequency 
symbol.  Then  the  system  randomly  chooses  P  combinations  out  of  2*  possible 
combinations  and  measures  the  PAPR  of  the  P  OFDM  symbols  composed  by  the  P 
combinations.  Then,  the  OFDM  symbol  that  has  the  lowest  PAPR  is  chosen  and 
transmitted.  As  P  increases,  the  suppressed  PAPR  decreases  and  both  side  information 
and  system  complexity  increase.  The  PAPR  of  the  OFDM  symbol  composed  by  all  WH 
sequences,  randomly  combined  sequences,  and  all  GC  sequences  are  lldB,  ll.SdB,  and 
12.5dB,  respectively  for  M=8,  N=16,  K=16,  and  D-128  where  the  probability  that  the 
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PAPR  is  higher  than  this  number  is  about  lO""*.  Then,  the  suppressed  PAPR  by  this 
scheme  is  8.5dB  for  P=8  according  to  their  simulation  results.  In  this  case,  the  side 
information  causes  an  error  floor  around  BER  equal  to  4x10"^.  Therefore,  the  major 
drawbacks  of  this  scheme  are  that  PAPR  suppression  performance  is  weak,  the  power 
cost  is  high,  and  it  can  not  be  used  for  basic  MC-CDMA  (M=l). 

In  summary,  the  MC-SS  scheme  can  achieve  low  PAPR  but  doubles  the 
bandwidth,  the  quasi -orthogonal  sequence  scheme  shows  good  PAPR  suppression  but 
must  transmit  MC-CDMA  signal  for  all  users,  and  an  extended  MC-CDMA  with 
selective  mapping  scheme  shows  weak  PAPR  suppression  performance  and  an  error  floor 
in  the  BER  plot.  In  this  chapter,  a  limiter  plus  APS  algorithm  scheme  that  was  used  in 
OFDM  will  be  applied.  It  does  not  change  any  bandwidth  and  suppresses  the  PAPR  from 
12dB  to  4dB  (N=D=128,  K=32)  at  the  cost  of  small  power  efficiency  loss  (IdB  in  an 
AWGN  channel). 


4.3  Spreading  Sequences  for  Synchronous  MC-CDMA 
A  set  of  spreading  sequences  for  synchronous  MC-CDMA  is  desirable  to  be 
orthogonal  to  one  another  because  symbol  and  chip  synchronization  can  be  achieved 
upon  receiving  a  down-link  signal.  A  commonly  used  spreading  sequence  set  is  a  set  of 
WH  sequences.  A  set  of  GC  sequences  is  another  alternative  [47].  A  set  of  orthogonal 
WH  sequences  can  be  recursively  obtained  by 
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C*'"  is  a  WH  sequence. 


In  a  similar  manner,  a  set  of  orthogonal  GC  sequences  can  be  recursively  obtained  by 
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,  An  and  Bn  are  NxN/2  matrices. 
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,    C*^'  is  a  GC  sequence. 


WH  sequences  generate  high  PAPR  for  small  number  of  active  users  and 
moderate  PAPR  for  large  number  of  users,  while  GC  sequences  generate  lower  PAPR  for 
small  number  of  active  users  and  higher  PAPR  for  large  number  of  users  than  WH 
sequences  as  can  be  seen  in  Figure.  4.4.  The  maximum  PAPR  is  measured  by  taking  the 
highest  PAPR,  while  the  high  1 0%  PAPR  is  measured  by  taking  the  average  PAPR  of  the 
640"'  highest  peak  symbol  out  of  6400  time  domain  symbols.  Randomly  generated 
(6400x1000)  time  domain  symbols  are  used  for  this  measurement.  The  worst  PAPR 
should  be  the  criteria  in  designing  a  linear  power  amplifier  because  the  MC-CDMA 
signal  should  not  be  distorted  for  any  number  of  active  users.  But  it  is  difficult  to  say  that 
the  GC  sequence  is  a  better  choice  because  the  PAPR  of  WH  sequence  is  higher  than  that 


93 


10         15        20  25 
Number  of  Active  users 


Figure  4.4:  The  discrete  PAPR  of  WH  and  GC  sequences  in  synchronous  MC-CDMA 


of  GC  sequence  only  for  one  active  user  case.  However,  if  a  limiter  plus  APS  algorithm 
scheme  is  applied  in  a  MC-CDMA  signal,  then  the  amount  of  signal  above  the  limiting 
level  or  the  distortion  amount  caused  by  the  limiting  decides  the  performance  of  the 
scheme.  Therefore,  the  high  10%  PAPR  is  more  important  criteria  than  the  highest 
PAPR.  Then,  a  set  of  GC  sequences  has  a  better  property  than  a  set  of  WH  sequences 
because  the  worst  high  10%  PAPR  of  the  GC  set  (3.7dB)  is  much  lower  than  that  of  the 
WH  set  (6.5dB)  as  seen  in  Figure  4.4. 

If  it  is  analyzed  more,  the  amplitude  distribution  of  the  MC-CDMA  signal  with 
the  GC  set  is  approximated  by  a  Rayleigh  distribution  as  with  OFDM  signals.  So  the 
portion  of  the  signal  with  high  amplitude  is  very  small  and  the  portion  of  the  signal  with 
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Figure  4.5:  The  amplitude  distribution  of  IDFT  of  HW  and  GC  sequence 

average  amplitude  is  relatively  large  as  found  in  Figure  4.5.  However,  the  amplitude 
distribution  of  the  MC-CDMA  signals  with  the  WH  set  looks  like  exponentially 
decreasing  but  not  evenly  decreasing  function.  The  signals  in  the  large  bar  near  0  will 
cause  the  un-reducible  high  dynamic  range  and  the  signal  in  the  other  large  bar  between 
the  amplitude  1 .75  and  2.0  will  cause  big  distortion  and  large  power  efficiency  loss  in  the 
PAPR  suppression  process  by  the  limiter  plus  APS  algorithm.  Therefore,  we  prefer  to 
choose  a  set  of  GC  sequences  as  the  set  of  spreading  sequences  in  a  synchronous  MC- 
CDMA  system  for  effective  PAPR  suppression  by  the  limiter  plus  APS  algorithm.  Figure 
4.5  is  drawn  by  taking  a  histogram  of  the  amplitude  of  time  domain  symbols  of 
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synchronous  MC-CDMA  with  N=32,  K=8.  Randomly  generated  (32x50,000)  symbols 
are  used  for  this  histogram. 

4.4  Limiter  plus  APS  Algorithm  for  Synchronous  MC-CDMA 
The  PAPR  suppression  of  asynchronous  MC-CDMA  is  generally  performed  by 
finding  a  low  PAPR  generating  spreading  sequence.  However,  the  PAPR  of  a 
synchronous  MC-CDMA  signal  varies  as  the  number  of  active  users  changes,  so  just 
finding  a  low  PAPR  generating  spreading  sequence  is  not  enough  for  synchronous  MC- 
CDMA.  On  the  other  hand,  the  PAPR  characteristic  of  a  synchronous  MC-CDMA  signal 
with  GC  sequences  is  similar  to  that  of  OFDM  signals  in  that  the  distribution  of  PAPR  is 
similar  to  a  Rayleigh  distribution  so  the  portion  of  high  PAPR  signal  is  small.  Therefore, 
the  limiter  plus  APS  scheme  that  is  used  for  the  PAPR  suppression  of  OFDM  in  Chapter 
3  can  be  applied  for  the  PAPR  suppression  of  synchronous  MC-CDMA. 

The  simulation  process  of  the  PAPR  suppression  by  the  limiter  plus  APS  is  as 
follows.  The  discrete  time  domain  symbol  of  synchronous  MC-CDMA  is 

a{m)  =  IDFTC£,djSj{n)),  (4.6) 
j 

where]  is  the  active  users.  The  limiter  just  reduces  the  amplitude  of  every  a(m)  m=l,2,3, 
...  D  to  a  certain  level  5  (normalized  amplitude)  if  the  amplitude  is  above  the  level, 
otherwise  it  leaves  the  symbol  alone.  After  that,  the  PAPR  reduced  signal  ai(m)  is 
normalized  by  the  long-term  RMS  of  the  signal  ai(m).  Then  a'(m)  is  the  PAPR 
suppressed  signal  by  the  limiter  method  as  in  (3.11)  and  (3.12).  This  amplitude  limiting 
process  is  the  same  as  in  OFDM.  The  suppressed  signal  by  the  limiter  has  a  low  discrete 
PAPR.  However,  its  continuous  PAPR  will  be  increased  after  D/A  LPF.  The  PAPR  of  the 
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Figure  4.6:  An  OFDM  system  model  with  PAPR  suppression  by  limiter  plus  APS 


amplitude  limited  time  domain  symbols  is  efficiently  suppressed  in  terms  of  discrete 
PAPR  but  not  in  terms  of  continuous  PAPR.  Therefore,  the  APS  algorithm  is  adopted  to 
suppress  the  additional  PAPR  caused  by  the  pulse  shaping.  The  APS  algorithm  adjusts 
the  amplitude  of  the  amplitude  limited  time  domain  symbols  for  adaptive  continuous 
PAPR  suppression.  The  detailed  procedure  of  PAPR  suppression  by  the  APS  algorithm 
can  be  found  in  Chapter  2.  The  RRC  filter  with  roll-off  factor  0.35  and  length  8Ts  is  used 
for  pulse  shaping  in  this  simulation.  Appropriate  parameters  p  (suppression  rate)  and  y 
(limiting  level)  are  found  for  the  given  system  specs  allowing  1  dB  power  efficiency  loss 
in  an  AWGN  channel.  Then,  a  normalizing  constant,  dt,  which  is  a  long-term  RMS  value 
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of  amplitude  adjusted  symbols,  is  found  for  each  iteration  of  the  algorithm.  The 
suppressed  signal  is  multiplied  by  the  reciprocal  of  dt  after  every  iteration  to  normalize 
the  total  signal  power.  This  normalization  process  removes  block  processing  and  gives 
better  BER  performance  in  OFDM  and  synchronous  MC-CDMA,  while  it  gives  small 
BER  performance  degradation  in  tc/4  QPSK  PAPR  suppression.  The  power  of  each  MC- 
CDMA  symbol  is  not  exactly  normalized  because  the  dt  is  a  long  term  RMS  value.  But 
the  long-term  power  of  the  MC-CDMA  signal  is  normalized.  The  complexity  reduction 
schemes  of  the  APS  algorithm  in  Section  2.4  can  be  applied  here.  Then,  the 
computational  complexity  of  the  algorithm  of  4  iterations  can  be  reduces  up  to  3Mips. 
Thus,  the  PAPR  suppressed  synchronous  MC-CDMA  signal  that  is  modified  from  (4.3) 
can  be  expressed  as 


where  LPA  stands  for  limiter  plus  APS  algorithm  scheme  and  j  is  the  active  users. 

Randomly  generated  BPSK-MC-CDMA  symbols  (32x1000  or  128x1000  bits)  are 
used  in  this  semi-analytic  simulation.  The  discrete  PAPR  and  the  continuous  PAPR  of  the 
original  signal,  the  amplitude  limited  signal  and  the  fully  suppressed  signal  are  measured 
in  32-BPSK-MC-CDMA  and  128-BPSK-MC-CDMA  as  seen  in  Table  4.1  and  Table  4.2. 
The  PAPR  of  the  original  signal  is  not  the  theoretical  value  but  the  measured  value  in  this 
simulation.  The  power  loss  is  difficult  to  be  measured  because  it  deviates  highly  for  a 
demodulated  active  user.  Therefore,  we  need  to  despread  the  received  signal  for  all  active 
users  and  measure  the  power  loss  by  taking  the  average  of  the  power  loss  of  each  active 
user.  The  g-loss  is  the  power  loss  in  an  AWGN  channel  and  measured  by  using  (3.13) 
and  (3.15)  when  A(k,n)  is  substituted  by  dk.  The  f-loss  is  the  power  loss  in  a  Rayleigh 
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Table  4.1.  PAPR  suppression  of  32-MC-CDMA  (GC)  by  limiter  plus  APS 


K 

6 

nrc 

P 

Y 

Original 
PAPR  (dB) 

Discrete/ 
Continuous 

After  limiter 
PAPR  (dB) 
Discrete/ 
Continuous 

After  APS 
PAPR  (dB) 
Discrete/ 
Continuous 

g- 
loss 

(dB) 

f- 

loss 
(dB) 

4 

1.4 

.880 

1.1 

1.4 

9.0  dB 

3.5  dB 

3.3  dB 

1.0 

.25 

10.1 

6.4 

3.5 

8 

1.4 

.872 

1.2 

1.5 

9.5 

3.5 

4.1 

1.0 

.20 

10.5 

6.4 

3.8 

16 

1.4 

.866 

1.2 

1.6 

10.1 

3.6 

4.1 

1.0 

.16 

10.9 

6.6 

4.3 

32 

1.5 

.900 

1.2 

1.75 

10.9 

4.0 

4.2 

1.0 

.11 

11.8 

6.8 

5.0 

Table  4.2:  PAPR  suppression  of  128-MC-CDMA  (GC)  by  limiter  plus  APS 

K 

6 
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P 
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Original 
PAPR  (dB) 

Discrete/ 
Continuous 

After  limiter 
PAPR  (dB) 
Discrete/ 
Continuous 

After  APS 
PAPR  (dB) 
Discrete/ 
Continuous 

g- 
loss 
(dB) 

f- 

loss 
(dB) 

16 

1.4 

.865 

1.0 

1.4 

ll.OdB 

3.6  dB 

3.6  dB 

1.0 

.26 

12.0 

6.5 

3.7 

32 
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1.1 

1.5 

11.5 

3.6 
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12.1 
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1.0 
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12.7 

7.0 

4.4 
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1.5 

.896 

1.0 

1.75 

10.9 

4.0 

4.2 

1.0 

.12 

12.0 

6.9 

5.1 

fading  channel  and  measured  by  using  (3.20)  and  (3.21)  when  A(k,n)  is  substituted  by  dk. 
The  PAPR  suppression  performance  has  small  difference  for  the  number  of  active  users 
and  the  processing  gain.  The  average  PAPR  suppression  gain  is  about  7dB  in  32-BPSK- 
MC-CDMA  and  128-BPSK-MC-CDMA  with  the  GC  spreading  sequences  at  the  cost  of 
l.OdB  power  loss  in  an  AWGN  channel.  The  power  loss  in  a  Rayleigh  fading  channel  is 
almost  negligible  (less  than  0.3dB).  The  suppressed  PAPR  of  synchronous  MC-CDMA 
by  the  limiter  plus  APS  is  only  3.5-5. IdB  in  continuous  PAPR. 
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4.5  Summary 

Synchronous  MC-CDMA  transmits  a  signal  from  a  base  station  to  many  active 
subscribers.  The  signal  is  the  sum  of  many  OFDM  signals  so  it  has  very  high  PAPR.  The 
PAPR  property  of  synchronous  MC-CDMA  is  similar  to  that  of  OFDM  when  GC 
spreading  sequences  are  adopted.  Therefore,  some  of  PAPR  suppression  schemes  for 
OFDM  can  be  applied  to  synchronous  MC-CDMA  signal  with  GC  spreading  sequences. 
Thus,  GC  spreading  sequences  and  the  limiter  plus  APS  scheme  is  adopted  for 
synchronous  MC-CDMA.  The  suppressed  PAPR  by  the  limiter  plus  APS  scheme  (4-5dB) 
may  be  a  little  bit  higher  than  that  of  MC-SS  (1.6-6dB)  [38]  [45].  However,  the  limiter 
plus  APS  algorithm  does  not  lose  any  bandwidth  efficiency,  while  MC-SS  doubles  the 
bandwidth. 


CHAPTER  5 

PAPR  SUPPRESSION  IN  ASYNCHRONOUS  MC-CDMA 
In  this  chapter,  the  PAPR  property  of  asynchronous  MC-CDMA  is  described.  It 
turns  out  that  the  PAPR  of  asynchronous  MC-CDMA  depends  on  the  spreading 
sequence.  Povopic  [43]  proposed  a  Zadoff-Chu  sequence  for  a  low  PAPR  generating 
spreading  sequence.  The  characteristic  of  the  Zadoff-Chu  sequence  is  investigated  in  this 
chapter.  Then,  the  PAPR  suppressed  Zadoff-Chu  (PSZ)  sequence  is  developed  by  using 
the  APS  algorithm  to  suppress  the  PAPR  more.  The  characteristic  of  the  PSZ  sequence  is 
also  investigated  and  compared  with  that  of  other  spreading  sequences. 

5.1  PAPR  Problem  in  Asynchronous  MC-CDMA 
The  PAPR  problem  in  asynchronous  MC-CDMA  is  more  critical  than  in 
synchronous  MC-CDMA.  The  PAPR  of  the  asynchronous  MC-CDMA  signal  is 
obviously  lower  than  that  of  the  synchronous  MC-CDMA  signal  because  the 
asynchronous  signal  (uplink)  is  the  signal  of  only  one  user,  while  the  synchronous  signal 
(downlink)  is  the  sum  of  the  signals  of  all  active  users.  However,  the  power  efficiency  is 
more  important  in  asynchronous  MC-CDMA  than  in  synchronous  MC-CDMA  because 
the  asynchronous  signal  is  transmitted  from  a  portable  unit  that  has  limited  power,  while 
the  synchronous  signal  is  transmitted  from  a  base  station  that  usually  has  enough  power. 
Thus,  in  the  asynchronous  MC-CDMA  system,  lower  PAPR  is  required  to  increase  the 
efficiency  of  the  power  amplifier  so  that  it  decrease  the  battery  size  or  increase  the  talk 
time. 
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Figure  5.1 :  Asynchronous  MC-CDMA  signal  generator  model 
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Figure  5.2:  The  IDFT  based  asynchronous  MC-CDMA  signal  generator  model 


As  mentioned  in  Section  4.1,  the  PAPR  of  asynchronous  MC-CDMA  is  fixed  by 
the  assigned  spreading  sequences,  while  the  PAPR  of  synchronous  MC-CDMA  is 
changing  according  to  the  number  of  active  users.  In  detail,  the  PAPR  of  an  asynchronous 
MC-CDMA  signal  is  fixed  but  different  according  to  the  assigned  spreading  sequence.  It 
means  that  each  spreading  sequence  in  a  set  of  spreading  sequences  causes  different 
PAPR  in  an  asynchronous  MC-CDMA  system.  For  example,  the  PAPR  of  a  set  of  WH 
spreading  sequences  CN=64)  varies  roughly  from  1  IdB  to  18dB  and  that  of  a  set  of  Gold 
sequences  (N=63)  varies  from  4dB  to  lOdB.  The  worst  PAPR  caused  by  the  set  of 
sequences  decides  the  power  efficiency  of  the  linear  amplifier.  Therefore,  the  PAPR 
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problem  of  asynchronous  MC-CDMA  can  be  solved  by  finding  out  a  set  of  good 
spreading  sequences  that  causes  low  PAPR  in  the  worst  case.  There  are  also  some  other 
properties  to  be  considered  for  a  set  of  good  spreading  sequences  such  as 
corsscorrelation.  These  will  be  discussed  in  Section  5.3.  The  transmitter  of  the 
asynchronous  MC-CDMA  is  simple,  while  the  receiver  is  complicated  because  the 
received  signals  of  many  active  users  interfere  with  one  another  with  different  time 
delays  and  phase  distortions.  The  system  model  of  the  transmitter  is  as  seen  in  Figure  5.1 
and  Figure  5.2.  The  detailed  simulation  model  including  demodulation  will  be  shown  in 
Section  5.4. 

5.2  Previous  Work  on  Asynchronous  MC-CDMA  PAPR  Suppression 
There  was  a  lot  of  research  on  finding  sequences  that  generate  low  PAPR  in 
multi-carrier  system  [39]-[42]  as  mentioned  in  Section  4.2.  Schroeder,  Shapiro-Rudin, 
and  Newman  sequences  are  representative.  But  they  are  just  a  single  sequence  not  a  set  of 
sequences  that  has  low  crosscorrelation,  hence  they  cannot  be  used  as  a  spreading 
sequence  in  a  multiple  access  system. 

However,  it  is  possible  to  generate  a  set  of  orthogonal  spreading  sequences  based 
on  these  chirp-like  poly  phase  sequences.  A  set  of  GC  sequences  is  the  example.  T.  Ho 
[48]  showed  that  GC  sequences  have  less  than  6dB  PAPR  and  suggested  the  sequences  as 
a  set  spreading  sequences  for  asynchronous  (up-link)  MC-CDMA  system.  The  GC 
sequences  are  orthogonal  to  one  another  as  mentioned  in  Section  4.3.  Actually,  a  large 
subclass  GC  sequences  with  length  2"  is  called  a  Shapiro-Rudin  sequence  [49].  However, 
the  BER  performance  of  the  GC  sequence  is  not  shown  in  [48],  so  it  is  difficult  to  say 
how  good  the  spreading  sequence  is. 
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The  major  contribution  on  finding  a  set  of  good  spreading  sequences  in 
asynchronous  MC-CDMA  has  been  done  by  B.  Popovic  [43]  [45].  It  was  shown  that  the 
time  domain  crosscorrelation  function  between  the  spreading  sequences  is  not  a  proper 
interference  measure  for  the  asynchronous  MC-CDMA.  So  a  spectral  correlation  function 
was  introduced  as  a  new  criterion  for  an  interference  measurement  to  find  spreading 
sequences  with  good  BER  performance.  Well-known  Gold  sequences.  Orthogonal  Gold 
sequences,  and  Zadoff-Chu  sequences  have  good  characteristic  in  this  criterion,  so  they 
all  have  good  BER  performance.  Among  them,  the  Zadoff-Chu  sequences  [50]  have 
about  6dB  continuous  PAPR  that  corresponds  to  crest-factor  2  (=6dB  crest  factor)  and 
lower  dynamic  range  than  any  other  sequences.  According  to  our  measurement,  the 
Zadoff-Chu  sequence  (N=31)  has  OdB  discrete  PAPR  before  pulse  shaping  and  has  3.4dB 
continuous  PAPR  with  the  RRC  filter  (length  12Ts,  roll-off  0.35)  and  6.1dB  continuous 
PAPR  with  the  RRC  filter  (length  20Ts,  roll-off  0.0)  in  the  worst  case.  The  MC-CDMA 
signals  with  Gold  or  Orthogonal  Gold  sequences  have  about  1  OdB  continuous  PAPR  in 
the  worst  case.  Therefore,  he  evaluated  the  Zadoff-Chu  sequences  as  the  optimum  choice 
for  spreading  sequences  in  the  asynchronous  MC-CDMA  system  taking  crest  factor, 
dynamic  range,  and  BER  performance  into  account. 

5.3  Zadoff-Chu  Sequence 
There  are  three  major  requirements  for  a  set  of  good  spreading  sequences  of  DS- 
CDMA.  They  are  low  autocorrelation,  low  crosscorrelation,  and  the  same  total  power. 
Let  C(t)  be  a  spreading  signal,  then  usually  it  takes  the  form  of 

C{t)=  f^C,p{t-kT,),  (5.1) 

A =-00 
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where  is  a  periodic  sequence,  p(t)  is  a  pulse  shape,  and  Tc  is  a  chip  interval.  In  order 
to  achieve  code  synchronization  or  to  combat  multipath  interference  it  is  desirable  to  be 
able  to  distinguish  C(t)  from  C(t+T)  for  0<t<Ts,  where  Ts  is  a  period  of  C(t)  or  a  symbol 
interval.  Then  the  following  continuous  time  periodic  autocorrelation  function  is 
desirable  to  be  minimized  for  t^^O. 

RcA^)  =  Yl'C{t)C\t  +  T)dt.  (5.2) 

In  (5.2)  the  pulse  shape  term  is  not  a  big  factor,  then  the  following  discrete  time  periodic 
autocorrelation  function  (n^^O)  is  desirable  to  be  minimized  for  code  synchronization. 

0cci")  =  ^Y,Cik)C\k  +  n).  (5.3) 

In  order  to  minimize  MAI  in  the  DS-CDMA  system,  which  is  the  major  factor  in  BER 
performance,  it  is  desirable  that  the  following  continuous  time  periodic  crosscorrelation 
function  of  spreading  sequences  is  minimized  for  i^^j. 

Rc,cji^)  =  Y^i^C,it)C]it  +  T)dt,  (5.4) 

where  i  and  j  are  the  active  user  numbers  and  Ci  and  Q  are  the  spreading  sequences  for 
the  users.  Because  of  the  same  reasons  in  the  autocorrelation  case,  the  following  discrete 
time  periodic  crosscorrelation  function  should  be  minimized  for  i^y. 

^c,.n("'>  =  VT^C^ik)C]{k  +  n),  (5.5) 

where  C,=  [Ci(0),  Ci(l),  Ci(2),  ...  ,  Ci(N-l)]  and  C  =  [Cj(0),  Cj(l),  Cj(2),  ...  ,  Cj(N-l)] 
are  a  pair  of  spreading  sequences  among  a  set  of  the  spreading  sequences. 
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The  other  requirement  for  a  set  of  good  spreading  sequences  is  that  each 
spreading  sequence  in  a  set  of  spreading  sequences  should  have  the  same  total  power, 
although  every  sequence  element  does  not  need  to  have  the  same  power.  A  spreading 
sequence  of  DS-CDMA  is  generally  a  sequence  of  1  and  -1.  Gold  sequences  and  WH 
sequences  are  representative.  But  a  sequence  does  not  need  to  have  1  and  -1  as  its 
elements  to  be  a  spreading  sequence.  If  a  spreading  sequence  has  good  despreading 
properties,  it  doesn't  matter.  For  example,  a  Gold  sequence  is  all  1  or  -1  elements  so 
obviously  each  element  has  the  same  power  (|Ci(n)|  =  |Ci(m)|)  and  each  sequence  has  the 

same  total  power  (C^C,  =C^C-  ).  However,  each  element  does  not  need  to  have  the 

same  power.  The  Huffman  sequence  [45]  in  Section  4.2  is  an  example.  Its  elements  have 
different  powers  (1,4,8,6,-8,4,-1).  But  it  is  just  a  sequence  not  a  set  of  sequences.  So  it 
may  not  be  very  good  example  because  the  same  total  power  requirement  of  sequences  in 
a  set  of  spreading  sequences  is  not  valid  for  a  single  sequence.  Another  good  example  is 
the  IDFT  of  Gold  spreading  sequences.  The  MC-CDMA  signal  with  Gold  sequences  is 
equivalent  to  the  DS-CDMA  signal  with  the  IDFT  of  Gold  sequences.  Each  element  of 
the  IDFT  of  Gold  sequence  is  a  complex  number  and  has  different  power  but  the  total 
power  of  each  sequence  is  the  same  because  the  IDFT  process  does  not  change  the  total 
power  of  the  sequence.  The  IDFT  of  Gold  sequences  can  spread  and  despread  a  signal 
correctly,  even  though  its  elements  do  not  have  the  same  power.  Thus,  the  major 
requirements  of  asynchronous  DS-CDMA  spreading  sequences  are  the  low 
autocorrelation,  the  low  crosscorrelation,  and  the  same  total  power. 
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The  major  requirements  for  MC-CDMA  spreading  sequences  are  a  little  bit 
different.  Let  Sk(l),  Sk(2),  ...  ,Sk(D)  be  a  spreading  sequence  in  Figure  5.2,  then  its  IDFT 
is  Ck(l),  Ck(2),...A(D)  where 

1  In 
Q  (m)  =  1DFT[S,  («)]  =  7=  Z    (")  ^^PO  —  nm) .  (5.6) 

Then  the  code  synchronization  of  the  MC-CDMA  system  does  not  depend  on  the 
autocorrelation  of  the  spreading  sequence  ( 6^  ^  {n) )  but  depends  on  the  autocorrelation  of 

the  IDFT  sequence  ( 6^.  ^.  {ri) )  because  the  IDFT  sequence  is  the  time  domain  sequence. 
Because  of  the  similar  reason,  the  MAI  of  MC-CDMA  does  not  depend  on  the  time 
domain  crosscorrelation  function  of  spreading  sequence  ( 6^.  («) )  but  depends  on  the 
following  spectral  correlation  function  according  to  Popovic  [45]. 

X.A^)  =  Y^S,ik)S;ik)expU2^,r) ,         where  0<t<T,  fk=fk.,+l/T.  (5.7) 

It  has  similar  characteristic  to  the  time  domain  crosscorrelation  function  of  the  IDFT 
sequences  (i?^..^-.(r)).  So  it  can  be  substituted  for  the  discrete  crosscorrelation  of  the 

IDFT  sequences  {G-.-.{ri)).  The  same  total  power  requirement  of  the  spreading 

sequences  is  also  necessary  for  MC-CDMA  because  the  total  power  is  not  changed  by 
IDFT  process.  Thus,  the  major  requirement  of  asynchronous  MC-CDMA  is  the  low 
autocorrelation  of  IDFT  sequences,  the  low  crosscorrelation  of  IDFT  sequences,  and  the 
same  total  power  of  the  sequences. 

A  Zadoff-Chu  sequence  is  not  a  sequence  with  elements  of  1  and  -1  but  a  poly- 
phase sequence  with  unit  power  (e''^,  0<6<27:).  A  set  of  length  N  Zadoff-Chu  sequences 
can  be  obtained  by  the  following  equation. 
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2  m 

Si(k)=exp(; — {—  +  qk)),  for  N  even,  (5.8) 

N  2 

expO-^(^^  +  ^^)),  for  N  odd, 

A'  2 

where  N  is  the  processing  gain  and  a  prime  number,  i  =  1 ,2,  ...  ,N,  and  q  is  any  integer. 
The  Zadoff-Chu  sequences  are  the  special  case  of  the  generaHzed  chirp-like  sequences 
having  the  ideal  periodic  autocorrelation  function  and  the  optimal  periodic 
crosscorrelation  function  [43]  [51].  Therefore,  it  has  good  properties  as  a  spreading 
sequence  for  DS-CDMA.  A  sample  of  the  Zadoff-Chu  sequence  is  as  follows  when 
N=31,q=l,and  i=l. 

S,(k)   =[S,(0),S,(1),...,S,(N-1)]  (5.9) 
=  exp(jx[0,  0.4054,  1.0134,  1.8242,2.8376,  -2.2295,  ...  ,-0.2027]). 
The  IDFT  of  this  sequence  is 

C,(k)  =[C,(0),C,(1),...,C,(N-1)]  (5.10) 
=  exp(jx[0.5574,  0.1520,  -0.4560,  -1.2668,  -2.2802,  2.7869,  ...  ,0.7601]). 
A  wonderful  property  of  the  Zadoff-Chu  sequence  is  that  the  amplitude  of  the  IDFT 
sequence  of  Zadoff-Chu  is  all  the  same,  which  means  that  the  discrete  PAPR  of 
asynchronous  MC-CDMA  signal  is  always  OdB.  This  is  the  ideal  characteristic  of 
spreading  sequence  in  the  PAPR  suppression  of  MC-CDMA. 

The  autocorrelation  property  of  Zadoff-Chu  sequence  is  ideal.  The  discrete  time 
periodic  autocorrelation  function  6s,s(n)  is  0  for  n^^O  and  1  for  n=0  in  a  period.  The 
autocorrelation  property  of  the  IDFT  of  the  Zadoff-Chu  sequence  is  also  ideal.  Thus, 
Gc.c(n)  is  0  for  n^^O  and  1  for  n=0.  Therefore,  the  Zadoff-Chu  sequence  has  an  ideal 
autocorrelation  property  to  achieve  code  synchronization  or  to  combat  multipath 
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interference  both  in  DS-CDMA  and  MC-CDMA.  The  mean  of  Gs,s(n)  of  the  Gold 
sequences  is  0.1378  for  n^^O  and  1  for  n=0  in  a  period,  while  6c,c(n)  for  the  Gold 
sequences  is  0  for  n^^O  and  1  for  n=0.  On  the  other  hand,  the  crosscorrelation  property  of 
Zadoff-Chu  sequences  is  as  good  as  that  of  Gold  sequences.  The  discrete  time  periodic 
crosscorrelation  function  of  Zadoff-Chu  (^v/s/^^^^     0.1796  for  i^j  and  1  for  i=j.  The 

crosscorrelation  of  the  IDFT  of  Zadoff-Chu  ( 0-^   {n) )  is  also  0. 1 796  for  i^y  and  1  for  i=j. 

The  mean  of  the  crosscorrelation  of  the  Gold  sequence  is  0.1378  for  i^^j  but  the  maximum 
of  the  crosscorrelation  is  0.1967.  The  mean  crosscorrelation  of  the  IDFT  of  the  Gold 
sequence  is  0.1765  for  i^^j  and  the  maximum  of  the  crosscorrelation  is  0.1776.  Therefore, 
the  MAI  of  Zadoff-Chu  and  that  of  Gold  sequence  are  similar  both  in  DS-CDMA  and 
MC-CDMA.  Thus,  the  BER  performance  of  Zadoff-Chu  sequences  and  that  of  Gold 
sequences  are  not  expected  to  be  much  different. 

A  covariance  matrix  also  shows  that  the  Zadoff-Chu  sequences  have  good 
crosscorrelation  properties  when  there  is  no  time  delay  between  users.  Let  the  30x31 

matrix  ^3,  be  a  set  of  Zadoff-Chu  spreading  sequences  with  processing  gain  N=31  and 

the  matrix  C3,  be  the  matrix  of  its  IDFT  sequences.  The  maximum  number  of  users  is  30 
in  this  case.  Then  a  30x30  covariance  matrix  of  Zadoff-Chu  is 

,  where  aij=-L,  N=31.  (5.11) 
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Figure  5.3:  DS-CDMA  Analytical  average  BER  for  K=4  active  users 


The  diagonal  elements  are  all  1  for  any  spreading  sequence  and  the  off-diagonal  elements 
show  the  corsscorrelation  of  two  different  sequences.  The  crosscorrelations  aj j  are  all 
0.18  for  Zadoff-Chu  when  N=31,  but  they  vary  from  0.03  to  0.29  for  Gold  sequences. 
The  covariance  matrix  is  the  identity  matrix  (a;  j=0)  for  a  set  of  orthogonal  sequence  such 
as  HW  sequences.  This  tells  us  that  both  the  Zadoff-Chu  and  the  IDFT  of  the  same  are 
expected  to  have  a  pretty  good  corsscorrelation  property  when  there  is  no  time  delay, 
which  results  in  good  BER  performance  even  in  synchronous  CDMA  systems. 

Figure  5.3  shows  the  analytical  BER  performance  of  asynchronous  DS-CDMA 
for  several  spreading  sequences  including  Zadoff-Chu  sequences,  which  is  based  on  the 
standard  Gaussian  approximation.  The  crosscorrelation  parameter  defined  by  Pursley  [52] 
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Figure  5.4:  MC-CDMA  Analytical  average  BER  K=4  Active  users  [45] 


is  the  key  parameter  in  the  standard  Gaussian  approximation  based  BER.  The 
crosscorrelation  parameter  of  the  Zadoff-Chu  sequences  is  close  to  those  of  random 
sequences.  For  example,  the  parameter  ^(1)  is  0  for  random,  48  for  Zadoff-Chu,  and 
-194  for  Gold  sequences.  The  BER  is  low  when  ^(1)  is  close  to  zero.  Therefore,  the  BER 
performance  of  Zadoff-Chu  is  close  to  the  BER  of  the  random  spreading  sequence  in  a 
DS-CDMA  system.  However,  the  BER  of  Figure  5.3  is  not  very  accurate  because  the 
number  of  active  users  is  too  small  for  Gaussian  approximation  to  be  valid.  But  it  is  close 
enough  to  compare  the  BER  performance  between  different  sequences.  The  BER 
performance  in  a  MC-CDMA  system  is  found  in  [45]  and  shown  in  Figure  5.4.  The  BER 
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Figure  5.5:  Continuous  PAPR  of  Zadoff-Chu  in  asynchronous  MC-CDMA 

performance  of  Zadoff-Chu  is  better  than  that  of  Gold  sequences  with  low  PAPR  and 
ideal  autocorrelation  property  in  MC-CDMA.  Therefore,  it  can  be  concluded  that  the 
Zadoff-Chu  sequence  is  the  best  choice  for  MC-CDMA  system  in  terms  of  PAPR  and 
BER. 

The  discrete  PAPR  of  Zadoff-Chu  sequences  is  OdB,  but  the  continuous  PAPR 
after  pulse  shaping  is  not  so  low.  The  continuous  PAPR  of  length  N=3 1  Zadoff-Chu 
sequences  is  2.7-6.4dB  in  ideal  asynchronous  MC-CDMA  systems  and  2.1-4.5dB  in 
IDFT  based  asynchronous  MC-CDMA  with  RRC  pulse  roll-off  0.2.  The  continuous 
PAPR  is  different  for  each  sequence  in  a  set  of  sequences  as  seen  in  Figure  5.5.  The 
worst  PAPR  is  the  PAPR  of  the  sequences  because  a  linear  amplifier  in  a  portable  unit 
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should  be  designed  to  be  working  for  any  sequence.  This  high  continuous  PAPR  of  an 
asynchronous  MC-CDMA  signal  can  be  suppressed  by  the  APS  algorithm.  The  scheme  is 
very  simple  and  easy  to  implement. 

5.4  APS  Algorithm  for  Asynchronous  MC-CDMA 
The  Zadoff-Chu  sequence  is  the  best  known  spreading  sequence  for 
asynchronous  MC-CDMA  systems  in  terms  of  PAPR.  However,  better  spreading 
sequences  can  be  found  by  applying  the  APS  algorithm  to  the  asynchronous  MC-CDMA 
system  with  Zadoff-Chu  sequences.  The  amplitude  limiter  is  designed  to  suppress  the 
PAPR  of  the  discrete  signal,  while  the  APS  algorithm  is  designed  to  suppress  the 
continuous  PAPR  caused  by  pulse  shaping.  The  APS  algorithm  can  be  used  without  the 
limiter  where  there  is  discrete  PAPR,  but  it  requires  a  lot  of  computational  complexity 
without  the  limiter  as  mentioned  in  Chapter  3.  Therefore,  the  limiter  is  not  required  in 
asynchronous  MC-CDMA  because  the  discrete  PAPR  is  already  OdB  and  the  high 
continuous  PAPR  is  caused  by  pulse  shaping.  The  main  difference  between  the 
asynchronous  MC-CDMA  system  and  the  OFDM  or  the  synchronous  MC-CDMA  system 
is  that  the  asynchronous  MC-CDMA  system  has  fixed  amplitude  data  symbols  (Sk(n)) 
even  though  the  phase  of  the  symbol  can  be  changed  by  the  data  (dk).  From  (4.2),  a 
synchronous  MC-CDMA  signal  is 

1  " 

a{t)  = —i='^^d -S ^{n)exp{j27rfj) ,        where  j  are  active  users.  (5.12) 

There  is  only  one  user's  signal  in  asynchronous  MC-CDMA  signal,  so  that  an 
asynchronous  MC-CDMA  signal  for  user  k  is 

1  " 

ai.(')  =  ^Y.'^,S,(n)expij2;f„t),  where  k  =  1,2,3,       U.  (5.13) 

V  D  ;,=1 
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The  user  capacity  (maximum  number  of  users)  U  is  N-1  in  the  Zadoff-Chu  sequences 
while  it  is  N+2  in  the  Gold  Sequences.  The  data  dk  is  generally  a  PSK  symbol  that  only 
changes  the  phase  of  the  OFDM  symbol  ak(t).  Therefore,  for  a  given  spreading  sequence, 
Sk(n),  the  PAPR  of  a  MC-CDMA  signal  ak(t)  is  the  same  for  any  data  dk.  The  IDFT  based 
asynchronous  MC-CDMA  signal  for  user  k  can  be  written  as 

a,  it)  =  LPF{IDFT{d,  S,  («))]  =  LPF[d,  Q  («)] .  (5.14) 
The  number  of  DFT  points  D  is  always  equal  to  the  processing  gain  N  in  a  basic  MC- 
CDMA  system. 

The  continuous  PAPR  depends  on  the  IDFT  sequence  Ck(n)  and  the  LPF.  So  the 
PAPR  suppression  is  performed  by  adjusting  the  amplitude  of  Ck(n).  The  APS  algorithm 
finds  the  amplitude  adjustment  factors  adk(n)  and  multiplies  them  to  the  Ck(n).  Then,  the 
PAPR  suppressed  MC-CDMA  signal  for  user  k  is 

a\  (0  =  LPF[d,adi,in)C,{n)],  where  n=l,2,3,  ...  ,  D.  (5.15) 

It  should  be  noticed  that  the  APS  algorithm  does  not  need  to  be  performed  at  the  time  of 
transmission.  It  needs  to  be  performed  only  one  time  before  transmission  for  a  given  LPF 
because  the  Ck(n)  is  always  the  same  for  the  specific  user  k.  This  means  that  a  set  of  new 
PAPR  suppressed  spreading  sequences  is  bom  by  taking  the  DFT  of  the  PAPR 
suppressed  IDFT  sequences.  The  new  sequence  S''k(n)  can  be  expressed  as 

5;  («)  =  DFT[ad,  («)Q  («)]  =  DFT[ad,  (n)IDFT[S,  («)]] .  (5.16) 
We  call  this  new  spreading  sequence  as  the  PAPR  suppressed  Zadoff-Chu  (PSZ) 
sequence.  The  PSZ  sequence  is  slightly  different  according  to  the  adopted  pulse  shaping 
LPF  because  the  APS  algorithm  finds  proper  adjustment  factors  based  on  the  pulse 
shaped  continuous  signal  generated  by  the  LPF.  The  PSZ  sequence  also  can  be  slightly 
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varied  for  PAPR  suppression  degree  even  though  the  pulse  shaping  fiher  is  the  same.  The 
PAPR  is  suppressed  as  long  as  the  PSZ  has  almost  the  same  characteristics  (such  as  BER 
performance)  with  the  original  Zadoff-Ghu  sequence. 

The  despreading  code  in  the  receiver  is  the  original  Zadoff-Chu  sequence  in  the 
synchronous  MC-CDMA  system  because  there  is  no  specific  PAPR  suppressed 
despreading  sequence  for  a  user.  However,  the  PSZ  sequence  can  be  used  for  despreading 
because  we  know  what  specific  PSZ  sequence  is  used  in  spreading  for  the  user.  The  BER 
performance  is  better  when  the  PSZ  sequence  is  used  for  despreading  than  when  the 
original  Zadoff-Chu  sequence  is  used  for  despreading.  According  to  one  of  our  semi- 
analytical  simulation  result,  the  power  efficiency  loss  is  very  low  (-.01  dB)  when  the 
Zadoff-Chu  is  used  for  despreading  but  the  power  efficiency  loss  does  not  exist  when  the 
PSZ  is  used  for  despreading,  while  2-3dB  continuous  PAPR  suppression  is  achieved  in 
both  case. 

The  semi-analytical  simulation  model  is  shown  in  Figure  5.6.  A  discrete 
asynchronous  MC-CDMA  signal  of  user  number  k  is 

a,{m)  =  IDFT{d,S,{n))  =  d,C,{n).  (5.17) 
The  discrete  signal  of  each  user  is  distorted  by  the  independent  time  delay  and  phase 
distortion  through  the  independent  channels.  These  time  delay  and  phase  distortion 
variables  are  simplified  to  discrete  variables  (O^jik)  with  equal  probability.  Then  the 
distorted  discrete  signal  of  user  number  k  is 

a,'{m)  =  0,a,{m-r,),  (5.18) 
where  a  random  phase  variable  6k  is  exp(g7:/4),  a  random  time  variable  Xk  is  Ts(i+h/4), 
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Figure  5.6:  Semi-analytic  simulation  model  to  evaluate  BER  of  asynchronous  MC- 

CDMA 


g  =  1,2,  ...  ,8,  h  =  1,2,3,4,  and  i=0,l,2,  ...  ,U-1.  Then  the  received  signal  in  a  base  station 


IS 


r{m)  =  ^aj'{m)  =  '^Ojaj{m-T.),  where  j  are  active  users.  (5.19) 

The  demodulation  is  performed  for  a  user  number  k.  The  perfect  time  and  phase 
synchronization  and  ideal  detection  are  assumed  in  this  simulation.  The  BER  for  a  user  is 
different  for  the  spreading  sequence  assigned  to  the  user,  so  the  average  BER  should  take 
the  average  of  BER  of  all  possible  users.  Therefore,  the  average  BER  of  the 
asynchronous  BPSK  MC-CDMA  system  in  an  AWGN  channel  that  is  calculated  by  the 
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Eb/No  (dB) 

Figure  5.7:  Average  BER  performance  of  PSZ  and  Zadoff-Chu  sequences  (N=31,K=4,8) 
mean  of  Q-function  modified  form  (3.15). 

hE 

BER  =  mean[mean[Q{cl^ ,real{di, ,')  — -)]]  (5.20) 

where  k  are  all  possible  active  users,  time  index  t  =  1,2,3,  •••  ,10^,  d^t  are  the  original 
data  of  user  k  (BPSK  symbols,  ±1)  and  dk,t'  are  the  demodulated  data  of  user  k.  Thus, 
dk,treal(dk,t')  are  the  distances  to  the  decision  boundary  (imaginary  axis)  from  the 
demodulated  symbols.  It  gives  the  correct  BER  even  though  the  distance  is  negative. 
There  is  an  error  floor  on  BER  plot  when  the  distance  is  negative  because  the  mean  of  Q- 
functions  never  goes  to  zero  even  though  Eb/No  goes  to  infinity.  One  million  bits  of  data 
are  used  in  this  semi-analytical  simulation.  The  average  BER  performance  of  the  MC- 
CDMA  (N=31,  K=4,8)  with  Zadoff-Chu  sequences  and  that  of  the  MC-CDMA  with  the 
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Table  5.1:  PAPR  Suppression  Performance  in  asynchronous  MC-CDMA  system 


XT 

K 

RRC 
Roll-off 

Window 
Length 

Y 

Max 
PAPR 

Suppressed 
Max  PAPR 

Max 
Distance 

Power 

Loss 

31 

4 

0.00 

20Ts 

1.55 

6.14dB 

3.95dB 

0.11 

-.004dB 

31 

4 

0.10 

20Ts 

1.45 

5.39 

3.33 

0.12 

-.006 

31 

4 

0.20 

I2Ts 

1.35 

4.53 

2.71 

0.12 

-.004 

31 

4 

0.35 

12Ts 

1.25 

3.38 

2.27 

0.10 

-.005 

128 

16 

0.00 

20Ts 

1.55 

6.19 

3.88 

0.12 

-.004 

128 

16 

0.20 

12Ts 

1.35 

4.59 

2.64 

0.12 

-.009 

31 

4 

0.00 

20T, 

1.40 

6.14 

3.13 

0.16 

-.010 

31 

4 

0.20 

12Ts 

1.30 

4.53 

2.53 

0.16 

-.011 

PSZ  sequences  is  shown  in  Figure  5.7.  The  PSZ  sequences  are  designed  for  RRC  (roll-off 
0)  and  have  3dB  PAPR  suppression  capability.  There  is  no  power  efficiency  loss  caused 
by  using  the  PSZ  sequences.  It  is  interesting  that  the  BER  of  PSZ  is  a  few  hundredth  of 
dB  better  than  that  of  Zadoff-Chu.  It  is  hard  to  tell  that  the  PSZ  has  better  BER  because 
the  difference  is  so  small,  but  at  least  we  can  say  that  there  is  no  power  efficiency 
degradation  caused  by  the  PAPR  suppression  of  PSZ  in  asynchronous  MC-CDMA. 

Table  5.1  shows  the  semi-analytical  simulation  result  of  PAPR  suppression  in  an 
asynchronous  MC-CDMA  system.  N  is  the  processing  gain,  K  is  the  number  of  active 
users,  and  y  is  the  PAPR  suppression  cutoff  level  in  the  APS  algorithm.  The  suppression 
rate  P  in  the  APS  algorithm  is  all  0.8  but  not  so  important  here  because  the  number  of 
iterations  is  high  enough.  The  number  of  iterations  is  30  in  this  simulation.  The  max 
PAPR  means  that  the  PAPR  of  a  user  whose  PAPR  is  higher  than  any  other  users. 

The  max  distance  is  the  mean  distance  between  the  Zadoff-Chu  sequences  and  the 
PSZ  sequences  of  the  user  whose  distance  is  longer  than  any  other  users'.  This  indicates 
how  much  the  PSZ  sequences  are  different  from  the  Zadoff-Chu  sequences  when  each 
element  of  the  Zadoff-Chu  has  unit  power. 
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max  distance  =  max  [  mean  (Sj'^'^(n)-Sj™-^'"'(n))] ,  (5.21) 

where  j=l,2,  ...  ,U,  n=l,2,  ...  ,N.  Therefore,  when  y  is  lower,  the  suppressed  max  PAPR 
is  lower  and  the  max  distance  is  longer.  Generally,  this  longer  maximum  distance  causes 
higher  power  efficiency  loss,  but  it  does  not  in  this  simulation.  The  last  two  rows  of  the 
Table  5.1  are  the  cases  of  the  lower  y  and  the  lower  suppressed  PAPR  and  the  longer 
distance.  But  the  power  efficiency  loss  does  not  increase  as  expected.  It  means  the  new 
PSZ  sequences  have  good  BER  performance  in  an  asynchronous  MC-CDMA  system. 

The  power  efficiency  loss  is  positive  if  the  data  is  spread  by  the  PSZ  sequences 
and  despread  by  the  original  Zadoff-Chu  sequences  as  in  synchronous  MC-CDMA.  The 
power  loss  is  higher  when  the  max  distance  is  longer  in  the  case.  However,  the  power 
loss  is  still  pretty  small  (-.01  dB).  The  spectral  efficiency  is  better  and  the  original  PAPR 
is  higher  when  the  roll-off  factor  of  the  RRC  pulse  shaping  filter  is  lower.  The  original 
max  PAPR  of  RRC  roll-off  0  is  6.4dB  in  ideal  case  but  the  simulation  result  is  6.14dB 
when  20Ts  length  rectangular  windowed  RRC  roll-off  0.00001  pulse  shape  is  adopted. 
The  suppressed  PAPR  is  3.13dB  for  RRC  roll-off  0.0  and  2.53dB  for  roll-off  0.2  by  using 
PSZ  (PAPR  suppressed  Zadoff-Chu)  sequences,  while  the  original  PAPR  is  6.14dB  for 
RRC  roll-off  0.0  and  4.53dB  for  roll-off  0.2  by  using  Zadoff-Chu  sequences. 

The  power  efficiency  loss  in  a  Rayleigh  fading  channel  (f-loss)  is  also  a  negative 
value  but  the  absolute  value  is  smaller  than  that  in  AWGN  channel  (g-loss).  The  f-loss  is 
-.005  when  the  g-loss  is  -.010  for  roll-off  0.0  and  the  f-loss  is  -.004  when  the  g-loss  is  - 
.011  for  roll-off  0.2.  The  absolute  value  is  smaller  because  the  amplitude  adjustment 
process  is  less  sensitive  in  a  Rayleigh  fading  channel.  However,  all  of  these  power  loss 
values  are  approximately  zero  and  we  can  say  there  is  no  power  efficiency  loss. 
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The  characteristics  of  PSZ  are  numerically  searched  as  follows.  Let's  call  the  PSZ 
sequences  that  suppresses  the  PAPR  from  6.14dB  to  3.95dB  as  PSZl  and  the  PSZ 
sequences  that  suppresses  the  PAPR  from  6.14dB  to  3.13dB  as  PSZ2.  The  mean  discrete 
autocorrelation  of  the  IDFT  sequences  of  PSZl  (9c.c(n))  is  0.0049  for  n^^O  and  that  of 
PSZ2  is  0.01 18,  while  that  of  Zadoff-Chu  is  0.  If  it  compare  to  the  mean  0s,s(n)  of  Gold 
sequence  (0.1378),  the  autocorrelation  property  of  the  IDFT  sequences  of  PSZ  is  very 
good.  The  mean  of  the  discrete  crosscorrelation  of  IDFT  sequences  of  PSZl  {6^.^^..{n))  is 

0.1794  and  the  maximum  of  them  is  0.1796  for  \^].  The  mean  of  the  discrete 
crosscorrelation  of  IDFT  sequences  of  PSZ2  is  0. 1 786  and  the  maximum  of  them  is 
0.1802  for  x^].  These  are  similar  or  slightly  better  than  those  of  Zadoff-Chu  (all  0.1796) 
and  this  result  is  consistent  with  the  semi-analytical  simulation  result  of  BER 
performance  in  that  the  BER  of  PSZ  is  slightly  better  than  that  of  Zadoff-Chu.  Therefore, 
the  PSZ  sequences  have  the  same  total  power,  the  low  autocorrelation  of  IDFT 
sequences,  and  the  low  crosscorrelation  of  IDFT  sequences  as  Zadoff-Chu  sequences 
have.  Moreover,  the  PAPR  of  PSZ  sequences  is  2-3dB  lower  than  that  of  Zadoff-Chu 
sequences  in  general  asynchronous  MC-CDMA  systems  and  there  is  no  virtual  cost  for 
the  improvement. 

5.5  Summary 

The  PAPR  of  asynchronous  MC-CDMA  depends  on  the  spreading  sequence. 
Povopic  [43]  suggested  a  Zadoff-Chu  sequence  as  the  optimal  spreading  sequence  for 
asynchronous  MC-CDMA  system  taking  PAPR,  dynamic  range,  and  BER  performance 
into  account.  However,  the  PAPR  of  the  Zadoff-Chu  sequence  is  still  more  than  6dB  in 
terms  of  continuous  PAPR.  Therefore,  a  PSZ  sequence  is  developed  by  using  the  APS 
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algorithm  to  suppress  the  PAPR  more.  The  characteristic  of  PSZ  sequences  is  similar  to 
Zadoff-Chu  sequences.  The  PAPR  of  PSZ  is  2-3dB  lower  than  that  of  Zadoff-Chu  and 
there  is  no  power  efficiency  loss  in  BER.  Thus,  a  new  low  PAPR  generating  spreading 
sequence  for  asynchronous  MC-CDMA  is  introduced. 


CHAPTER  6 
CONCLUSIONS  AND  FUTURE  WORK 

In  this  chapter,  the  key  contributions  of  this  dissertation  to  the  existing  body  of 
research  will  be  summarized.  Possible  areas  of  future  research  will  also  be  suggested. 

6.1  Conclusions 

In  this  dissertation,  PAPR  suppression  schemes  for  tc/4  QPSK,  OFDM,  and  MC- 
CMDA  are  investigated  and  developed.  The  goal  of  PAPR  suppression  is  to  increase  the 
efficiency  of  linear  power  amplifiers,  thus,  to  increase  the  power  efficiency  of  a  system. 
PAPR  suppression  is  one  of  the  efforts  to  find  a  bandwidth  and  power  efficient 
modulation.  Therefore,  PAPR  suppression  should  not  sacrifice  the  bandwidth  efficiency. 

n/4  QPSK  is  the  one  of  the  best  bandwidth  and  power  efficient  linear  modulation. 
To  increase  power  efficiency  of  n/4  QPSK,  the  APS  algorithm  is  adopted.  The  algorithm 
suppresses  the  PAPR  of  ti/4  QPSK  with  RRC  (roll-off  0.35)  from  3.2dB  to  1.7dB.  The 
computational  complexity  is  reduced  from  23Mips  to  3Mips  at  cost  of  0.2dB  less  PAPR 
suppression.  The  look-up  table  can  substitute  for  the  APS  algorithm  to  remove  high 
computafional  complexity  at  the  cost  of  a  few  tenth  of  dB  less  PAPR  suppression.  The 
hybrid  method  of  the  look-up  table  and  the  ASP  algorithm  suppresses  the  PAPR  more 
than  the  look-up  table  does  at  cost  of  some  complexity. 

OFDM  is  a  very  effective  and  popular  modulation  for  high  data  rate  transmission 
but  has  very  high  PAPR.  The  discrete  PAPR  of  OFDM  caused  by  the  IDFT  procedure 
can  be  suppressed  by  the  amplitude  limiter  and  the  addifional  PAPR  caused  by  pulse 
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shaping  can  be  suppressed  by  the  APS  algorithm.  Thus,  the  Hmiter  plus  APS  scheme  (4 
iterations)  suppresses  the  continuous  PAPR  of  128-QPSK-OFDM  with  RRC  (roll-off 
0.35)  from  12.7dB  to  5.2dB  without  any  bandwidth  efficiency  loss  when  IdB  average 
power  loss  is  allowed  at  10"^  SER.  The  scheme  is  better  than  any  other  schemes  in  BPSK- 
OFDM  and  QPSK-OFDM  when  the  number  of  sub-carriers  is  large  (>32).  The  average 
power  loss  in  a  Rayleigh  fading  channel  is  about  one  third  of  the  average  power  loss  in  an 
AWGN  channel. 

The  PAPR  property  of  Synchronous  MC-CDMA  with  GC  sequences  is  similar  to 
that  of  OFDM.  Therefore,  the  limiter  plus  APS  scheme  (4  iterations)  can  suppress  the 
PAPR  of  synchronous  MC-CDMA  (N=128,  K=32)  with  RRC  (roll-off  0.35)  from  12. IdB 
to  3.9dB  without  any  bandwidth  efficiency  loss  when  IdB  average  power  loss  is  allowed 
at  a  10"^  SER.  The  average  power  loss  in  a  Rayleigh  fading  channel  is  about  one  fourth  of 
the  average  power  loss  in  an  AWGN  channel. 

The  PAPR  of  Asynchronous  MC-CDMA  is  determined  by  the  set  of  spreading 
sequences  adopted.  Zadoff-Chu  sequences  (N=31)  generate  OdB  discrete  PAPR  and 
6.4dB  continuous  PAPR  signal,  while  gold  sequences  (N=31)  generate  about  9dB 
continuous  PAPR.  The  continuous  PAPR  of  Zadoff-Chu  for  RRC  (roll-off  0.2)  is  4.5dB. 
A  set  of  PSZ  sequences  is  found  by  using  the  APS  algorithm  to  suppress  the  PAPR.  The 
continuous  PAPR  of  PSZ  is  3. IdB  for  RRC  (roll-off  0.0)  and  2.5dB  for  RRC  (roll-off 
0.2).  It  should  be  noted  that  there  is  no  cost  in  terms  of  power  efficiency,  bandwidth 
efficiency,  or  computational  complexity  to  use  the  PSZ  sequence. 
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The  amplitude  limiter  and  the  APS  algorithm  and  their  application  suppress  the 
PAPR  of  a  signal  without  any  bandwidth  expansion  by  amplitude  adjustment  before  pulse 
shaping.  They  perform  better  than  any  other  previous  schemes  in  most  cases. 

6.2  Future  Work 

The  PAPR  suppression  schemes  for  asynchronous  MC-CDMA  can  be  applied  in 
general  asynchronous  DS-CDMA.  The  discrete  PAPR  of  asynchronous  DS-CDMA  with 
gold  spreading  sequences  (N=31)  is  known  as  OdB  but  the  continuous  PAPR  for  RRC 
(roll-off  0.2)  is  3.7-5.7dB.  Therefore,  it  is  desired  to  suppress  the  PAPR  even  in  DS- 
CDMA.  The  APS  algorithm  may  be  applied  to  gold  spreading  sequences  as  it  applied  to 
Zadoff-Chu  sequences  in  Chapter  5.  Then,  a  new  PAPR  suppressed  gold  (PSG)  sequence 
may  be  found. 

Other  than  this,  the  amplitude  limiter  and  the  APS  algorithm  can  be  applied  to  any 
high  PAPR  generating  modulation  that  has  D/A  pulse  shaping  process.  Currently  popular 
asymmetric  digital  subscriber  line  (ADSL)  adopts  a  modulation  format  called  discrete 
multi-tone  (DMT).  DMT  has  a  similar  structure  to  OFDM  and  a  high  PAPR  problem,  so 
it  might  be  another  good  candidate  for  these  schemes. 


124 


REFERENCES 

[I]  S.  L.  Miller,  R.  J.  O'Dea,  "Peak  Power  and  Bandwidth  Efficient  Linear  Modulation," 
IEEE  Transactions  on  Communications,  Vol.  46,  No.  12,  Dec.  1998,  pp.  1639-1648 

[2]  G.  J.  Foschini,  R.  D.  Gitlin,  S.  B.  Weinstein,  "On  the  Selection  of  a  Two- 
Dimensional  Signal  Constellation  in  the  Presence  of  Phase  Jitter  and  Gausian  Noise," 
Bell  System  Technical  Journal,  Vol.  52,  No.  6,  Jul.-Aug.,  1973,  pp.  927-965 

[3]  M.  k.  Simon,  J.  G.  Smith,  "Hexagonal  Multiple  Phase-and-Amplitude-Shift-Keyed 
Signal  Set,"  IEEE  Transactions  on  Communications,  Vol.  COM-21,  No  10,  Oct.  1973, 
pp.  1108-1115 

[4]  D.  H.  Morals,  K.  Feher,  "The  Effects  of  Filtering  and  Limiting  on  the  Performance  of 
QPSK,  Offset  QPSK,  and  MSK  Systems,"  IEEE  Transactions  on  Communications,  Vol. 
COM-28,  No.  12,  Dec.  1980,  pp.  1999-2009 

[5]  M.  U.  Chang,  M.  C.  Austin,  "Prefiltered  Quadrature  Phase-Shift  Keying:  A 
Bandwidth-and-Power-Efficient  Modulation,"  Journal  on  Spacecraft  and  Rockets,  Vol. 
18,  No.  3,  May-Jun.  1981,  pp.279-284 

[6]  M.  C.  Austin,  M.  U.  Chang,  "Quadrature  Overlapped  Raised-Cosine  Modulation," 
IEEE  Transactions  on  Communications,  Vol.  COM-29,  No.  3,  Mar.  1981,  pp.  237-249 

[7]  T.  Le-Ngoc,  K.  Feher,  H.  P.  Van,  "New  Modulation  Techniques  for  Low-Cost  Power 
and  Bandwidth  Efficient  Satellite  Earth  Stations,"  IEEE  Transactions  on 
Communications,  Vol.  COM-30,  No.  1,  Jan.  1982,  pp.  275-283 

[8]  J-S  Seo,  K.  Feher,  "SQAM:  A  New  Superimposed  QAM  Modem  Technique,"  IEEE 
Transactions  on  Communications,  Vol.  COM-33,  No.  3,  Mar.  1985,  pp.  296-300 

[9]  F.  Amoroso,  "The  Bandwidth  of  Digital  Data  Modulation,"  IEEE  Communication 
Magazine,  Vol.  18,  Nov.  1980,  pp.  13-24 

[10]  F.  deJager,  C.  B.  Dekker,  "Tamed  Frequency  Modulation,  A  Novel  Method  to 
Achieve  Spectrum  Economy  in  Digital  Transmission,"  IEEE  Transactions  on 
Communications.  Vol.  COM-26,  No.  5,  May  1978,  pp.  534-542 

[II]  Y.  Akaiwa,  Y.  Nagata,  "Highly  Efficient  Digital  Mobile  Communications  with  a 
Linear  Modulation  Method,"  IEEE  Journal  on  Selected  Areas  in  Communications,  Vol. 
SAC-5,  No.  5,  Jun.  1987,  pp.  890-895 


125 


[12]  C-L  Liu,  K.  Feher,  "Noncoherent  Detection  of  tiM-QPSK  Systems  in  a  CCI-AWGN 
Combined  Interference  Environment,"  IEEE  Vehicular  Technology  Conference,  Vol.  I, 
1989,  pp. 83-94 

[13]  Y.  Yamao,  S.  Saito,  H.  Suzuki,  T.  Nojima,  "Performance  of  7r/4  QPSK  Transmission 
for  Digital  Mobile  Radio  Applications,"  IEEE  Global  Telecommunications  Conference 
and  Exhibition,  Vol.  1,  1989,  pp.443-447 

[14]  S.  L.  Miller,  R.  J.  O'Dea,  'Adaptive  Peak  Suppression  for  Power  and  Bandwith 
Efficient  Linear  Modulation',  IEEE  Global  Telecommunication  Conference  (Mini- 
Conference),  Nov.  1997,  pp.  1-5. 

[15]  P.  A.  Barker,  "Phase-Modulation  Data  Sets  for  Serial  Transmission  at  2,000  and 
2,400  Bits  per  Second,"  AIEE  Transactions  on  Communications  Electronics,  Jul.  1962, 
pp.  166-171 

[16]  S.  Ariyavisitakul,  T-P  Liu,  "Characterizing  the  Effects  of  Nonlinear  Amplifiers  on 
Linear  Modulation  for  Digital  Portable  Radio  Communications,"  IEEE  Global 
Telecommunications  Conference  and  Exhibition,  Vol.  1,  1989,  pp. 448-453 

[17]  Antonio  S.  Fedon,  "University  Research  in  Peak  Power  and  Bandwidth  Efficient 
Modulation,"  University  Research  Report,  University  of  Florida,  1996. 

[18]  B.  R.  Saltzberg,  "Performance  of  an  Efficient  Parallel  Data  Transmission  System," 
IEEE  Transactions  on  Communication  Technology,  Vol.  COM-15,  No.  6,  Dec.  1967,  pp. 
805-811 

[19]  S.  B.  Weinstein,  P.  M.  Ebert,  "Data  Transmission  by  Frequency-Division 
mutiplexing  Using  the  Discrete  Fourier  Transform,"  IEEE  Transactions  on 
Communication  Technology,  Vol.  COM-19,  No.  5,  Oct.  1971,  pp.  628-634 

[20]  W.  Y.  Zou,  Y.  Wu,  "COFDM:  an  Overview,"  IEEE  Transactions  on  Broadcasting, 
Mar.  1995,  Vol.  41,  No.  1,  pp.  1-8. 

[21]  T.  A.  Wilkinson,  A.  E.  Jones,  "Minimization  of  Peak-to-Mean  Envelope  Power 
Ratio  of  Multicarrier  Transmission  Schemes  by  Block  Coding,"  IEEE  Vehicular 
Technology  Conference,  July  1995,  Vol.  2,  pp.  825-829. 

[22]  R.  D.  J.  Van  Nee,  "OFDM  Codes  for  Peak -to-Average  Power  Reduction  and  Error 
Correction,"  IEEE  Global  Telecommunications  Conference,  Nov.  1996,  Vol.  1,  pp.  740- 
744. 

[23]  M.  Friese,  "OFDM  Signal  with  Low  Crest-Factor,"  IEEE  Global 
Telecommunications  Conference,  Nov.  1997,  Vol.1,  pp.290-294. 

[24]  R.  W.  Bauml,  R.  F.  H.  Fischer,  J.  B.  Huber,  "Reducing  the  Peak-to-Average  Power 
Ratio  of  Multicarrier  Modulation  by  Selected  Mapping,"  Electronic  Letters,  Oct.  1996, 
Vol.  32,  No.  22,  pp.  2056-2057. 


126 


[25]  S.  H.  Muller,  J.  B.  Huber,  "OFDM  with  Reduced  Peak-to-Average  Power  Ratio  by 
Optimum  Combination  of  partial  Transmit  Sequences,"  Electronic  Letters,  Feb.  1997, 
Vol.  33,  No.  5,  pp.  368-369. 

[26]  S.  H.  Muller,  J.  B.  Huber,  "A  Comparison  of  Peak  Power  Reduction  Schemes  for 
OFDM,"  IEEE  Global  Telecommunication  Conference,  Nov.  1997,  Vol.  1,  pp.  1-5. 

[27]  C.  Tellambura,  "Phase  optimization  criterion  for  reducing  peak-to-average  power 
ratio  in  OFDM,"  Electronic  Letters,  Jan.  1998,  Vol.  34,  No.  2,  pp.  169-170. 

[28]  X.  Li,  L.  Cimini,  "  Effects  of  Clipping  and  Filtering  on  the  Performance  of  OFDM," 
IEEE  Communications  Letters,  May.  1998,  Vol.  2,  No.  5,  pp.  131-133. 

[29]  R.  Van  Nee,  A.  Wild,  "Reducing  the  Peak-to-Average  Power  Ratio  of  OFDM"  IEEE 
Vehicular  Technology  Conference,  1998,  Vol.  3,  pp.  2072-2076 

[30]  R.  O'Neill,  L.  B.  Lopes,  "Performance  of  Amplitude  Limited  Multitone  Signals," 
IEEE  Vehicular  Technology  Conference,  1994,  pp.  1675-1679 

[31]  N.  Dinur,  D.  Wulich,  "Peak  to  Average  Power  Ratio  in  Amplitude  Clipped  High 
Order  OFDM,"  IEEE  Military  Communications  Conference  MILCOM,  1998,  Vol.  2,  pp. 
684-687. 

[32]  D.  Wulich,  L.  Goldfeld,  "Reduction  of  Peak  Factor  in  Orthogonal  Multicarrier 
Modulation  by  Amplitude  Limiting  and  Coding,"  IEEE  Transactions  on 
Communications,  Vol.  47,  No.  1,  Jan.  1999,  pp.  18-21 

[33]  L.  W.  Couch  II,  "Digital  and  Analog  Communication  System,"  Prentice  Hall,  1997, 
Fifth  Edition,  pp.  18,  160. 

[34]  G.  Fettweis,  A.  Bahai,  K.  Anvari,  "On  Multi-Carrier  Code  Division  Multiple  Access 
(MC-CDMA)  Modem  Design,"  IEEE  Vehicular  Technology  Conference,  Vol.  3,  Jun. 
1994,  pp.  1670-1674 

[35]  E.  Sourour,  M.  Nakagawa,  "Performance  of  Orthogonal  Multicarrier  CDMA  in  a 
Multipath  Fading  Channel,"  IEEE  Transaction  on  Communications,  Vol.  44,  No.  3,  Mar. 
1996,  pp.  356-367 

[36]  S.  Hara,  R.  Prasad,  "Overview  of  Multicarrier  CDMA,"  IEEE  Communication 
Magazine,  Dec.1997,  pp.  126-133 

[37]  A.  Bahai,  B.  Saltzberg.  "Multi-Carrier  Digital  Communications,"  Kluwer  Academic 
/  Plenum  Publishers,  New  York,  1999,  pp.  210-213 

[38]  V.  Aue,  G.  Fettweis,  "Multi-Carrier  Spread  Spectrum  Modulation  with  Reduced 
Dynamic  Range,"  IEEE  Vehicular  Technology  Conference,  Apr.  1996,  pp.  914-917 


127 


[39]  M.  Schroeder,  "  Synthesis  of  Low-Peak-Factor  Signals  and  Binary  Sequences  with 
Low  Autocorreation,"  IEEE  Transaction  on  Information  Theory,  Jan.  1970,  pp.  85-89 

[40]  S.  Boyd,  "  Multitone  Signals  with  Low  Crest  Factor,"  IEEE  Transaction  on  Circuits 
and  Systems,  Vol.  CAS-33,  Oct.  1986,  pp.  1018-1022 

[41]  E.  Van,  J.  Schoukens,  J.  Renneboog,  "Peak  Factor  Minimization  of  Input  and  Output 
Signals  of  Linear  Systems,"  IEEE  Transaction  on  Instrumentation  and  Measurement, 
Vol.  37,  No.  2,  Jun.  1988,  pp.  207-212 

[42]  D.  Gimlin,  C.  Patisaul,  "On  minimizing  the  Peak-to-Average  Power  Ratio  for  the 
Sum  of  N  sinusoids,"  IEEE  Transaction  on  Communications,  Vol.  41,  No.  4,  Apr.  1993, 
pp.  631-635 

[43]  B.  Popovic,  "Spreading  Sequences  for  Multi-Carrier  CDMA  System,"  lEE 
Colloquium  (Digest),  No.  129,  May  1997,  pp.  8/1-8/6 

[44]  M.  Golay,  "Complementary  Series,"  IRE  Transactions  on  Information  Theory,  Vol. 
7,  Apr.  1961,  pp.  82-87 

[45]  B.  Popovic,  "Spreading  Sequences  for  Multicarrier  CDMA  Systems,"  IEEE 
Transaction  on  Communications,  Vol.  47,  No.  6,  Jun.  1999,  pp.  918-926 

[46]  L.  Freiberg,  A.  Annamalai,  V.  Bhargava,  "Crest  Factor  Reduction  Using  Orthogonal 
Spreading  Codes  in  Multi-Carrier  CDMA  systems,"  IEEE  International  Symposium  on 
Personal,  Indoor  and  Mobile  Radio  Communications,  Vol.  1,  1997,  pp.  120-124 

[47]  H.  Ochiai,  H.  Imai,  "OFDM-CDMA  with  Peak  Power  Reduction  Based  on 
Spreading  Sequences,"  IEEE  International  Conference  on  Communications,  Vol.  3,  1998, 
pp.  1299-1303 

[48]  T.  Ho,  V.  Wei,  "Synthesis  of  Low-Crest  Waveforms  for  Multicarrier  CDMA 
System,"  IEEE  Communication  Theory  Mini-Conference,  1995,  pp.  131-135 

[49]  M.  Popovic,  "Synthesis  of  power  Efficient  Multitone  Signals  with  Flat  Amplitude 
Spectrum,"  IEEE  Transaction  on  Communications,  Vol.  39,  No.  7,  Jul.  1991,  pp.  1031- 
1033 

[50]  D.  Chu,  "Polyphase  Codes  with  Good  Periodic  Correlation  Properties,"  IEEE 
Transaction  on  Information  Theory,  Vol.  IT- 18,  Jul.  1972,  pp.  531-532 

[51]  M.  Popovic,  "Generalized  Chirp-Like  Polyphase  sequences  with  Optimum 
Correlation  Properties,"  IEEE  Transaction  on  Information  Theory,  Vol.  38,  No.  4,  Jul. 
1992,  pp.  1406-1409 

[52]  M.  Pursley,  "Performance  Evaluation  for  Phase-Coded  Spread-Spectrum  Multiple- 
Access  Communication  -  Part  I:  System  Analysis,"  IEEE  Transactions  on 
Communications,  Vol.  25,  No.  8,  Aug.  1977,  pp.  795-799 


BIOGRAPHICAL  SKETCH 
Young-Seo  Park  was  bom  in  Seoul,  Korea,  in  1969.  He  graduated  from  Hyundai 
High  School,  Seoul,  Korea,  in  February  1988.  In  1990,  he  served  in  the  Korean  army.  He 
received  the  B.S.  degree  in  electrical  engineering  in  February  1994  from  Hanyang 
University,  Seoul,  Korea.  In  June  1995,  he  came  to  the  United  States  to  study  electrical 
engineering  at  the  University  of  Florida.  He  received  M.S.  degree  in  electrical  and 
computer  engineering  in  December  1 997  from  the  University  of  Florida.  In  August  2000, 
he  received  his  Ph.D.  degree  in  electrical  and  computer  engineering  from  the  University 
of  Florida. 


128 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fuUy  adequate,  in  scope  and  quality, 
as  a  dissertation  for  the  degree  of  Doctor  of  Philosophy. 


Scott  L.  Miller,  Chair 

Associate  Professor  of  Electrical  and 

Computer  Engineering 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope  and  quality, 
as  a  dissertation  for  the  degree  of  Doctor  of  Philosophy 


Lecm  W.  Couch  H,  Co-chair 
Professor  of  Electrical  and  Computer 
Engineering 

I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope  and  quaUty, 
as  a  dissertation  for  the  degree  of  Doctor  of  Philosophy 


Tan^Weng 
Assistant  Professor  of  Electrical  and 
Computer  Engineering 


I  certify  that  I  have  read  this  study  and  that  in  my  opinion  it  conforms  to 
acceptable  standards  of  scholarly  presentation  and  is  fully  adequate,  in  scope  and  quality, 
as  a  dissertation  for  the  degree  of  Doctor  of  Philosophy. 


•1/ 


Randy  Carter/ 
Professor  of  Statistics 

This  dissertation  was  submitted  to  the  Graduate  Faculty  of  the  College  of 
Agricultural  and  Life  Sciences  and  to  the  Graduate  School  and  was  accepted  as  partial 
fulfillment  of  the  requirements  for  the  degree  of  Doctbj^xjfPhiios^Aj 

August  2000  /(/    ""^--^f  [y 

M  J.  6hanian 

Dean,  College  of  Engineering 


Winfred  M.  PhUlips 
Dean,  Graduate  School 


